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ABSTRACT 



An acoustic resonance based technique using a free-free bar has been extended to 
investigate the complex (storage and loss) moduli of non-magnetic materials having circular 
cross section. Using this technique, the bar can be selectively excited in three independent 
vibrational modes, i.e., torsional, flexural, and longitudinal modes. The torsional mode 
yields the shear modulus. Either the flexural or longitudinal mode can be used to obtain 
Young's modulus. These resonant modes can be tracked continuously by means of a 
phase-locked-loop (PLL) as the temperature (and resonant frequency) of the rod is 
changed. The in-phase amplitude of the receiver output of the electrodynamic transducer is 
proportional to the quality factor, Q, of the material. It can be used to continuously track 
the loss tangent ( = HQ ) of the material as a function of temperature and frequency. 
Results for complex shear modulus and Young's modulus were obtained for a castable 
epoxy type PR 1592 and complex shear modulus for polymethyl methacrylate (PMMA) and 
Uralite 3130. Over the temperature and frequency range that was accessible, a clear 
viscoelastic transition was observed in both the storage modulus and loss tangent curves of 
PR 1592. 
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THESIS DISCLAIMER 

The reader is cautioned that computer programs developed in this research may not 
have been exercised for all cases of interest. While every effort has been made, within the 
time available, to ensure that the programs are free of computational and logic errors, they 
cannot be considered validated. Any application of these programs without additional 
verification is at the risk of the user. 
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I. INTRODUCTION 



A. MOTIVATION 

It is important that the mechanical properties of elastomeric or rubber-like materials be 
accurately measured and their dependence on temperature, static pressure and other 
ambient parameters be determined for many fields of science and engineering. These 
rubber-like materials are being used in various applications such as materials for anti- 
vibration mounting, hydrophone designs and other acoustics applications. It is known that 
the static and dynamic moduli of plastics can differ substantially [Ref. 1] and in many 
situations, it is the dynamic moduli that are the appropriate moduli which determine the 
quantity of interest such as the dissipative characteristics of vibration isolation material, the 
resonant frequencies of a linear mechanical system and the sensitivity of many transduction 
mechanisms. 

For most applications, the elastic moduli at the frequencies of intended operations are 
of interest rather than the static modulus. Nevertheless, manufacturers' specifications for 
elastic constants of castable polymers are not particularly useful as they are usually 
determined by static techniques and rarely contain more than one modulus if available. For 
an isotropic, homogeneous material, all the elastic properties can be completely 
characterized by just two moduli. Therefore, development of a technique that can 
accurately determine these dynamic properties is essential for applications engineers and 
designers and is the topic of this research. 

B. TECHNIQUE 

In general, the dynamic moduli of a material can be determined by forced vibration test 
methods. That is, directly measuring the force and resulting displacement or acceleration 
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and using Hooke's law. These can be broadly classified as resonant and non-resonant 
methods, each method having its advantages and disadvantages [Ref. 2]. Resonant based 
techniques have the advantage of higher signal-to-noise ratio because at resonance, the 
response of the sample is quality factor, Q, times higher than the response off resonance. 
This is an important consideration with a high loss material. Further, as the fundamental 
measurement is frequency, this suggests that one can obtain extremely high precision with 
a relatively inexpensive instrument such as a frequency counter. The resonant technique 
used in this research has been described and discussed in detail in a paper written by 
Garrett [Ref. 3], using a transducer consisting of coUs of magnet wire placed in the 
magnetic field created by a pair of permanent magnets. Using identical transducers for 
driver and receiver (pick-up), it is able to selectively excite a single rod-shaped sample of 
circular or elliptical cross-section into torsional, flexural and longitudinal resonant modes. 
After obtaining the resonant frequency, knowing the mass density, p , and dimensions of 
the sample, the storage modulus can be obtained from the fact that the bar resonance is 
proportional to the appropriate wave speed. The complex moduli of the rubber-like 
material can be obtained together with the measurement of the quality factor, Q , or the free 
decay time, T . A commercially available device can be obtained from Reference 4. 

C. HISTORY 

The resonant bar technique, as it will be called from now onward, is a refinement of 
one developed by Barone and Giacomini [Ref. 5] to study the modes of vibration of bars of 
different cross sections. It was used as a teaching laboratory experiment at University of 
California Los Angeles by Professor Isadore Rudnick and currently, it is being used in 
acoustic laboratory courses at Naval Postgraduate School to demonstrate modes of bar, and 
recently, by Wetterskog, Beaton and Serocki [Ref. 6], to determine the dynamic moduli 
and their temperature dependence for numerous sample of materials for fiber-optic 
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hydrophone applications. Improvements to the resonant bar technique and continuous 
resonance tracking to yield both storage and loss moduli as a function of temperature and 
frequency are described in this thesis. 
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II. THEORY 



A. DYNAMIC MODULI DETERMINATION 

A uniform, isotropic, cylindrical rod of a homogeneous solid with a diameter, d, and 
length L, with L> X» d, so that radial motion can be neglected, will propagate three 

independent waves. These modes will exhibit resonances at appropriate frequencies, 
depending on the boundary conditions imposed on the ends of the rod. Applying the 
boundary conditions to a slender rod with unrestricted ends, i.e., free-free boundary 
conditions, the solutions to the resonant modes can be obtained [ Ref. 7: pp. 57-76, Ref. 
3]. 

1 . Non-dispersive Modes 

The displacement, <^, and angular displacement, associated with the 
longitudinal and torsional modes satisfy an ordinary second-order wave equation [Ref. 8 : 
pp.94, 112]. For a free-free boundary condition, the resonances are harmonically related 
and correspond to an integral number of half-wavelength contained within the length of the 
rod. It should be noted that the assumed boundary condition does not take the added mass 
of the transducers and their adhesive into account. Nevertheless, the effect is not 
ordinarily significant since the additional mass is rarely more than a few percent of the mass 
of the rod. 

a. Longitudinal Vibration - Dynamic Young's Modulus 

The phase speed of the longitudinal waves, , is given by Young's 

modulus, E, and the mass density, p [Ref. 7 ; p. 59] 
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The normal-mode frequencies of the bar are then given by 




( 2 ) 



where n is the mode number of the vibration corresponding to the number of nodes in the 
standing wave. When the wavelength, A , begins to become comparable with the lateral 
dimension of the bar, the assumption of slenderness fails and the normal-mode frequencies 
deviate increasingly from a harmonic progression. Nevertheless, for a slender rod, these 
two equations can be used to solve for Young's modulus of the bar 



b . Torsional Vibration • Dynamic Shear Modulus 

For a cylindrical rod, the phase speed for the torsional waves, cj, is 
given by the shear modulus, G , and the mass density, p [Ref. 8 : p.l 1] 



If the rod was elliptical, with major and minor radii, a and b, respectively, the speed is 
modified by the multiplicative factor [ lab /{a^ + b^ )] [Ref. 9]. Applying the free-free 
boundary condition leads to a series of harmonic modes with frequencies given by 




(3) 




(4) 




( 5 ) 



These two equations can be used to solve for the shear modulus of the bar 



( 6 ) 



T ^ 

G = 4pLH^) 

where n is again the mode number of the vibration. By obtaining the resonance frequency 
of both the torsional and longitudinal modes of a homogeneous, isotropic rod of known 
mass and dimensions, the complete set of elastic constants ( i.e., bulk modulus, Poisson’s 
ratio, etc.) can be determined as only two independent moduli are required to completely 
characterize such a system. 

2 . Dispersive Mode - Flexural Mode 

The flexural phase speed, Cf, is given as [Ref. 7 : p. 71] 

CF=n7lfKCL, (7) 

where /cis the radius of gyration given by [Ref. 7 : p.69] 

K'^ = (^)| dS 

' , ( 8 ) 
where S is the cross-sectional area of the rod, and z is the distance of an element above the 
neutral axis in the direction of flexure. Therefore the flexural wave phase speed varies with 
the square root of the frequency. A slender rod thus exhibits dispersion for flexural mode. 
The application of free-free boundary condition gives a series of modes that are not 
hannonic. The frequency of the nth overtone, tn , is given by [Ref. 7 : p. 75] 

= 3 0112,4.9994,7,9, 11... 

8L2 , (9) 

for a rod of circular cross-section, K = d 1 4, where d is the diameter of the rod. Solving 

the above equations to obtain Young's modulus 

E ^ . /I = 3.01 12, 4.9994, 7... 

. ( 10 ) 
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B . VISCOELASTIC MODEL OF SOLID 
I . Viscoelasticity 

Many solids exhibit primarily elastic effects when subjected to low levels of 
strain and obey Hooke's law. Under a low amplitude dynamic force, there is a 
corresponding deformation such that in the linear limit, the resulting strain is proportional 
to the magnitude and in-phase with the applied stress. The imparted energy is recoverable 
and not dissipated as heat. The ratio of the applied stress to the resulting normalized 
deformation or strain is equal to the elastic modulus. The modulus of these materials may 
be independent of frequency over a large range of frequencies. On the other hand, many 
liquids show appreciable viscous effects. The stress and strain are always 90° out-of-phase 
under infinitesimal rates of strain, and all of the shear energy transferred to the liquid is 
dissipated as heat. If the properties of a material fall between an ideal Hookean solid and 
an ideal Newtonian fluid, which is the case for most polymers, when a dynamic stress is 
applied, some of the energy input wiU be stored and some of the energy input will be 
dissipated. The material may recover part of its deformation when the stress is removed. 
Under sinusoidal oscillating stress, the strain is neither exactly in-phase with the applied 
stress nor 90° out-of-phase but it is somewhere in between. Materials whose behavior 
show such characteristics are called viscoelastic. 

. 2 . Simple Mechanical Model 

The simplest mechanical model of a viscoelastic system is one spring ( elastic ) 
and one dashpot ( viscous ), either in series or in parallel [Ref. 10 : pp. 16-18]. A series 
element is called a Maxwell element and a parallel element is called a Voigt element. 
Nevertheless, a combination of these elements are needed to accurately represent the 
viscoelastic properties of a material. 

For a Voigt element, the force, F, can be written in terms of the viscous element, 
R, the elastic element, k, and the displacement, as 
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F = {k + jcoR ) 

where,; = V-1 , hence, the equivalent complex modulus can be represented as 

k* =k ■¥ j(oR ^ 



( 11 ) 



( 12 ) 



which has a real part and imaginary part and a relaxation time constant, T , and a phase 
angle, 5 , given by 




(13) 



S =tan-^ 

k . 



(14) 



In this simple model, the real component is independent of frequency, but it may be a 
function of temperature. The imaginary part is definitely a function of frequency and it 
increases without bound with increasing frequency. Obviously, this model is inadequate 
to accurately describe the behavior of a viscoelastic material under dynamic stress. 

By adding an elastic element in parallel with a Maxwell model, as in Figure 1(a), 
which has a equivalent circuit as shown in Figure 1(b), the equivalent impedance can be 
expressed as 

Z = — = + ^2^2 

“ ytu k2-\-j(oR2^ (15) 

where u is the velocity equal to d^dt . Then the dynamic force can be written as 



F =Zu 




+ 



k2 +j(OR2 






= {ki 



aP-Rlk2 

+ ^ + J 

kl + ( 0 >/? 2)2 



kl + (a>/?2)^ 



)l 



(16) 



Hence the equivalent complex modulus is 

, * , oP'R^k2 . coR2kl 

k = ki+ — + j — 

kl + icoR2f kl + (o)/?2)^ . (17) 

From this expression, it can be seen that the real part and the imaginary part depend on 
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frequency, and if /?2 / ^2 is a function of temperature, k depends on temperature also. 
The frequency dependence of k' (real part of fc*) and k" (imaginary part of /:*) is shown in 
Figure 2. Other similar models can be found in Bland [Ref. 11 : pp.l 14-115]. 




k 

2 

jw 



R 

2 

AA 1 

(b) 

Figure 1. The simple viscoelastic model: 

(a) . The mechanical equivalent circuit. 

(b) . The electrical equivalent circuit. 



k 

1 



jw 



This model, with a single relaxation time constant, is similar in approach to that 
applied by Rudnick in describing the low temperature liquid ^He [Ref. 12 ]. Using this 
model, some preliminary calculations were done. Nevertheless, as it has been stated 
above, the results for polymeric elastomers deviate substantially. This can be expected as 
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the behavior of the viscoelastic polymeric material is far more complex than this simple 
single relaxation time model can fully describe. 




Figure 2. Frequency dependence of k' and k" 



3. Complex Modulus, Loss Tangent and Quality Factor 

In a dynamic measurement, the characteristic elastic modulus of a viscoelastic 
material can be represented as a complex quantity, G* , [Ref. 10 : p.32] 



G* = G'+jG" 
tan 5 = 

G' . 



(18) 

(19) 



The storage modulus, G , is defined as the stress in-phase with the strain in a 
sinusoidal shear deformation divided by the strain. It is a measure of the energy stored and 
recovered per cycle. 
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The loss modulus, G \ is defined as the stress 90° out-of-phase with the strain 
divided by the strain. It is a measure of the energy dissipated or lost as heat per cycle of 
sinusoidal deformation. 

The loss tangent ,tan 5, is thus a measure of the ratio of energy lost to energy 
stored in a cyclic deformation process. It is also equal to the reciprocal of quality factor, Q. 
The quality factor of a system can be defined by [Ref. 7 : p.l6] 

e=x 

( 20 ) 

where / is the resonant frequency and Af is the full bandwidth over which the average 
power has dropped to one-half its resonance value. In terms of relaxation time of the 
system, t, the quality factor can be expressed as [Ref. 7 : p.l6] 

Q = (21) 

or in a simple harmonic oscillator, as a function of mass, m, and mechanical resistance, R, 
as 

Q _ 2itfm 

R . (22) 

The quality factor can also be expressed as [Ref. 13] 

Q _ 

, (23) 

where is the energy stored in the resonance and is the average rate of energy 
dissipation. 

If u is the speed of the vibration, F is the amplitude of the constant applied force 
and m is the mass, at resonance, 

u = £ 

R, (24) 

and thus. 



11 



(25) 



R = 



_ F 



u , 



combining equations (22) and (25), we arrive at 



Q=Afu 

where A is a constant of proportionality 



(26) 



and hence, 



A _ liun 

■ ~r, 



fu 



4 . Dynamic Behavior of Viscoelastic Material 

From equation (17), the storage modulus, k ' , can be expressed as 



k = k\ ■¥ 






(27) 



(28) 



k2 + {(oRiY ( 29 ) 

As (D—¥0,k'—^kj, and as (D k' k] + k2- As stress is applied at low frequency, 

the delayed response of the material occurs within the period of the stress reversal. Near 
equilibrium is achieved at all time. However, as the frequency of the applied stress 
increases, a point wUl be reached when the response of the material cannot "keep up" with 
the stress and its deformation reduces. Consequently, as the frequency of applied stress 
increases, the storage modulus increases substantially. Moreover, energy loss, which 
can be measured by the phase angle between the stress and strain ( or the loss tangent ), 
behaves differently. At low frequencies, where the phase angles are small, the energy 
losses are small. When frequencies are high, the strain cannot response fast enough to the 
applied stress, so energy loss remain small too. Thus, the energy loss is highest between 
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the two frequency limits, when the phase angle and the strain amplinide assume relatively 
large values. When temperature increases, the mobility of molecules increases, and the 
strain can track the dynamic stress more closely. With lower temperatures, stiffening of 

the material occurs. Therefore, the same effect can be achieved by either of the following : 

1 . Increase of temperature or decrease of frequency. 

2 . Decrease of temperature or increase of frequency. 

This correspondence allows one to experimentally determine the behavior of 
materials at frequencies higher than typical apparatus may accommodate by decreasing the 
sample temperature. Likewise, an experiment can obtain very low frequency performance 
predictions by measuring the response of sample at elevated temperatures. 

Factors like molecular weight, composition of the elastomer, pressure, etc, 
have their effects on the dynamic properties of the viscoelastic materials but they are not 
included in this study. [Ref. 14] 

As the dynamic response of a viscoelastic material depends on the temperature 
and frequency [Ref. 10 : Chap. 11 ] , it appears to be most convenient if one of the two 
variables is being held constant while the other varies. Nevertheless, only a small range of 
viscoelastic behavior can be observed over the experimentally accessible frequency range, 
therefore a complete characterization of the material and its viscoelastic transition is only 
possible through measurements over a limited range of frequencies at various temperatures. 
This is the most popular approach to characterize viscoelastic materials in the field today, 
and it is based on the method of reduced variables developed by Williams, Landel and 
Ferry [Ref. 10 : pp. 294 - 320]. The treatment yields a master curve relating a chosen 

dynamic property such as storage modulus to a reduced frequency through a temperature- 
frequency shift factor, aj . Each material has its own expression for this temperature- 

frequency factor. The development of this temperature-frequency factor is discussed in 
depth in the original paper of Williams, Landel and Ferry [Ref. 15 ]. 
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In the resonant bar technique, it is required that measurement be made at discrete 
frequencies which correspond to the normal-modes of the bar. When the temperature of 
the bar varies, the resonant frequency changes also, and it is tracked by a phase-locked- 
loop (PLL) which will be discussed in Chapter in. The two variables of interest, i.e., 
temperature and frequency, are coupled and this presents a challenging problem in 
presenting and analyzing the experimental data. This problem is rather unique. In the more 
established technique using an accelerometer as a receiver and a shaker table as a driver 
[Ref. 16], the excitation frequency can be chosen independently of the temperature. 
Therefore, it can be seen that the resonant bar technique presents a simple solution to the 
problem of characterization of material under dynamic stress but with variables that are 
coupled together. In contrast, the direct measurement method using accelerometer is able to 
isolate the variables but it is mathematically more complex due to the large mass loading 
induced by the accelerometer and the attachment to the shaker table. To verify the 
apjplication of this technique in the investigation, an effort was made to select a material that 
has been investigated by authorities in this research area, in particular, by Capps of the 
Naval Research Laboratory - Underwater Sound Reference Detachment (NRL-USRD). 
The experimental data obtained in this investigation are converted to reduced frequency 
using the empirical formulae given in Capps' report [Ref. 16 ]. 
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III. EXPERIMENTAL APPARATUS 



A. INTRODUCTION 

The resonant bar theory is based on simple wave equations. Experimentally, it is also 
easy to set up, very repeatable due to the free-free boundary condition, and easy to operate. 
Moreover, with the electrodynamics transduction scheme, one can selectively excite 
torsional, flexural, or longitudinal modes using the same sample, transducer pair, and 
experimental set up and thus it is extremely versatile. There is no other known apparatus 
that can excite all three modes with the same set up, same transducers and with virtually no 

mass loading effects. Basically, the set up consists of the following components : 

1 . Test sample. 

2. Two transducers mounted on the sample. 

3 . Test jig for supporting the sample and providing the magnetic field for the 
transducer cods. 

4 . Electronic instmmentation and computer/controller. 

B. SAMPLE 

1 . Preparation of Materials 

With the exception of PMMA, the samples used in the study each consisted of 
two potting components, the resin (part A) and the hardener (part B). They are mixed 
together in proportions recommended by the manufacturers. The epoxies E-CAST® 
F82/215, PR1592 and EN9, were heated and degassed at 50'C before mixing. Epoxies 
PR1570 and EN5, were mixed and degassed at room temperature. The epoxies are heated 
and degassed in an oven equipped with a vacuum pump. The data sheets for the epoxies 
are included in Appendix B. 
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2. The Casting of Sample 

The cylindrical rod shaped samples are typically 30 cm in length and 1.27 cm 
(0.5") in diameter, so that a large Ud ratio is achieved. The apparatus can accommodate 
samples of dimensions with diameters of 1.2 ± 0.5 cm and lengths of 30 + 10 cm. The 
samples are made by pouring the prepared elastomers into a mold. The molding tubes for 
these samples are thick walled (1/8") Teflon® tubing with inner diameter of 0.5 ± 0.05". 
Teflon® is used as the mold material due to its release properties. It does not require any 
mold release agent in the preparation of a sample. 

3. Curing 

After the polymer is poured into the mold, the epoxy is cured according to 
manufacturer's recommendation. When the epoxy is properly cured, the sample is pushed 
out from the Teflon® tubing. A post cure (at 50‘C for 24 hours) is carried out so that the 
properties of the sample reach their final stable state. Figure 3 shows the sample cast from 
epoxy EN9 with attached transducer coils and thermistor. 

C. TRANSDUCER COILS 

The transducers that act as the driver and receiver of the resonant bar are identical in 
construction. Each consists of No. 32 gauge copper wire with a nominal diameter of 0.2 
mm which are approximately one meter in length, and wound into a 10-tum coil. Each coil 
weighs about 0.56 gram. TTiey are glued to the ends of the sample by a general purpose 
plastic cement ( GC Electronics 10-324 ). The ends of the coil are stripped of their 
insulating layer, and are then coated with solder to assure good conductivity. After each 
coil is prepared, a continuity check is done to ensure that it is capable of driving and 
receiving of signals. The coil typically has a resistance of about 1.1 12. If a large signal is 
desired, the turn ratio can be increased by a factor of N, thus increasing the the signal by 
at the expense of added mass of coil. If this added mass becomes significant, a simple 



16 




effective length correction factor can be included in the analysis of the result [Ref. 3]. The 
elastic properties of the plastic cement are not expected to affect the observed resonance and 
in-phase magnitude of the bar. 



FIgu re 3. Photograph of a prepared EN9 material sample. 



D. APPARATUS 

The apparatus used in this research, shown in Figure 4, was constructed "in-house" 
based on a commercial apparatus [Ref 4]. The apparatus is made up of an excitation 
system, an adjustment system and a suspension system. The electrodynamic excitation 
system has two pairs of permanent magnets. They provide the necessary magnetic field, 
so that a varying sinusoidal electrical signal can be converted into a mechanical driving 
force or torque and the reciprocal effect. The gap between the magnetic poles is 
approximately two centimeters and it is adjustable to accommodate samples of different 
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diameters. The magnetic field strength in the gap is approximately 2.4 ±0.1 KOe (0.24 
Tesla) over a temperature range of -17'C to 85°C. The temperature dependence of the 
magnets was determined to be -2.3 Oe/'C. One pair of the magnets is attached to a sliding 
platform to accommodate samples of different lengths. 

The suspension system for the material sample consists of a pair of "X" supports with 
rubber bands in a crisscross manner. The rod-shape sample is then placed on the two 
supports so that a free-free boundary condition can be realized. The "X" supports can be 
raised or lowered vertically and they can slide horizontally. This is useful as they can 
support the material sample at any desired points and heights. The smallest effect from the 
suspension supports is achieved when the rod is placed so that the vibrational nodes rest on 
the "X" supports. This is not necessary for high loss materials due to the fact that the 
suspension losses are relatively low as has been quantified in this research (see Chapter IV 
Section C). 

E. ENVIRONMENTAL CHAMBER 

An environmental chamber is used to control the temjjerature of the material sample 
under test. The model used is BHD-408 Bench-Top Temperature and Humidity Test 
Chamber manufactured by Associated Environmental System. It has a internal dimensions 
of 24" ( Height ) x 24" ( Width ) x 24" ( Depth ). The temperature range of the system is 
-17°C to 85*C and the temperature control stability is ±0.5°F (approx. ±0.3'C). It has a 
dual loop controller capable of independent operation, or it can be controlled by a host 
computer through a IEEE 488 Interface. For the controller, a computer like the HP9000 
series computer can be used. For our application, the local control mode is used. The 
detailed specification of the chamber is provided in Appendix C . 
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Figure 4. Resonance based moduli measurement apparatus. 



F. SELECTION AND EXCITATION OF A DESIRED MODE 

The material sample is positioned on the apparatus so that the ends of the sample and 
the transducer coils are in the vicinity of the maximum magnetic field strength during the 
dynamic moduli measurements. The mode of excitation can be selected by the orientation 
of the coils within the magnetic field. Figures 5, 6, and 7 shows the relative arrangement 
of the cod and magnetic field to produce the torsional, flexural and longitudinal vibration 
respectively. A detailed discussion can be found in Garrett's paper [Ref. 3]. 
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Figure 5, Torsional vibration excitation and detection scheme 

(a) . Arrangement for torsional vibration, the arrows indicate the 

direction of the electromagnetic forces on the coil for a given phase 
of current. 

(b) . The torsional frequency response of F-28 
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(b) 



Figure 6. Flexural vibration excitation and detection scheme 

(a) . Arrangement for flexural vibration, the arrows indicate the direction 

of the electromagnetic forces on the coils for a given phase of 
current. 

(b) . The flexural frequency response of F-28. 
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(b) 

Figure 7. Longitudinal vibration excitation and detection scheme 

(a) i^angement for longituduial vibration, the arrows indicate the 
direction of the electromagnetic forces on the coils for a given phase 
of current. 

(b) TTie longitudinal frequency response of F-28. 
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(,. ELECTRONIC INSTROMENTATIO" 



1 . Room Temperature Dynamic Moduli Determination 

The instrumentation required for determination of dynamic moduli at room 
temperature is shown in Figure 8. The HP3562A Dynamic Signal Analyzer is used as a 
signal source as well as a receiver, so that the frequency response of the material sample 
can be obtained. The excitation signal used is a swept-sine wave with a driving output 
from the amplifier (HP467A) of typically two to three volts peak-to-peak. Typical 
frequency response are shown in Figure 5(b), 6(b), and 7(b). 




Figure 8. Block diagram of the instmmentation used at room temperature. 

The HP3562A has a built-in curve fit function which is used in the study 
extensively. Using this pole-zero plot curve fit function, the quality factor, resonant 
frequency, and loss factor can be calculated. The theoretical basis for obtaining the 
resonance frequency and quality factor from the pole-zero plot is discussed in Brown, Tan 
and Garrett's paper in Appendix A. 
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2 . Temperature and Frequency Dependence of the Dynamic Moduli 



The instrumentation for this measurement is shown in Figure 9. The temperature 
of the material sample is varied by an environmental chamber. Using a lock-in amplifier 
and a voltage controlled oscillator (VCO) in a phase-locked-loop configuration, the 
resonant frequency of a selected mode of the material sample is tracked. When the 
temperature of the specimen changes, its resonance frequency changes due to the changes 
in the modulus. There is also a contribution to the change in frequency due to the change 
in the length of the sample via the coefficient of thermal expansion, but this effect is small 
by comparison and is neglected in this study. Using a HP9836C computer as a system- 
controller through the HPIB, readings of the temperatures, resonant frequencies of the 
sample and the in-phase component of the lock-in amplifier can be recorded. A listing of 
the control program is included in Appendix D. 

The sample rod is placed on the apparatus in the same manner as in the dynamic 
moduli determination at room temperature. The temperature of the rod is monitored by a 
thermistor (HP0837-0164) which is compatible to the HP3456A Digital Voltmeter. The 
voltmeter is set up to read the thermistor temperature in degrees Celsius ('C) with a 
resolution of 0.00 UC. 

The procedure of this measurement is to identify the modes of the bar at room 
temperature first, as discussed previously, and then select a given mode for automatic 
tracking with temperature. This is done by tuning the voltage controlled oscillator (VCO) 
manually to resonance with the integrator shorted and the error signal feedback path open 
so that there is no error signal presented to the voltage control (feedback) point. The phase 
shifter is then adjusted so that there is zero output from the quadrature signal channel of the 
lock-in amplifier. The integrator can then be op>ened and the control loop completed. If the 
resonance "mns away", the signal (either transmitter or receiver) needs to be inverted. This 
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Figure 9. Block diagram of the instrumentation used to track the dynamic moduli 

dependence on temperature and frequency. 
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is most easily accomplished by shifting the phase an additional ± 180* if a lock-in analyser 
is being used as the mixer/phase shifter/low pass filter. It may also be necessary to adjust 
the filter time constant and amplifier gain if the signal oscillates. Once the control loop is 
locked and stable, the temperature may varied. [Ref.3] 

After setting up the PLL, the environmental chamber is programmed to run a 
temperature profile as shown in Figure 10. The long time interval is used to achieve a 
quasi-equilibrium state, so the material sample is in equilibrium with the uniform 
temperature environment. 

The nature of the experimental set up calls for automation. Thus, in their 
material selection study [Ref. 6], Beaton has written a BASIC program using the HP9000 
series 200 computer Model 23 6C to acquire inputs of temperature and frequency data. The 
program is modified to include the input from the in-phase component of the lock-in 
amplifier and to extend the duration of the measurement from 40 minutes to 20 hours. 

To facilitate the graphical presentation of data, the data are converted from the 
HP data format to Macintosh Cricket Graph™ format. Using a software package called 
MacTerm for Macintosh II and software written by Professor Steven R. Baker of NFS 
(Appendix E) for HP9000 series computer, the data recorded in HP BOAT file are 
transported from the HP9000 computer via a RS 232C port to the Macintosh II computer. 
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Figure 10. Typical temperature profile for the dynamic mcxiuli measurement 
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IV. SYSTEM VERIFICATION 



A. INTRODUCTION 

To investigate the complex moduli of a material using the resonant bar technique, it is 
essential to quantify the losses due to the suspension system. As mentioned earlier, the 
computations of the moduli are based on a free-free boundary condition. This assumption 
is valid when the suspension has a negligible effect on the response of the rod. 

B. INITIAL INVESTIGATION 

An existing sample rod of E-CAST® F82/215 that was used in a previous experiment 
involving fiber-optics (see Aptpendix A) was initially chosen to investigate the losses due to 
the suspension system. The quality factor, Q, which is a measure of the energy stored to 
the energy lost per acoustic cycle (Chapter II, equation (23)) was measured, for various 
support positions (Figure 11) and two types of suspensions: rubber bands and fishing line. 
As the suspension positions one and three illustrated in Figure 1 1 correspond to the nodal 
position of mode two and three respectively, it is expected that the losses are lower than the 
non-nodal placement positions. 

The experimental set up is shown in Figure 8, Chapter IE. The HP3562A Dynamic 
Signal Analyzer is used to generate a swept sine signal which is used as the reference input 
(channel 1) for the HP3562A. This signal is amplified by HP467A power amplifier and 
then applied to the driver coU. The electrical signal is converted to a mechanical driving 
force through the magnetic field and the transducer coil electrodynamic interaction. The 
response of the material sample is then converted to an electrical signal through the 
reciprocal electrodynamic transduction mechanism. The received signal is amplified and 
filtered by an Ithaco 1201 pre-amplifier and input to the second channel of HP3562A. 
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Plotting the ratio of these two inputs as a function of frequency, the frequency response of 
the sample material can be obtained directly using the HP3562A. Typical frequency 
response curve is shown in Figure 5(b). Using the built-in curve fit function of the 
HP3562A, the quality factor and resonant frequency of each resonant mode can be 
determined from the pole-zero table generated by HP3562A (see Appendix A). 

From Table fra), it can be seen that quality factor of mode one and mode two is almost 
identical, but mode three and mode four differ from these two modes by approximately 5% 
and 10% respectively. Using the grand average as the reference, only mode four deviates 
by about 10%. It is noted that, suspension position three produces the highest quality 
factor in the third mode as expected, but suspension position one produces the lowest 
quality factor in the second mode instead of the highest which is not the expeaed result. 

From Table 1(b), response of mode three and four is equal within 5%. The quality 
factor of mode one is the lowest and for mode two it is the highest. Except for mode one, 
the other modes are less than 18% different from the grand average. Suspension position 
three has the highest quality factor within the third mode as expected and it is significantly 
higher than the average quality factor within the mode. Susp>ension position one does not 
produce the highest quality factor in second mode as it is expected, although it is the second 
highest among the quality factors obtained in the same mode. 

The experiment was repeated and the similar observations were made. Thus, within 
the experimental uncertainty, we may conclude that the suspension loss contribution is 
small and it will be substantial only if a very high quality factor material is under 
investigation. 
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Position 



1 



2 



3 



4 



5 



6 



7 



8 



9 




Figure 1 1 . Support positions used in suspension loss measurements 

and indicate the nodal positions of second and third mode respectively 
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TABLE 1. QUALITY FACTOR OF E-CAST® F82/215 OBTAINED FROM 
DIFFERENT SUSPENSION POSITIONS IN TORSIONAL 
VIBRATION 

(a) RUBBER BAND SUSPENSION 



Position 


1st Mode 


2nd Mode 


3rd Mode 


4th Mode 


1 


22.20 


21.79 


21.01 


19.83 


2 


22.83 


22.29 


21.48 


19.92 


3 


23.09 


22.51 


21.76 


19.71 


4 


22.12 


23.02 


21.41 


19.42 


5 


22.89 


22.47 


21.20 


20.37 


6 


22.01 


22.19 


21.24 


20.67 


7 


21.94 


22.32 


21.18 


20.45 


8 


22.49 


22.40 


21.17 


20.44 


9 


22.68 


22.28 


21.27 


20.48 


Average 


22.4 ± 0.4 


22.4 ± 0.3 


21.3 ±0.2 


20.1 ±0.4 



Grand average: 22 ± 1 
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TABLE 1. QUALITY FACTOR OF E-CAST® F82/215 OBTAINED FROM 
DIFFERENT SUSPENSION POSITIONS IN TORSIONAL 
VIBRATION 

(b) FISHING LINE SUSPENSION 



Position 


1st Mode 


2nd Mode 


3rd Mode 


4th Mode 


1 


15.47 


22.97 


20.89 


19.20 


2 


12.57 


22.84 


20.77 


19.21 


3 


10.55 


22.73 


21.59 


19.45 


4 


10.31 


23.12 


21.46 


18.88 


5 


14.71 


22.38 


20.43 


19.31 


6 


12.95 


22.30 


20.12 


20.10 


7 


17.78 


21.51 


19.77 


20.10 


8 


14.72 


21.92 


19.09 


19.79 


9 


12.05 


22.04 


20.60 


19.61 


Average 


13±2 


22.4 ± 0.5 


20.5 ± 0.8 


19.5 ± 0.4 



Grand average: 19 ±4 
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C. SUSPENSION LOSS DETERMINATION USING AN ALUMINIUM 

SAMPLE 

The total losses of a forced vibration system in the resonant bar technique is the sum of 
the intrinsic losses due to the sample's material properties, suspension losses due to the 
suspension system of the apparatus and radiation losses. Radiation loss is considered to be 
negligible and therefore, to be able to measure the intrinsic loss of a sample material, an 
upper bound must be placed on the suspension losses. 

For a high Q material, the resonant peak is very sharp and its -3 dB bandwidth is very 
small. If the Q of the material is of the order of 1000, then the required sweep time, 
oscillator stability, and sample temperature stability become a problem both with 
instrumentation and the patience of the experimentalist. Therefore, we can not simply 
replace the E-CAST® F82/215 by the 6061T6 aluminium rod and repeat the swept sine 
frequency response measurements using the HP3562A and expect to obtain accurate curve 
fit parameters (poles and zeroes) for the quality factor. 

The half-power point measurement of Q also depends on the frequency resolution and 
is not suitable due to the temperature stability of the resonance frequency. Another 
approach that "exploits" the high Q characteristics of the system was used. As the 
reciprocal of frequency is time, one can use a time measurement technique that has a long 
time record, and thus an equivalent high frequency resolution to measure the Qofa low 
loss material. Free decay measurements is one of such techniques that involves the 
measurement of successive attenuated amplitudes, from which the damping and hence the 
quality factor can be obtained. Since a long time record means large amount of data can be 
generated, regressional analysis can be performed to minimize any random errors due to 
decaying amplitude variation in the determination of the exponential decay time, T. 
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A block diagram of instrumentation for the free decay measurement is shown in Figure 
12. The sample is driven at one of the resonant frequencies and the tone burst excitation is 
abruptly removed. This is achieved by means of the trigger pulse switch which 
disconnects the driving signal to the transducer coil and at the same time provides a trigger 
pulse to the external trigger input of the Nicolet digital storage scope which records the 
response of the material sample. Using a digital storage scope and linear regression 
analysis, the quality factor of the sample can be determined. A typical graph of the 
measurement is shown in Figure 13. The results are summarized in Table 2 where the 
uncertainties in Q result from the uncertainty of the slope as determined by the least square 
method. Based on the average of multiple measurements, the systematic error is larger than 
the random error when the quality factors are this large and the losses are correspondingly 
small. From these results, it can be seen that the difference of the quality factors between 
the "best" and "worst" support positions varies from 24% to 43% for the first three modes. 

It can therefore be concluded that for a low loss ( high Q ) system, the suspension loss 
makes a significant contribution to the total losses of the suspension system/sample 
combination as currently configured. Using the worst case analysis, the suspension loss 
is about 43% for a system with aQ of about 1300. For a high loss (low Q ) system, such 
as the F82/215 sample, that has a Q of about 20 (and most viscoelastic materials with Q's 
typically less than ten, see Chapter V, Table 5), the suspension loss contributes less than 
1% to the total losses of the entire system. 
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Figure 12. Block diagram of the free decay measurement of aluminium rod 




TIME [sec] 



Figure 13. Free decay measurement of Q 
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TABLE 2. QUALITY FACTOR OF 6061T6 ALUMINIUM ROD OBTAINED 

FROM TWO DIFFERENT SUSPENSION POSITIONS USING 
FREE DECAY METHOD . 





1st Mode 


2nd Mode 


3rd Mode 


Best 

support 


1306 + 0.8% 
1338 ±0.4% 


868.7 ± 0.6% 
830.5 ± 0.9% 


1222 ± 1.4% 
1437 ± 0.8% 
1792 ± 1.3% 


Average ± o 


1320 ± 16 


850 ± 20 


1500 ±300 




Worst 

Support 


751 ±0.4% 


644.9 ± 0.4% 
653.5 ± 0.5% 


968 ± 5% 
940 ± 3% 
934 ± 4% 


Average ± o 


751 


650 ±5 


950 ± 18 


Relative difference 


43% 


24% 


36% 
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V. DATA 



A. INTRODUCTION 

This study consists of three main areas. The first is the verification of the resonant bar 
technique to determine the elastic properties of a sample material called E-CAST® F82/215. 
The second portion involved quantifying the suspension losses using E-CAST® F82/215 
and 606 1T6 aluminium bar sample. The results of the suspension losses were reported in 
Chapter IV. The third section is the investigation into the temperature and frequency 
dependence of the elastomeric materials: Uralite 3130, PR1592 and plexiglass or 
polymethyl methacrylate (PMMA). Uralite was chosen because of its use in fiber-optic 
hydrophones [Ref.l7], PR1592 because of an accessible NRL-USRD report [Ref.l6], and 
PMMA because of its wide use as a "standard" material in viscoelastic investigation. 

B. PRELIMINARY INVESTIGATIONS 

To gain experience with the resonant bar technique, some preliminary measurements 
were made on E-CAST® F82/215 to determine its dynamic moduli at room temperature. 
The result, which is summarized in Table 3, where p is the mass density, G is the shear 
modulus, and E is the Young's modulus, agreed with those obtained in Reference 6 which 
used a similarly constnicted resonant bar apparatus. 

TABLE 3 . DYNAMIC MODULI OF E-CAST® F82/215 





P 


G 


E 


a 




[kg/m^] 


[GPa] 


[GPa] 




Ref 5 


1120 


0.91±0.06 


2.5±0.2 


0.37 


Present study 


1120+20 


0.96±0.06 


2.70±0.05 


0.40 



37 



A fiber optic interferometer was constructed by embeding a fiber which comprised of 
one of the legs of a Michelson interferometer, into a E-CAST® F82/215 rod. In this study, 
the interferometric output is used as second means to determine the resonant frequency of 
the rod. It is being investigated for smart skin applications to monitor the state of strain of 
struts. The results were reported in the proceedings of the SPIE Conference "Fiber Optic 
Smart Stractures and Skins HI" in San Jose, September 1990. That paper is included as 
Appendix A. 

C. TEMPERATURE AND FREQUENCY DEPENDENCE OF 

VISCOELASTIC MATERIALS 

1 . Test Matrix 

In this study, Uralite 3130, PR 1592 rods cast in Teflon® tubing and PMMA 
rods were used. The torsional mode is the easiest mode to excite and detect, and thus most 
of the data collected were on the shear modulus. For the PR1592, data was also collected 
for the flexural resonant mode which depends on the Young's modulus. The test matrix is 
as follows 

TABLE 4. TEST MATRIX FOR VISCOELASTICITY MEASUREMENTS 





PR 1592 


PMMA 


U3130 


Torsional 


x/ 


x/ 




Flexural 




— 


— 


Quality Factor 






— 



2 . Data Presentation for PR1592, PMMA and U3130 

The resonant bar technique produced data with the two variables, frequency and 
temperature, coupled together. That is, neither was held constant through a measurement 
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run. This was dis assed in Chapter II. The temperature-frequency dependence 
measurement of the first four torsional modes of PR 1592 have been investigated in this 
study. Graphs produced from the data are shown in Figure 14 and Figure 15. Note that 
the points were taken at arbitrary time intervals, in this case, every five minutes. The 
resonance is tracked continuously, thus the discrete data points which are presented could 
have been obtained with different data density. It can be seen from Figure 14 that as 
frequency increases, the shear modulus increases in a quadratic manner. As the higher 
order mode is examined, the slope of the response curve reduces at each frequency. In 
Figure 15, the higher modes have a higher value of shear modulus at a constant 
temperature. Further, it can be observed that at the highest temperature, the shear modulus 
of aU the four modes approaches a constant value asymptotically as expected in a single 
relaxation time constant system discussed in Chapter n. The temperatures were not large 
enough to obtain the high temperature asymptotic behavior. 

The experimental data was then plotted in three dimensional space to visualize the 
relationships between shear modulus, temperature, and frequency. To compute the loss 
tangent, the quality factor, Q , of specimen was measured as a function of temperature. In 
order to present the data in the conventional format, Capps' [Ref. 16 : p. 210] and Ferry's 
[Ref. 10 : p.316] results were used to determine the temperature-frequency shift factor 
necessary to obtain the master curve. 

a. Three Dimensional Plots 

The effects of temperature and frequency on the behavior of a viscoelastic 
material can be most easily visualized in a three dimensional plot. The following sections 
wUl discuss the various data obtained to construct the three dimensional plots. 



(1) Shear Modulus. The three dimensional plots are shown in Figure 16, 
Figure 17 and Figure 18 for PR1592, PMMA and Uralite 3130 respectively. In Figure 16, 
all the data from the lowest four torsional modes of PR1592 are presented. 




Figure 14. Shear modulus of PR1592 plotted against resonance frequency (temperature 
is not constant). The numeric indicates the mode number. 
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Figure 15. Shear modulus of PR1592 plotted against temperature. 

It can be seen that all these are smooth, continuous curves. As 
temperature increases, the resonant frequency decreases and the shear modulus decreases 
correspondingly. Both Uralite 3130 and PR 1592 are "rubbery" materials at room 
temperature. The curvature of the projection of their three dimensional modulus graphs in 
the frequency-temperature plane is opposite to that of the PMMA which is in the "glassy" 
state at the temperature of this study. When the higher resonant mode is used, the resonant 
frequency increases and the span of frequency also increases, however, the amplitude of 
the in-phase output of the lock-in amplifier decreases. Lastly, it can be seen that, the shear 
modulus is not a simple function of temperature and frequency. 
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Figure 16. Shear modulus of PR 1592 as a function of frequency and temperature. 



42 




Figure 17. Shear modulus of PMMA as a function of frequency and temperature. 
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Figure 18. Shear modulus of U3 130 as a function of frequency and temperature. 
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SHEAR MODULUS [Pa] 




(2) Loss Tangent. The three dimensional plots of the loss tangent of 
PR1592 and PMMA are shown in Figure 19 and Figure 20 respectively. For PR1592, the 
presence of glass-rubber transition is evident. The graph of loss tangent has a maximum. 
There is no local maximum in the loss tangent plot of PMMA for the range of frequencies 
and temperatures tested. 

From Chapter II equation (28), it is shown that the ratio of quality 
factor, Q , to product of vibrational speed of the end of the bar and resonance frequency, 
/«, is a constant. Since the in-phase output, , of the lock-in amplifier is proportional to 
this speed, equation (28) can be rewritten as 

Q_ ^ . 

Vinf , (30) 

where is a constant of proportionality. Noting that. 



tan 5 = 




and combining equations (30) and (31), the loss tangent can be expressed as 



(31) 



tan 5 = 



1 

bVinf. 



(32) 



By measuring the quality factor, Q, emd in-phase voltage at a set of 
temperatures and frequencies, the loss tangent empirical constant b can be obtained at the 
start of a test and then the loss tangent of the sample material can be determined 
continuously. One set of results used to obtain a value for the constant b is shown in Table 
5. For this example, the constant has a value of 64 milliseconds per volt (msec/volt). The 
standard deviation of the 13 measurements is 8 msec/volt so the uncertainty is 2.3 
msec/volt or about ±4%. The reader should be cautioned that this constant is a function of 
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the position of the coils and the gain of the pre-amplifier and lock-in amplifier. Thus, it 
must be determined for each sample and instrtiment setting used. 

(3) Young's Modulus. The plot of Young’s modulus of PR1592 is 
shown in Figure 21 . It was obtained using the flexural mode of the sample material. 
b . Master Curves 

The conventional way to present viscoelastic behavior is to use the master 
curve. The master curves are produced from reduced frequency for the shear modulus or 
loss tangent. The reduced frequency is calculated from data of Capps [Ref. 16 : p. 210] 
and Ferry [Ref. 10 : p.316]. 

The general form of the temperature-frequency shift factor, aj, based on the 
WLF equation [Ref. 10 : chap. 1 1], for a great variety of polymers can be expressed as 



log aj = 



-8.86 (T-T,) 

101.6 +r-r. 



(33) 



where T is the absolute temperature in degrees Kelvin K, and , is the reference 
temperature in K. The equation is applicable at temperature about 50* C above the glass 
transition temperature. However, as discussed in Ferry [Ref. 10 : chap. 1 1] , it is better to 
determine specific values of the two constants in this equation for each sample. Therefore, 
to explore this particular area, constants from the measurements of Capps for PR 1592 and 
Ferry for PMMA were used. Their expressions are as follows 



, -12.9 (T- 283.15 ) 

107 +r- 283.15 , 

for PR1592 and 



(34) 
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Figu re 1 9. Loss tangent of PR1592 as a function of frequency and temperature . 
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LOSS TANGENT 




Figure 20. Loss tangent of PMMA as a function of frequency and temperature. 
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LOSS TANGENT 



TABLE 5. 



DETERMINATION OF THE LOSS TANGENT 
EMPIRICAL CONSTANT " b" , 



T 

Temperature 

[C] 


f 

Frequency 

[Hz] 


in 

In-phase voltage 
[DCV] 


Q 

Quality factor 


Q/fV . 
in 

[msecAH 


-17.138 


801.9 


-0.1125 


5.523 


-61.22 


-10.124 


678.4 


-0.1048 


3.825 


-53.80 


-0.224 


515.2 


-0.1140 


4.378 


-74.53 


0.085 


511.5 


-0.1147 


2.913 


-49.65 


9.321 


398.7 


-0.1372 


2.927 


-53.51 


19.430 


306.5 


-0.1842 


3.415 


-60.49 


29.480 


259.4 


-0.2452 


3.847 


-60.48 


39.342 


231.6 


-0.3142 


5.082 


-69.84 


49.096 


217.1 


-0.4092 


6.304 


-70.97 


58.929 


201.7 


-0.5401 


7.805 


-71.66 


69.546 


195.0 


-0.7048 


9.701 


-70.58 


79.561 


190.0 


-0.9350 


12.41 


-69.86 


84.277 


188.1 


-1.0700 


13.60 


-67.58 



Average: -64 ± 8 msecA^ 




Figure 21. Young’s modulus of PR 1592 as a function of frequency and temperature. 
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VOCnVG'S AfO£>ULC/S fPaJ 



( 35 ) 



log aj = 



-21.5(7-211 ) 
43.1 +7-211 , 



for PMMA. The parameter, aj, was not independently determined in our study since we 
did not have sufficient time to use samples of different length. Samples with different 
lengths would allow moduli at the same temperature but with different frequencies, which 
were sufficiently well separated to determine the temperature-frequency shift factor with 
adequate precision. 

(1) Shear Modulus and Loss Tangent. The plots of the master curve for 
the shear modulus of PR 1592 and PMMA are shown in Figure 22 and Figure 23 
respectively. It should be noted that, for a one order-of-magnitude change of frequency 
and a temperature range of 100' C (approximately 30% change in absolute temperature), 
there is a ten-order-of magnitude change of the corresponding reduced frequency. This is 
rather extraordinary, and it demonstrates the strong dependence of the dynamic modulus of 
viscoelastic materials on temperature. 

From Figure 22, it can be clearly seen that there is a transition from 
mbbery plateau to glassy zone for PR 1592. This is evident from the maximum in the loss 
tangent curve which is a characteristic of a glass-rubber transition region. By inspection, 
for PR 1592, the transition has occured at a reduced frequency of about 100 Hz, assuming a 
reference temperature of lO'C. 

However, looking at the data obtained from the PMMA (plexiglass) 
sample, only the glassy zone and part of the glass-rubber transition region are observed. 
This is also evident in the loss tangent curve in which the loss tangent increases 
monotonicaUy as the reduced frequency decreases. 
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SHEAR MODULUS [Pa] 



0.3 




Figure 22, Master curve of PR 1592 based on the measurement of the first torsional 
mode. 
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LOSS TANGENT 



SHEAR MODULUS [Pa] 




REDUCED FREQUENCY [Hz] 



Figure 23. Master curve of PMMA based on the measurement of the first torsional 
mode. 
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LOSS TANGENT 



TABLE 6. COEFFICIENTS OF CURVE FIT FOR PR1592 TORSIONAL MODE 



Coefficient 


1st mode 


2nd mode 


3rd mode 


4th mode 


aO 


7.5324 


7.5210 


7.5607 


7.5529 


al 


0.13067 


0.13592 


0.13998 


0.13972 


a2 


1.0742e-2 


1.1154e-2 


1.1289e-2 


1.1633e-2 


a3 


-1.162e-3 


-1.159e-3 


-1.949e-3 


-1.757e-3 


a4 


-1.679e-5 


2.309e-5 


5.790e-5 


1.876e-5 


Correlation 

Coefficient 


1.000 


1.000 


1.000 


1.000 



TABLE 7. COEFRCIENTS OF CURVE FIT FOR LOSS TANGENT (TORSIONAL 
MODE) 





PR 1592 


PMMA 


Coefficient 


1st Mode 


2nd Mode 


1st Mode 


aO 


-0.63667 


-0.60388 


9.0104 


al 


7.2174e-2 


7.8122e-2 


4.2034 


a2 


-1.8544e-2 


-2.1478e-2 


0.61220 


a3 


1.6625e-3 


1.3403e-3 


3.9509C-2 


a4 


-1.1907e-4 


-4.5023e-5 


9.8409e-4 


Correlation 


0.998 


0.998 


0.999 


Coefficient 
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The master curves for PR 1592 are fitted to a fourth order polynomial 
function given in equations (36) and (37) and the coefficients for the shear modulus is 
tabulated in Table 6. The coefficients of the loss tangent of PR1592 and PMMA are given 
in Table 7. 



LogG = Logfr , n = 1 , 2 , 3, 4. 



(36) 



Log (tan 5)= ^ a„( LogfT , n = 1,2,3 ,4. 

n 



(37) 



By inspection of Figure 22, we see that the location of the inflection 
point in shear modulus occurs at approximately the maximum of the loss tangent as would 
be expected for a simple relaxation model discussed in Chapter n. Differentiation of the 
fourth order polynomial fits to the shear modulus and loss tangent yields an inflection point 
at log ajf= 2.85 for shear modulus and a maximum at log a-jf= 2.65 for loss tangent. 
This is considered to be in excellent agreement for the span of log a which is from -3 to 
7, over ten order-of-magnitude. 

(2) Young's Modulus. As for the shear modulus, the master curve for 
Young's modulus uses the same fourth order polynomial curve fit. The result is shown in 
Table 8. 
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TABLE 8. 



COEFFICIENTS OF CURVE FIT FOR PR1592 FLEXURAL 
MODE - YOUNG'S MODULUS 



Coefficient 


1st mode 


Capps' result 


aO 


8.0219 


7.71248 


al 


0.12331 


0.12531 


a2 


9.985e-3 


2.1141e-2 


a3 


-8.876e-4 


-1.8216e-3 


a4 


-4.27 le-5 


-2.1893e-5 


Correlation 

Coefficient 


1.000 


- 



These results are plotted in Figure 24 in which Capps' curve was 
overlaid on the regression curve obtained from the flexural experiment. We may conclude 
that the deviation between the two curves is greater than the experimental uncertainty of the 
data. Therefore, no definite conclusion can be made from this comparison as there are too 
many possible reasons for the deviation. The samples are not be identical in every aspect, 
namely, in preparation and curing and hence the assumption that aj is the same in our 
sample may not be true. It might be advised to obtain a sample previously tested by a more 
conventional method such as Capps' [Ref. 16] and collect data using the technique outlined 
in this study. 
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Log (Reduced Fre<i) 



Figure 24. Comparison of experimental data with data from Capps. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 



The purpose of the research reported in this thesis was the extension of the free-free 
resonant bar technique described by Garrett [Ref. 3], which measured the shear and 
Young's moduli of elastic materials, to include the measurement of complex moduli of 
viscoelastic materials, particularly those materials which are important in the fabrication of 
hydrophones. 

A. LOSS MODULUS MEASUREMENT TECHNIQUE 

Since much of the utility of the free-free resonant method is due to its simplicity, both 
in the acquisition and analysis of the data based on measurement of the frequencies of the 
resonance modes of the sample, a simple extension of the method was devised for 
determining the loss tangent. This was accomplished by using the in-phase output of the 
phase sensitive detector to determine loss while the quadrature (out-of-phase) component 
was integrated and feed back to a voltage controlled oscillator to maintain the system at 
resonance. The validity of the relation used to convert in-phase output voltage to loss 
tangent (or its inverse, the resonance quality factor, Q ) was verified by obtaining pole-zero 
curve fits to the entire resonance response curves as a means of independently deterrnining 
the loss tangent. The experiment dependent (but temperature independent) conversion 
factor was found to be constant to within a few percent for the samples which were 
investigated. 

Ehie to the interest in the measurement of loss modulus of polymeric samples, the 
contribution to the loss due to the sample suspension system, e.g., rubber bands, were 
quantified by measuring the change in quality factor as a function of position of the sample 
suspension points. This was done with a low quality factor (=20) elastomeric sample 
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using a pole-zero curve fit and with a high quality factor (=1000) metallic sample using a 
free oscillation decay time measurement. Both techniques unambiguously established that 
the suspension losses were minimized when the sample was suspended at vibration nodes 
and that the maximum effect on the measured loss modulus due to the suspension system 
losses for elastomeric samples with quality factors less than 20 was less than 1%. For 
samples with intrinsic quality factors of order 1000, the suspension losses could in the 
worst case be as large as 40% as the support system is currently configured. 

B . MODULI OF VISCOELASTIC SAMPLES 

The technique was used successfully to measure the elastic moduli (both shear and 
Young's) and loss tangent for three elastomeric samples (Uralite 3130, PR1592, and 
PMMA) over a temperature range of -17*C to +85”C. The measured data for PR1592 and 
PMMA were converted to master curves using the frequency-temperature shift factors for 
those materials from the published literature. Over that temperature range, the PR 1592 
sample exhibited a complete viscoelastic transition with a peak in the loss tangent and an 
inflection point in the shear modulus which were observed at the same value of reduced 
frequency to within experimental error as would be expected from a simple viscoelastic 
model. 

The absolute value of the Young’s modulus for PR 1592 differed from the master 
curve for that substance measured by NRL-USRD by an amount which exceed 
experimental error. The cause of this discrepancy has not yet been determined. One 
possible reason is that the samples were not of identical composition (molecular weight, 
cure time, etc.). We hope that this discrepancy can be resolved in the future by an 
exchange of samples between NPS and NRL-USRD. 
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C. ADVANTAGES OF THE FREE-FREE RESONANT TECHNIQUE 

One conventional technique for making the measurement of storage and loss moduli 
involves the longitudinal excitation of a rod-shaped sample at one end by an electrodynamic 
shaker table and the detection of the rod response at the other end by a piezoelectric 
accelerometer. The complex transfer function thus obtained is used to derive the Young's 
modulus of the sample. This is accomplished by solving a complex transcendental 
equation which results from the fact that the mass of the accelerometer is a non-negligible 
fraction of the sample mass and the fact that the data obtained is, in general, at frequencies 
below resonance. The conventional technique also only obtains the Young's modulus. A 
different set-up, sample, and transducers are required to measure the shear modulus. 

The free-free resonant bar technique is far simpler and all of the advantages discussed 
by Garrett [Ref.3] for this technique when used in the measurement of the shear modulus 
also apply to the measurement of the loss modulus. These include (1) the simplicity of the 
conversion of resonance frequency to modulus since the mass loading due to the 
transducers is small; (2) the fact that the transducers, consisting only of magnet wire, are so 
inexpensive that they can be left on the sample; and most importantly, (3) the ability to 
selectively and strongly excite all three modes of the bar (longitudinal, torsional, and 
flexural) and detect them with a high value of signal-to-noise using the same transducer. 
This ability to measure all three modes means that (4) both the Young's and shear moduli 
are available from the same apparatus and (5) there are two modes which yield values for 
Young's modulus so there is an intrinsic self-consistency check. The use of a phase 
sensitive detection scheme allows (6) the resonance to be tracked automatically and 
continuously by using the quadrature signal in a feed back circuit to maintain the bar at 
resonance as the temperature of the sample is varied. This thesis established that the in- 
phase signal can also be used to (7) automatically and continuously track the loss tangent as 
a function of temperature. 
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D. DISADVANTAGE OF THE FREE-FREE RESONANT TECHNIQUE 

At the present time there seems to be only one significant disadvantage to this modulus 
measurement method in comparison with the conventional transfer function method 
described by Capps [Ref. 16]. This problem arises when one attempts to determine the 
temperature-frequency shift factor that is required to present the modulus data in the master 
curve format. In its present implementation, the free-free resonant technique described in 
this thesis provides modulus data at resonance frequencies which necessarily change with 
changing temperature. The determination of the shift factor is simplified if the modulus is 
obtained at fixed frequencies and a variety of temperatures. 

In principle, the free-free resonant technique can provide the modulus at several 
frequencies for each temperature since more than one resonant mode can be observed (four 
modes were tracked for torsional vibrations of PR 1592) but at this time we have not been 
successful in developing an algorithm which can extract the shift parameter from these 
multiple modes. One could, of course, obtain a variety of frequencies at all temperatures 
by using samples of different length but this would increase the measurement complexity 
by requiring more than one sample. 

E, RECOMMENDATION FOR FURTHER STUDY 

Based on the conclusions reported in this chapter, the two most interesting and 
important follow-on research problems would have to be (1) the derivation of an algorithm 
for the extraction of the frequency-temperature shift factor for the modulus measured by the 
temperature dependance of the fundamental resonance and its next few overtones and (2) a 
"round robin" calibration of "standard" samples using both the free-free resonant technique 
and the conventional transfer function technique. 
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APPENDIX A. CONFERENCE PAPER PRESENTED IN SPIE 
CONFERENCE "FIBER OPTIC SMART STRUCTURES AND 
SKIN III" IN SAN JOSE, SEPTEMBER 1990. 
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Nondestructive Dynamic Complex Moduli Measurements Using a 
Michelson Fiber Interferometer and a Resonant Bar Technique 

David A. Brown, Beng-Hock Tan, and Steven L. Garrett 
Naval Postgraduate School, Physics Department, Code PH/Gx 
Monterey, CA 93943 

ABSTRACT 

An optical fiber, used as a leg of a Michelson interferometer, is cast into a long cylindrical bar of E- 
CAST F-82 epoxy^. The bar can be selectively excited in any of its lowest flexural, torsional, and 
longitudinal modes. The interferometer is used to detect the resonant modes of the "free-free" bar and 
from these modes both tne Young's and shear elastic moduli are determined. The complex modulus is 
determined by measuring the quality factor (Q) for each resonant mode. The measurement technique is 
entirely nondestructive and yields results of the two independent moduli with the same transducers. 

1. INTRODUCTION 

4 

This article involves the use of fiber-optic sensors in the non destructive resonant determination of 
elastic moduli of epoxies. A "free-free" bar is selectively excited in its flexural, torsional, and 
longitudinal vibrational modes using a transducer consisting of .coils of magnet wire placed in the 
magnetic field created by a pair of permanent magnets, the details of which are covered in an article 
entitled "Resonant Acoustic Determination of Elastic Moduli"^. A commercially available apparatus^ 
was used to excite the modes of a bar which had embedded in it an optical fiber acting as the sensing leg 
of an interferometer. The bar is placed on a pair of soft rubber bands so that the ends are free to move. 
The phase modulations induced in the optical fiber are converted interferometrically into intensity 
modulations at the coupler location. The output of the interferometer at the photodetector is monitored to 
accurately determine the resonant frequency (and its overtones).of the particular mode of excitation of the 
bar. The square of the frequencies of the flexural and longitudinal resonant modes are proportional to 
the Young's modulus and the square of the frequency of the torsional modes is proportional to the shear 
modulus. The quality factor, or Q, of the resonant modes is equal to the ratio of the real to imaginary 
parts of the complex moduli and the inverse of the characteristic loss tangent, tan 5. Since the 
measurement technique is resonant, the signal to noise ratios are typically very high and the modulus that 
is obtained is a dynamic complex modulus at the frequencies corresponding to the fundamental bar 
resonance and its overtones. The modes are also determined by a second coil electrodynamically as 
described in reference 2. 
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2. EXCITATION AND DETECTION 



The differential Lorentz, dF, force produced on a segment of wire, dl, carrying a current, I, in a static 
magnetic field, B, is given by 

dF=IdrxB. (1) 

Longitudinal, torsional, or flexural forces can be generated in order to selectively excite each of the three 
vibrational modes with the apparatus pictured in Figure 1. The particular mode excited depends on the 
relative positioning of the wire coils carrying the current, I, and the direction of the magnetic field. 
Typically the magnetic field direction and strength created by the pair of permanent magnets as well as 
the current driven through the coil of wire are constant and independent of firequency. When the 
frequency of the oscillator driving the wire coil is varied the bar is excited in its characteristic resonant 
modes of vibration. The detection of these modes is typically an easy task and can be accomplished by 
placing a second wire coil at the opposite end of the bar within a magnetic field created by a second pair 
of permanent magnets. The output of this coil of wire transducer is an EMF (voltage) which is 
proportional to the change in magnetic flux linking the coU and is given by 

V = -^| B • n dA (2) 

For a small segment of wire moving with velocity, u in a magnetic field B, the induced EMF is given by 




Figure 1. The MMD measurement apparatus^ used to electrodynamically excite and detect the flexural, 
longitudinal, and torsional modes of a bar in order to determine the shear and Young's elastic moduli. 
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The resonant vibrational modes can also be detected interferometrically due to the sinusoidally excited 
strain induced phase shifts in the leg of the interferometer. The output of a pigtailed single mode laser 
(Sharp LTOlO, 818 nm) is delivered to a 2 x 2 fiber optic spUtter/coupler where it is divided into two 
separate legs of a Michelson interferometer. One leg is cast axially a constant radial distance from the 
center of the rod and the other acting as a reference placed on the lab bench. The interferometer is 
cleaved at the end of the rod so that the coherent light reflects and combines interferometrically at the 
coupler which is characteristic of a Michelson interferometer. The resulting characteristic fringe pattern 
is observed with an oscilloscope from the output of a photodetector as illustrated in the measurement 
setup in Figure 2. 




Figure 2. Illustration of the measurement setup for the electrodynamic excitation and both electrodynamic 
and interferometric detection of the modes of a bar. 

3. THEORETICAL RESONANCE FREQUENCY 

Once the resonances have been determined the appropriate moduli, Young's, E, or shear, G, can be 
determined from the equations in this section provided that the dimensions tmd density of the bar can be 
measured. The resonance frequencies will also depend on the boundary conditions which for this 
apparatus is "free-free" corresponding to zero stress and zero moment at both ends. 

A uniform, cylindrical rod-shaped sample of a homogeneous, isotropic solid having circular cross- 
sectional diameter, d, and length, L, which is significantly greater than its diameter, will propagate three 
independent waves if their wavelengths, X, are much greater than d. The displacements associated with 
the longitudinal and torsional modes satisfy a partial second-order wave equation, and for a free-free 
boundary condition, the resonances are harmonically related. The Young’s modulus can be expressed 
in terms of these parameters as 

^ = (4) 

where n is a positive integer. Similarly, the free-free boundary condition leads to a series of harmonic 
torsional modes whose frequencies, are related to the shear modulus by 
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V n / . 



(5) 



The measurement of the flexural mode is not necessary but it does provide a second estimate for the 
Young's modulus and its fundamental frequency is typically an order-of-magnitude lower than the 
longitudinal modes. The flexural mode is also strongly excited by this electrodynamic transduction 
scheme so it is easy to observe. 

Unlike the torsional and longitudinal modes, the flexural waves of the bar obey a fourth-order 
differential equation and the flexural wave phase speed is dispersive. The application of "free-free" 
boundary conditions in this case leads to a series of modes which are overtones but not harmonics. The 
Young’s modulus of the rod can be expressed in terms of the flexural modes as 



£_ 1024 

7^2 ^2 _ 



( 6 ) 



where n takes on the values 3.0112, 4.9994, 7, 9, 11... 



4. CONCEPT OF A COMPLEX MODULUS 



4.1 The Loss Tangent and Quality Factor 



Many solids and most metals behave as nearly perfect elastic media and obey Hooke's Law. For the 
application of a dynamic force to these materials, there is a corresponding deformation that, in the linear 
limit, is proportional to the magnitude of, and in phase with, the applied stress. The ratio of the applied 
stress to the resulting normalized deformation or strain is therefore a characteristic property of the 
material which is termed the elastic modulus. The modulus of these metals may be independent of 
frequency over a large frequency range. If a material possesses sufficient damping characteristics, (or 
equivalently low Q’s) energy can be dissipated as heat in each cyclic deformation resulting in a strain that 
is out of phase an angle 5 with the applied stress. The characteristic elastic modulus can thus be 
represented as a complex quantity in phasor notation with magnitude G and phasor angle 5 ,or as an in 
phase, G’, and quadrature term, G" as"**^ 



G* = GeJ^°=G’ -t-jG" (7) 

where G"/G' = tan 5 g and j = V-1. The loss tangent is related to the mechanical Q by the relation: 



r;= U /U = ^ 

° Q (8) 

To illustrate the concept of the complex modulus further, let us st^ with the familiar equation of 



tan5Q=G"/G’ = 

le concept of the c 
motion for a damped driven oscillator 



M 




+ R 



^+K|(0 = F(t) 



( 9 ) 



where ^ is the displacement from equilibrium, M is the mass of the oscillator, R is the resistive 
mechanical damping coefficient, K is the spring constant or characteristic stiffness, and F is the applied 
force. If the applied force is sinusoidal of frequency O), the above equation is expressed as complex 
variables in the form 



-coM^ +K 



1+i 



coR 



K 



^ =f’ 



( 10 ) 



where the coefficient of ^ (the effective stiffness term) can also be expressed as a complex variable^ 
K*= Kjl+j^^J= K{l+jtan5). (11) 

At a driving frequency coinciding with resonance, co = cOq = fK/M, the tangent of the phase difference 



5 is inversely related to the quality factor, Q by 



1 ^ K ^ VkM ^ cOoM .q 

tan5 ^ ^ 

The quality factor also has the following common definitions 



( 12 ) 



Q= 



«u-"l 



CO,X 



(13) 



where col are the frequencies above and below resonance for which the average power is 3 dB down 

relative to its value at resonance and x is the decay modulus or characteristic time required for the free 
decay amplitude to decrease to 1/e of its initial value. The Q is also equal to the ratio of 2 tc times the 
mechanical energy stored in the oscillator to the energy dissipated per acoustic period.^ 



4.2 Quality Factor determination from pole-zero plots 

The Q, and thus complex modulus, can also be found from a pole zero plot of the mechanical 
admittance of the oscillator as will be illustrated in this section by measuring the frequency response of 
the output, eqn. (3), (induced EMF) of the detection transducer (which is proportional to velocity) to the 
input force, eqn. (1), (which is typically constant). 

Again considering the damped driven harmonic oscillator of frequency co, and noting that the complex 
velocity and displacement are related by u* =jco^*, the complex admittance u*/F* (which is equal to the 
inverse of the complex impedance) is given by 
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Y(jco) = 



( 14 ) 



J“ 



M 



/• R/- ^ K 

Oco) +-^CO)H 

M M 



The complex mechanical admittance is directly analogous to the the complex electrical admittance of 
an RLC circuit^ with the familiar electrical quandties voltage, current, electrical resistance, inductance, 
and capacitance (V, I, R, L, C) corresponding to the mechanical quantities force, velocity, mechanical 
resistance, mass, and compliance or inverse stiffness (F, u, R, M, 1/K) respectively. Letting jco = s, we 
can factor the admittance in pole-zero format as follows 



Y (s) = 



s 



M 



^ R . 

S H S-l- 

M 



IC 

M 



1 


s 


M 


_(s-y) (s- y*)_ 



(15) 



where Y,7* are complex conjugates. Noting cOq = VK/M, the characteristic roots are 




Thus the transfer admittance function has poles at y and y* and a single zero at co = 0. 



(16) 

Note that in this 



derivation we are only considering oscillatory solutions which correspond to R/(2M) < co^. The pole 
zero plot is illustrated in Figure 3 and the poles are located on a circle of radius cDq. 




Figure 3. Pole zero plot of complex admittance function. 



If we express the location of the poles as y = a + jb and y* = a-jb, then Q can be expressed as follows 

(17) 



-2a 
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and the resonance frequency is equal to 




where the first approximation is from the binomial expansion theorem and the second approximation is 
valid for low to moderate damping. 

6. MEASUREMENTS AND RESULTS 

A single optical fiber is used as the leg of a Michelson interferometer, and was cast into a long bar of 
E-CAST F-28 (1.24 cm diameter, 30.36 cm length, and density of 1 120 Kg/m^) at an off axis distances 
of approximately 6.0 mm. As the rod is sinusoidally driven with either longitudinal, transverse, or 
torsional stresses, strains will be produced in the leg of the interferometer which result in phase shifts 
relative to the unattached unstrained reference leg of the interferometer. 

The flexural, torsional, and longitudinal modes were clearly detected from the transducer coil output 
and interferometric phase modulated output. The first three torsional resonances occurred at 1485, 
3075, and 4684 Hz. The first two flexural resonances occurred at 182.3 and 514.7 Hz. The first three 
longitudinal resonances occurred at 2535, 5243, and 7712 Hz. These resonances were measured at an 
average room temperature of 24 °C. A typical interferometric output and electromagnetic coil output pair 
for the flexural fundamental resonance mode is illustrated in Figure 4. From eqn. (4-8), the complex 
shear and Young's modulus corresponding to these resonances are tabulated in Table 1 for the bar. 



INST TIMEl 





Fxd X 0.0 Sac 750fJ 



Figure 4. Typical interferometric output (Top) and corresponding transducer coil output (Bottom) at 
resonance for the fundamental flexural mode of the E-CAST bar. 
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Mode# 


Frequency 


Freq/Mode # 


Quality factor 


Loss factor 


Modulus 


[n] 


[fn] (Hz) 


[fn/n] (Hz) 


[Q] 


[5] (%) 


(Pa) 


1. Torsional 












1 


1485 


1485 


20.1 


4.97 




2 


3075 


1538 


19.8 


5.05 




3 


4684 


1561 


19.1 


5.24 


(shear modulus G*) 




Average 


1528 


19.7 


5.09 


0.935 X 10^ (1 -H j .0509) 


2. Flexural 












(3.0112)2 


182.3 


20.11 


30.0 


3.33 




(4.9994)2 


514.7 


20.59 


30.4 


3.29 


(Young's modulus Y*) 




Average 


20.35 


30.2 


3.31 


2.66 X 10^ (1 -t-j .0331) 


3. Longitudinal 










1 


2535 


2535 


21.6 


4.63 




2 


5243 


2622 








3 


7712 


2571 






(Young's modulus Y*) 


Average 




2576 


21.6 


4.63 


2.74 X 10^ (1 + j .0463) 



Table 1. Summary of the frequencies of the modes of vibration of the sample bar described in the text 



After using the Eff*3562 (dynamic signal spectrum analyzer) to sweep through the appropriate 
frequency range in order to excite each mode of vibration as illustrated in the measurement setup diagram 
in Figure 2, the resonances were determined from the displayed frequency response and verified by 
manually adjusting the oscillator for greater resolution. The quality factor, or Q's, for each mode were 
obtained by measuring the upper and lower 1/e (or -3dB) points on each side of resonance by adjusting 
the local oscillator. The Q's listed in Table 1 were obtained from this method. 

Since the HP3562 also possesses pole-zero curve fitting features, the resonances and quality factors 
were determined from the HP3562 generated curve fit parameters according to eqns. (17) and (18) 
derived earlier. It should be noted that the poles and zeros from the HP3562 generated plot are in Hertz 
and not in radians as suggested in eqns. (17) and (18). The frequency response of the detection wire 
coil transducer to the input drive for the torsional mode of vibration, the corresponding HP3562 curve fit 
response, and the list of poles and zeros are illustrated in Figure 5. The first four torsional resonances, 
which are harmonically related, are clearly identifiable. The frequency response of the flexural mode of 
vibration, the corresponding HP3562 curve fit response, and corresponding table of poles and zeros are 
illustrated in Figure 6. The corresponding resonance frequencies and quality factors determined from 
these poles and zeros are tabulated in Table 2. Due to temperature changes in the room, the measured 
resonance frequencies do differ for different measurements and thus more accurate temperature control 
by means of an environmental chamber is needed for measurement-to-theasurement repeatability. 

Having obtained the two independent moduli, the set of elastic constants for the isotropic material can 
be determined, i.e. the Poisson's ratio,v , is found to be 0.444 ffom4, v = [E/(2G)j - 1 , and using the 
average value for the Young's modulus of E = 2.7 x 10^ Pa and the shear modulus of G = 0.935 x 10^ Pa. 
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CURVE FIT FREQ RESP 




POLES S ZEROS lO 



1 


-37 . 944 


1^3 


1.4731k 


1 72 . 644 






2 


-70 . 1 lOS 


iLj 


3 . 07572k 


19 


. 2379k 






3 


-113. S72 


iLj 


4 . 69 109k 


1 . 


OB750kiij 


3 . 


. 57B5 Ik 


4 


- 162 .413 


±.j 


6 . 27439k 


- 1 . 


05 105k±_j 


3 . 


. 75323k 


5 








1 . 


0363Sk±_j 


6 . 


. 160 17k 


6 








- 1 . 


2 127 Iki^j 


6 . 


. 400 10k 



Figure 5. Frequency response of the transducer coil output for the torsional modes of the bar. The 
dotted line is a generated HP3562 frequency response with corresponding curve fit parameters (poles 
and zeros) listed below. 



FREQ RESP X-iea.5 Hz 

Yb—— 25.513 dB Ya—— 25.515 dB 




1 


-3 . 


, 50457 


±.j 


103.001 


15 1 


. 207 






2 


-9 . 


. 16000 


±.j 


515 . 069 


260 


. 023 


±J 


404 . 557 


3 










-67 . 


4909 


±j 


520 . 94 



Figure 6. Frequency response of the transducer coil output for first two flexural modes of the epoxy 
bar. The dotted line is a the HP3562 curve ft frequency response with corresponding curve fit 
parameters (poles and zeros) listed below. 
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Mode # 


Frequency 


Freq/Mode # 


Quality factor 


[n] 


[fn] (Hz) 


[fn/n] (Hz) 


[Q] 



Torsional 



1 


1474 


1474 


19.42 


2 


3077 


1539 


21.94 


3 


4692 


1564 


20.64 


4. 


6276 


1569 


19.32 




Average 


1536 


20.34 


Flexural 


(3.0112)2 


183.1 


20.11 


35.54 


(4.9994)2 


516.0 


20,59 


1SH4 




Average 


20.35 


36.86 



Table 2. Resonance frequencies and Q's determined from the pole zero 
plots in Figures 5 and 6 and the eqns. (17) and (18). 



6. CONCLUSIONS 

A resonant method for determining the complex Young's and shear modulus of elasticity has been 
presented that uses coils of wire to excite and detect vibrational modes of a bar. For detection, an optical 
interferometer, whose leg is cast in the material under investigation, can also be used to accurately detect 
the three resonant modes via strain induced phase modulations. 
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APPENDIX B. MATERIAL SPECIFICATION SHEETS 



1. URALITE3130 

2. PR1592 
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TECHNICAL DATA BULLETIN 

• No moisture sensitivity 

• Excellent hydrolytic stability {Reversion Resistance) 
• No TDI, no 4,4' •MethyIene*bis-(2*chloroaniline) 

• Fast cure, quick demolding 
• High abrasion resistance 
• Low Viscosity 

• Good Electrical Properties 

• Good Adhesion to many substrates 



DESCRIPTION 

Uralitd 3130 is a two component, natural amber or black urethane casting elastomer. This 
tough roctn temperature mixing and curing system has excellent handling characteristics. 
Uralite 3130 Is a middle of the hardness range, general purpose, versatile Uralite elastomer. 
It does not exhibit typical moisture sensitive characteristics of most urethane elastomers. 



USES: 



— Electrical and Electronic Encapsulating 

— Molds and Mold Facings 

— Metal forming pads 

— Flexible snakes 

— Gaskets 

— Abrasion resistant parts 






1 itb: 


Test Method 


Value 


(Metrics) 


Shore Hardness A/D 


ASTM D 224a68 


80-85/25-30 




Viscosity Part A — cps 


ASTM D 2393-71 


:^400 




Viscosity Part B — cps 


ASTM D 2393-71 


120 




Mixed Viscosity — cps 


ASTM D 2393-71 


2000 




Tensile Strength — psi (kg/cm’ ) 


ASTM D 412-68 


2750 


(193) 


Elongation — % 


ASTM D 412-68 


250 




Tear Strength — pli (kg/cm) 
Dielectric Strength — step 


ASTM D 624 Die C 


250 


(45) 


@77“F (25“C), volts/mil 
Dielectric Constant — @77"F (25* C) 
ICr Hz 
10^ Hz 

Volume Resistivity — @77* F (25*C) 


ASTM D 149-64 
ASTM D 150-54T 


240 

5.6 

7.2 




@1000V, ohm-cm 

Surface Resistivity — @77* F {25*0 


ASTM D 257-70 


1 X 10*^ 




1000V, ohms 


ASTM D 257-70 


2 X 10*^ 





(over) 
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URALITE® 3130 TECHNICAL DATA BULLETIN 



PROPERTIES: (Continued) 

t Insulation Resistance — @77® F (25®C) ohms 



Test Method 



Value (Metrics) 



after 28 days @ 95® F (35®C) 95% RH 


WE ATS612 


1 X JO** 




Pot Life - min. @ 77® F (25®C) 


ASTM D 2471-71 






Shrinkage • in/ln (mm/mm) 


ASTM D 2566-69 


0.0016 




Density 


ASTM D 792-66 






Cured Compound Ibs/in^ (g/cm^ ) 




0.039 


(1.07) 


Part A lbs. /gal. (g/cm^ ) 




8.58 


(1.03) 


Part B Ibs./gal. (g/cm^ ) 




9.15 


(1.10) 


Demolding Time — hrs @ 77® F (25®C) 




4 




Demolding Time — hrs @ 1 75® F (79® C) 




1 




Complete Cure — days @ 77® F (25® C) 




2-4 




Complete Cure — hrs. @ 1 75® F (79® C) 




2-3 




Color 




Amber or Black 




Ratio (By Weight); 








Part A 




100 




Part B 




30 




Ratio (By Volume): 








Part A 




100 




Part B 




28 




physical properties are representative of 


typical values obtained by tests conducted 


in the i 



ical Products Division laboratory. 



STORAGE: 

Uralite 3130 should be stored In a cool dry area. Avoid temperatures above 90® F and below 
65® F. Always blanket Uralite 3130 with dry nitorgen or 8440 Inert Blanketing Gas and reseal 
container after use. 



SURFACE PREPARATION SUGGESTIONS: 

Porous materials, such as plaster and wood, must have all surfaces that come in contact with 
Uralite 3130 well-sealed with a sealer which is compatible with urethane (acrylic sealers are 
suggested). After sufficient drying time, approximately thirty (30) minutes after last coat, the 
final surface preparation consists of the application of a release agent such as Partingkote® 
8302 (wiping off excess) to accomplish a complete and uniform release coating. 



MIXING: 

Ratio: Parts by Weight Part A — 100 Ratio: Parts by Volume Part A — 100 

Part B — 30 Part B — 28 

^ Weigh both components into same container and stir slowly for 2-4 minutes, scraping the sides 
^and bottom of container periodically to include unmixed material which may adhere to these 
surfaces. Care must be taken to avoid whipping air into mixture. Pour the thoroughly mixed 
Uralite 3130 onto the prepared surface and allow to cure. Working pot life of Uralite 3130 is 
approximately 14 minutes. 



URALITE® 3130 TECHNICAL DATA BULLETIN - SHEET #2 



CURING: 

Near ultimate physical properties are normally attained after 2 days at room temperature 
(77* F). Curing of Uralite 3130 may be accelerated by heating for 1-2 hours at 175*F. 
Demolding can t 

PACKAGING: 

UNIT 

DESIGNATION 

12 Qt. Pack 
(pre-weighed) 

Pail Pack 

Drum Pack 

CAUTION: 

' FOR INDUSTRIAL USE ONLY. 

This product contains an isocyanate based prepolymer, amines and heavy metal catalysts 
which are harmful if swallowed. It does not contain toluene diisocyanate or 4,4' methylene 
bis-(2-chloroaniline). It may cause burns or skin irritation. Use only in a well ventilated area. 
Protect skin and eyes from contact and avoid inhalation of vapors. 

Should skin contact occur, wash with soap and water. For eye contact, flush with water 
immediately and obtain medical attention. If swallowed, drink water, induce vomiting and 
contact a physician immediately. 



e accomplished after 4 hours at room temperature. 



PART A 

1 2-One Qt. Cans 



PART B 
12-!/? Pint Cans 



UNIT 

NET Weight 

22 lbs. 13.12 ozs. 



5 Gal. Pail 
55 Gal. Drum 



2/2 Gal. Can 
30 Gal. Drum 



52 lbs. 
585 lbs. 



WARRANTY The following is made in lieu of all warranties, express or implied; Seller's only obligation shall be to replace such 
quantity of this product which has proven to not substantially comply with the data presonted in the Manufacturer s latest bulletin 
describing the product. In the event of the discov«ry of a non-conforming product. Seller shall not be liable for any property loss or 
damage, direct or consequential, arising out of the use of or the inability to use the product Before using user shall determine the^ 
suitability of the product for his intended use. and user assumes all risks and liability whatsoever in connection therewith. 
Statements relating to possible use of our product are not guarantees that Such use is free of f>atent infringement or is approved by 
any government agency. The foregoing may not be changed except by an agreement signed by an officer of seller 

HEXCEL CORPORATION. RESIN CHEMICALS GROUP 
20701 NOROHOFF STREET. CHATSWORTH. CALIFORNIA 91311 • ^^^882.3022 

807 Pn'nittd m U.S A 
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AEROSPACE/ELECTRONIC PRODUCTS • - 

PR-1592 

POTTING AND MOLDING COMPOUND 

Use 

Especially formulated as a high hardness molding compound for electrical cables and/or as a polling com- 
pound for electrical connectors where resistance to cold flow, fexibility, high tensile strWigth, and exposure 
from -70®F to 275°F for extended periods are required. 




Description 

PR-1592 is a noncracking, chemically curing, polyurethane compound which is supplied in Amber or Black and 
in two-part quantities or premixed and frozen in plastic cartridges. PR-1592 cures at 180®F to a lough, flexible, 
cold flow-resistant, high tensile strength rubber with a Shore A hardness of 85. PR-1592 will cure at 75®F, but 
the resultant physical properties will be lower than when heat cured. Properly cured PR-1592 is designed to 
withstand temperatures as high as 275®F for extended periods without blowing, deterioration or loss in elec- 
trical properties, and has excellent flexibility at -70®F. 



The viscosity of PR-1592 allows the material to be degassed easily when supplied in two part units and pro- 
vides for excellent flow of characteristics in encapsulating, molding and potting applications. 



PR-1592 must be used with a primer on metal, neoprene and polyvinyl chloride surfaces. 



Specifications 

PR- 1592 meets the requirements of MIL-M-24041C Category B 



Application Properties (Typical) 



Color 
Part A 
Part B 

Mixing Ratio 



Nonvolatile Content 



Dark Amber or Black 
Straw 

Part AiPart B 

53:100 by weight 

99% 



Viscosity (Brookfield Spindle @ 10 rpm) 
Two-part unit 200 poises 

Premixed, frozen 700 poises 

Application Life 
Two-part unit 



Time to 2500 poises 
Premixed, frozen 


2 hrs. 




75 “F 




Time to 2500 poises 


1 hr. 




75“F 




Mold Release Time 


2 hrs. 




u_ 

0 

0 

00 




Cure Time 


To 75 Shore A hardness 


7 days 




75-F 


or 




6 hrs. 




180*F 




To 85 Shore A hardness 


21 days 




75°F 


or 


(Ultimate cure) 


16 hrs. 




180*F 


e: '^6 


NOTE: The above times are at the temperature indicated. 




Therefore, it is necessary to allow time tor the mass ot 
malerial. molds, etc. lo reach ihe temperalure. 





Performance Properties (Typical) 

(Cured 16 hours at 180®fD 



Color 

Specific Gravity 
Hardness, Shore A 
Volume Shrinkage 
Tensile Strength 
Ultimate Elongation 
Tear Strength (Die C)' 
100% Modulus 

j 

Compression. Set (ASTM 
22 hrs. ® 158®F • 



Dark Amber or Black 
1.08 
85 
40/0 

6000 psi 
425% 

320 lbs. /in. . :f 
600 .psi 

395. Method B) 

40% 



Ozone Resistance Conforms 

(Tested in accordance with MSFC-SPEC 202B) 

Moisture* Absorption 2.4% 

(Tested in accordance with MIL-M-24041) 

Flame Resistance * ' * ’ 

Current overload; 55 amps DC 
applied 2 V 2 mins. 

through #16 wire No ignition 

Flame exposure; ‘ 

ASTM D 635-63 Self-extinguishing 



(Continued on Page 2) 



supeasEOES 
INARCH 1987 


PRODUCTS RESEARCH & CHEMICAL CORPORATION 


Oi-T.0:0 


DATE ISSUED 


5454 SAN FERNANDO ROAD 
POST OFFICE BOX 1800 


410-416 JERSEY AVENUE 


PR. 1592 


SEPTEMBER 1907 


GLENDALE. CALIFORNIA 91209 
AREA CODE (818) 240-2060 


GLOUCESTER CITY. NEW JERSEY 08030 
AREA CODE (609) 456-5700 
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Performance Properties (Continued Irom Page }) 



with MIL-E-5272A) 



Electrical Properties (Typical) 

(Cured 16 hfs. @ 180®F) 



Hydlolytic Stability, Hardness Change 


Arc Resistance 


150 seconds 


After 120 days @ 
158®F, 950/0 RH 


- 150/0 


Dielectric Strength 


300 volts /mil 


Dielectric Constant 


• 


Adhesion, Peel Slength 




At 1 KHz @ 75®F 


6.5 


Aluminum 




At 1 MHz @ 75®F 


4.6 


(Primed with PR-420) 


40 piw 


Power Factor 




Neoprene 




At 1 KHz 75 «F 


0.08 


(Primed with PR-1523-M) 


25 piw 


At 1 MHz <g> 75®F 


0.06 


PVC (MEK tackified) 


28 piw 


Insulation Resistance 


Fungus Resistance 


Non-nutrient 


At 75®F \ 


5 X 105 megohms 


(Tested In accordance 




At 250»F 


150 megohms 



Resistivity At 75 At 250 

Volume, ohm-cm 1 x 10’^ 4 x 10^ 

Surface, ohms 1 x 10^^ 8 x 10’° 



NOTE: The above appficalion, performance and electrical property values are typical lor the material, but are not Intended for use in 
specilicalions or lor acceptance Inspection criteria because ol variations in testing methods, conditions and configurations. 



Purchasing Data 

PRODUCT DESIGNATION PRIMER DESIGNATION 



When ordering this product, designate PR-1592 and 
color as follows: PR-1592 Amber or PR-1592 Black. 

NOTE: Reler to “Surface Preparalion” for primer requiremenl. 



STANDARD PACKAGING 






Two-Part Units 


Part A 


Pan B 


Units per 


Designation 


Container 


Container 


Case 


40 fl. OZ. unit 


1-pt. can 


1-qt. can 


9 


160 fl. oz. unit 


V 2 -gal. can 


1-gal. can 


4 



NOTE; The unit designates the total fluid ounce content ol Part 
A and Pari B (128 fluid ounces per gallon). Standard units are 
furnished with a premeasured quantity of Part A and Part B, in- 
dividually packaged. 



When ordering primer, designate PR-420, PR-1523-M 
or PR-1543 as required under "Surface Preparation." 



STANDARD PACKAGING 
PR-420 is packaged as follows: 

Part B Total 

Designation Container Contents 



# 7 kit V 2 -pt. can 7 fl. 02 . 

/M4 kit 1-pt. can 14 fl. oz. 



No. per 
Case 

16 

16 



NOTE: the kit number designates the total fluid ounce content of 
Part A and Part B. Kits are furnished with a premeasured quantity 
of Part A and Part B individually packaged. The kits are designed 
so that adequate space is available in the Part B container for ad- 
dition of Part A and mixing. 



Frozen Cartridges 


Approximate 


No. Cartridges 


PR-1523-M and PR-1543 are packaged ready to use 
in the following containers: 

^ No. per 


Designation 


Contents 


per Case 


Designation 


Container 


Contents 


Case 


2 V 2 OZ. cartridge 


2V- fl. OZ. 


72 


2 fl. OZ. 


2 fl. OZ. bottle 


2 fl. OZ. 


6 


6 oz. cartridge 


53/4 fl. OZ. 


36 


V 2 pint 


V 2 -pint bottle 


8 fl. OZ. j 


12 


12 oz. cartridge 


IIV 2 fl. OZ. 


36 


One pint 

1 


1-pint bottle 


16 fl. oz: ' 
# 


12 



Surface Preparation 

CLEANING 

Connectors or other metal surfaces must be free of grease, oil and wax In order to insure good adhesion. Use 
oil-free solvent applied with a small brush or oil-free cloths for cleaning (reclaimed solvents should not be 
used). Premixed cleaners are commercially available. Do not expose wire insulation and inserts to the cleaning 
solvent beyond the time necessary for adequate cleaning. ^ ^ * 



APPLICATION OF PRIMER 

For maximum adhesive strength between PR-1592 and the material to which it is to be bonded, the following 
surface preparations are required: • ' 

NOTE: Do not dip priming brush Into primer supply. To maintain an unconlaminaled primer supply, pour a small portion of primer into a 
clean container from which it should be used. Reseaf primer supply immediately alter portion has been removed. 

(a) Metal 

Metal must be primed with PR-420. Thoroughly mix 1 part of Part A with 6 parts of Part 8 by volume. Do not 
mix more than can be used within a 4 hour period. Brush a thin film of mixed PR-420 on all Inside surfaces of 
connectors and on wire, but not on the insulation. Let primer dry for 1 hour at 75®F. If primer becomes con- 
taminated, reclean primed surface lightly with methyl ethyl ketone and dry. Stripping the primer from the con- 
nector and repriming is not necessary. 
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(b) Neoprene 

To obtain good adhesion to neoprene insulation, the surface should be abraded with. a suitable abrasive to • 
remove grease, oil, wax or mold release. Remove rubber particles with a dry oil-free brush. 

Apply a liberal coat of PR-1523-M to the clean neoprene surface by brush and allow to dry for approximatefy 
30 minutes at room temperature. After 30 minutes drying time, wipe off excess PR-1523-M primer with a clean, 
gauze pad and start the potting or molding procedure. Drying time of PR-1523-M should not exceed 4 hours at 
room temperature before potting or molding. If primed surface becomes contaminated or potting or molding is 
not accomplished within 4 hours after application of PR-1523-M, buff neoprene and repeat priming procedure. 
NOTE; PR-1523-M Is hygroscopic and musi be kepi free of moislure. When PR-1523-M hydrolyzes, a dark grainy precipilale is formed 
decreasing Ihe primer usefulness. Malenal containing precipilale should be lesled lo determine lhal adhesion is salisfactory before using. 

(c) Polyvinyl Chloride 

To obtain good adhesion to polyvinyl chloride insulation, the surface should be made tacky with methyl ethyl 
ketone. The use of a primer may be necessary only with some formulations of polyvinyl chloride. Therefore, it 
is suggested that tests be made to determine the adhesive strength of PR-1592 to the polyvinyl chloride in 
question. 

Should a primer be required, then apply a thin coat of PR-1543 lo the lackified surface by brush and allow to 
dry 30 minutes at room temperature. If primed surfaces become contaminated before potting or molding, buff 
primed surface with a suitable abrasive and reapply a thin coat of PR-1543. 

It should be noted that there are many formulations of polyvinyl chloride. Therefore, it is suggested that before 
production quantities of PR-1543 are ordered, tests be made to determine the adhesive strength of the 
PR-1 543/PR-1 592 system to the polyvinyl chloride in question. 

NOTE: PR-1543 is hygroscopic and musI be kepi free of moislure. When PR-1543 hydrolyzes, a precipilale is formed decreasing the 
primer usefulness. 

(d) Teflon* and Other Fluorocarbons 

To obtain good adhesion to insulation made of Teflon and other fluorocarbon resins, it is essential that the in- 
sulation be etched or treated to provide a bondable surface. After neutralization of the etchant, in accordance 
with the manufacturer’s instructions, apply PR-1592 directly to the etched surface without primer. 

REPAIR 

To obtain good adhesion to previously cured PR-1592, the surface should be buffed with a suitable abrasive to 
remove grease, oil, wax or mold release. Remove rubber particles with a dry, oil-free brush. 

CAUTION: Do not use solvents for cleaning cured PR-1592. Apply new PR-1592 directly to buffed surface and 
cure as recommended. No primer is required. 

Mixing Instructions 

FOR TWO-PART UNITS 

Do not open containers until ready to use. 

Part A will solidify at room temperature. Prior to use, loosen lid and warm to 250® ± 10®F with thorough stirring. Do 
not heat over 260®F. When warming the material, use a thermometer to determine the actual material temperature. 
Liquefaction is complete when the material becomes smooth and uniform in appearance and loses all signs of 
graininess. Stirring is essential during liquefaction to provide a uniform material and to hasten melting. Care should 
be taken to dissolve all solidified Part A around the top of the container. Trace quantities of unliquefied Part A will 
cause premature solidification. Do not store Part A at temperatures exceeding 100®. 

Premeasured quantities may be divided into smaller quantities by using the following proportions: 

By weight 53 Parts A to 100 Parts B. 

NOTE: After removing a portion of Part B from a full container, moislure in Ihe air in Ihe empty portion of Ihe Part B container lends lo 
cause Ihe remaining material lo skin over during extended periods of standing. This material may be used by removing the skin. 

Where a dense compound free of voids is required, it is recommended that the mixed material be degassed 
before applications are made. Standard vacuum equipment may be used or, for small usages, the material may 
be degassed in a standard laboratory desiccator connected lo a vacuum pump. 

The following degassing procedure is recommended: 

1. Stabilize Part A and Part B at 75® to 85®F before mixing. 

2. Place Parts A and B in a clean, dry, metal container having at least twice the volume of the material to be 
degassed. Mix Part A and Part B thoroughly with a metal mixing paddle. NOTE: DO NOT USE WOOD. 

3. When Parts A and B are mixed thoroughly, degas the mixture until foaming subsides which is approximately 
10 minutes at a pressure of less than 5mm of mercury for a pint of material. Larger quantities will require 
slightly longer periods of degassing. 

4. When the material is to be applied by extrusion gun, it is suggested that after tranferring the degassed 
material into the extrusion gun cartridge the material be degassed at a pressure of less than 5mm of mercury 

'Registered trademark tor DuPont lelrafluoroelhylene resin. 
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for 2 minutes to remove any air which may have been entrapped during cartridge filling. If extra care is 
laken when filling cartridges with degassed material, such as flowing the material down the side of the car-^- 
tridge^ the material may not need to be degassed in the cartridge. 

NOTE; After mixing Paris A and 8? subsequeni operations should be accomplished as quickly as possible to minimize the reduction ol ap- 
plication life. 

After units of PR-1592 have been mixed and degassed, they may be frozen for storage under refrigeration. Use 
of a quick-freeze technique is recommended so as to minimize the amount of applicalion life that would be lost 
in a slower cooling procedure. One successful method is to immerse the filled cartridge in a slurry of dry ice 
and alcohol for 10 minutes. The temperature of the sealant will drop to approximately -70®F and the car- 
tridges may be transferred to a refrigeration storage unit maintained at -20®F or below. 

NOTE: Approximately 1 hour application life will be lost in the freezing and thawing process. 

It is suggested that tests be conducted to determine whether degassing is necessary^for the material to meet 
particular requirements. When it has been found that degassing is not necessary to meet requirements, mix 
entire contents of Part A with Part B, or use in the proportions as shown above, and thoroughly mix. Slow mix- 
ing by hand or with a mechanical mixer is recommended. A high-speed mechanical mixer will generate internal 
heat which reduces application life and will whip air into the mixture resulting in a porous material upon curing. 

THAWING OF FROZEN CARTRIDGES 

Frozen cartridges of PR-1592 must be thawed as follows: 

NOTE: It is essential that the thawing lime and temperature be controlled closely to obtain the maximum applicalion lile in the shortest 
thawing period. An increase in eilher thawing time or temperature will result in an incomplete thaw. 

Remove cartridge from storage and thaw for 30 minutes in a 100®±5®F heating block or water bath. If an 
oven is used lo thaw the material instead of the recommended water bath or heating block, the time and 
temperature will have to be determined because of the variations in oven controls and circulation. 

NOTE: If a water bath is used, be sure water does not enter cartridge. The recommended method ol preventing water from entering the 
cartridge is to use a metal sleeve, closed at one end, with a diameter just large enough to allow the cartridge to be inserted to the lip and 
then place the sleeve upright in the water bath. 

CAUTION: When thawing frozen cartridges of PR-1592, keep cartridges in an upright position, nozzle end down 
v/ith CaPlug in place, to prevent air thal may have been drawn into the cartridge during freezing from entering 
and being trapped in the compound. After thawing is complete, just before placing the cartridge in the gun, in- 
sert a thin piece of metal or spatula blade between plunger and cartridge wall to the shoulder of plunger. Then 
force plunger down to exhaust any air between compound and plunger face and continue until material just 
starts to extrude between plunger and cartridge wall. Remove metal or spatula while holding plunger down. 

Application Instructions 

rOR POTTING 

Application is accomplished by injecting the sealanl into the end bell of the connector which has been coaled 
previously with primer. Separate the wires evenly so that the potting compound will flow around all wires and 
soldered connections. 

For most potting applications, a Semco #440 nozzle with 1/8" diameter has been found mosi practical. How- 
ever, nozzles of smaller diameter are available. j 

When applying PR-1592 with an extrusion gun, the nozzle first should be inserted between the center wires so 
that the sealant will flow around the contacts. It may be necessary to reposition the nozzle for even coverage. 
Care should be taken to eliminate voids and air bubbles. It is good practice when filling the connector to keep 
the tip of the nozzle just above the swell of the material level, moving the nozzle up as the connector fills. The 
tip of the nozzle should not be held so as to allow it to become submerged nor should it be held so as to allow 
the material to fold. 

!n applying PR-1592, every effort should be made to eliminate voids and air bubbles. A small piece of metal, 
such as a welding rod, should be used to tamp or pack the compound around the wires and the base of the 
pins. In conjunction with the tamping, the connector should be tapped on a resilient surface or even * - ’ 
mechanically vibrated to facilitate the flow of the sealant into the small recesses. 

To reduce viscosity and enhance flow of compound, connectors or molds may be healed up to 140®F prior to 
filling. Do not heat the molding compound, as it will thicken and shorten application life. 

Molds may be prepared to the exact shape and form required in a particular connector application by utilizing 
various plastic materials commercially available. A form or mold may be made for the back of a connector by 
using masking or polyethylene tape. Reference Is made to the handbook. Installation Practices for Aircraft Elec- 
tric and Electronic Wiring, (Navy) NAVAER 01-1A-505 or (USAF) T.O. 1-1A-14 dated 30 June 1975. 

FOR MOLDING 

When PR-1592 is used as a molding compound, molds should be provided for forming the particular form and 
shape required. To facilitate the removal of PR-1592 from the mold after curing, spread a thin film of mold rer 
lease compound on mold surface to prevent adhesion. Metal molds may be coated with Teflon tetrafluoroethylene 
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resin to provide a permanent release film and eliminate the necessity of using a mold release compouQd. 

Inject PR-1592 into the mold, which has been coated previously with mold release and/or on previously 
prepared surfaces, and allow to firmly set before removing mold. It is recommended that injection holes be 
located in the bottom of the mold and air bleed holes located in the top of the mold to prevent air pockets in 
the mold. Flash may be trimmed with a sharp knife or razor blade. 



Mold Release and Cure 

The mold release time for PR-1592 is dependent upon the temperature, quantity of material, mold mass, ect. In 
molds containing less than 24 fluid ounces of PR-1592, the mold release time is approximately 2 hours at 
180®F. In molds containing 24 to 48 fluid ounces of PR-1592, the mold release time is approximately 3 hours 
at 180°F. Longer times will be required for the mold release of PR-1592 in larger molds. 

\ 

PR-1592 cures to a tough, resilient material having a 75 Shore A hardness in 6 tiours at 180®F. PR-1592 will 
cure at 75®F, but the physical properties will be lower than when cured at 180®F. 

Cleaning of Equipment 

Wash equipment, tools and brushes with methyl ethyl ketone immediately after use or before material cures. 
Use commercial stripping compounds to remove cured material. 

Safety Precautions 

PR-420 contains flammable and volatile solvents. Keep away from heat, sparks and flame. Proper precautions 
for working with flammable liquids should be followed as well as applicable safety precautions. 

storage Life 

Storage life of PR-1592 in two-part quantities is 12 months when stored at temperatures below 80'’F in original 
unopened containers. 

Storage life of premixed, frozen PR-1592 is at least 7 days when stored at -20®F or at least 30 days when 
stored at -40®F. 

Storage life of PR-420 is approximately 1 year when stored in original unopened containers below 80®F. 

Storage life of PR-1523-M and PR-1543 is approximately 6 months when stored in original unopened containers 
below 80®F. 

NOTE: PR-1523-M and PR-1543 are hygroscopic and must be kept free of moisture. When these primers hydrolyze, a precipitate is lormed 
decreasing the primers usefulness. Material containing precipitate should be tested to determine thal primer is satisfactory before using. 



Health Precautions 

Part B of PR-1592 and related primers, PR-420, PR-1523-M, and PR-1543, contain isocyanates and can pro- 
duce irritation or allergic reaction following contact with the skin, eyes or mucous membranes. Avoid all contact 
with the uncured components of these materials. In case of contact, promptly wash off with copious 
quantities of water and follow with a soap and water wash. Avoid breathing vapors. Use with adequate ventila- 
tion. Individuals with chronic respiratory problems or prior respiratory reactions to isocyanates shoulpj not be ex- 
posed to vapors. 

PR-420 contains mixed solvents and a lead compound. Use adequate ventilation or air-supplied respirator 
during application. 

In cases of extreme exposure or adverse reactions to any of the products mentioned above, remove affected personnel 
to fresh air and obtain immediate medical attention. If a rash develops, consult a physician. Ordinary hygienic principles 
such as washing hands before eating or smoking must be observed. For complete health and safety information, 
consult a Material Safety Data Sheet. * * ’ 



"PRC” Is a trademark of Products Research & Chemical Corporation, registered with the U.S. Patent Office 

All recommcrxjalions. stalemenis and lechntcal data contained herein are based on lesls we believe lo be reliable aixJ correct bul accuracy and completeness 
of said tests are not guaranteed and are nol lo be conslruded as a warranty, either expressed or implied User shall rely on his owninlomiation and lesls lo deter- 
mine suitability of the product for the intended use arxJ user assumes all nsk and liability resulting Irom his use ol the prtxJuct Seller's and manulacturer's sole 
responsibility shall be lo replace that portion of Ihe product ol this manulacturer whch proves lo be deleclive Neither seller nor manulacturer shall be liable lo the 
buyer or any third person for any injury, loss or damage directly or indirectly resulting Irom use of, or inability lo use, the product Recommendalior« or statements 
other than those conlamed in a written agreement signed by an officer of the manulacturer shall nol be binding upon Ihe manulacturer or seljer. 

Printed in U.SA. 
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APPENDIX C. ENVIRONMENTAL CHAMBER 
SPECIFICATION SHEETS 
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TE^yiPERATU RE/HUM I DiTY 
Ch \MBERS A 



ASSOCIATED ENVIRONMENTAL SYSTEMS 



These are all standard untft — readily available - many Irom stock. Now you can 
pck the unit to suit your needs and your budget . Bench Top and Floor Models — four 
temperature ranges to choose from — choice oi 6 working volumes — 15 different units in all 



BENCH TOP MODELS— BHD Series 

The 8H0-403. 405 or 408 (last two digits indicate cubic foot working volume) are bench lop 
units that will allow testing to meet known Mil Std Humidity Speciiicauons These units are 
basically Humidity chambers since they have a limited temperature range when compared to 
the Floor Models in order to properly control humidity a refrigeration system is required 

Simple and easy to operate The basic chambers have a digital dry bulb and a direct setting 
digital % RH controller (see optional extras for automatic programming and recording) Simply 
connect the power, 1/4" water supply line and 3/8" dram line and you are ready to operate 

Temperature Range; - 18X (0‘’F) to + 93"C ( + 200‘’F) 

Humidity Range; 207o to 95% ± 5% RH within the range of + 20'*C (68“F) to + 85“C (IRS^F) 
and IS limited by a +4®C (39*’F) dew pomt temperature, using standard controls Optional 
microprocessor and recorders are available Lower humidities available v/ith optional chemical 
drier 



FLOOR MODELS— HD SERIES 

This series of Humidity Chambers allows you to custom design a chamber to satisly your needs 
and your budget. Four working volumes to choose from 08. 16, 33 and 64 fi^—Three Tem- 
perature Ranges. 



Humidity Range: 20% to 95% RH within the range oi +20® C (68® F) to +85® C (185® F) and is 
limited by a +4°C (39® F) dew point temperature 

With 86651 Programmable controller. Control statwlily of ±3% RH. 



HOW TO ORDER BY MODEL NO. 



Example; MODEL HO • 264 



HO Series — 
Floor model 



200 Series 

(-70®Cto +177®C) 

64 cubic feet 



S C 

1 C = ”C 

f = ”F 

S = standard 

voltage & 
frequency. 

0«= optional 
nonstandard 
power 

requirements 



5e/ecf the temperature range: 

The 400 series has a temperature range of ~ 18®C to + 177®C (0®F to +- 350®F) 

The 500 series has a temperature range of - 4(TC to + 177®C ( - 4Q®F to +• 350‘’F) 

The 200 series has a temperature range of -70®C to +177®C (-100® F to +350®F) 

Select the working volume: 

Indicate it as the last two digits of the model number - Sizes; 08. J6. 33 or 64 ft.^. 



FEATURES (illustrated on page 4) 

• FAILSAFE- An ad|ustable high temperature failsafe 
IS standard on all Associated Environmental System 
Chambers It is lactory set at the ultimate high tern 
perat'jre capability of the chamber Upon receipt of 
the chamber, the customer should set the dry bulb 
controller to 10® above the highest test temperature 
and adjust the -failsafe. This will prevent the test 
chamber from a temperature overrun 

• FUSIBLE LINKS (in the heater circuit) - For over 
temperature protection 

• MULTIPANE VIEWING WINDOWS - In the door ol 
all chambers (optional on BHO Senes) 

• INTERNAL LIGHT - with external switch standard 
on all floor models optional on bench top units. 

• PILOT LIGHTS - to monitor proper function of vari- 
ous systems. 

• INTERIORS - All stainless steel type 304 heharc 
welded 

• EXTERIORS - Heavy gauge, cold rolled steel with 
two coats of textured enamel finish 

• INSULATION — High density, low K lactor, non-settling 
fiberglass. 

• VAPOR GENERATORS - When humidity is required 
the controller signals the vapor generator to induce 
water vapor into the chamber This method ol creat- 
ing high humidities is superior to the older methods 
of spraying water or passing warm air over water. 

• VAPOR GENERATOR FAILSAFE - All A E.S humid- 
ity chambers are equipped with a factory set failsafe 
control to shut down the vapor generator in the event 
of a no water situation 

• FREEZE PROTECTION - Since all AES Chambers 
have the capability of going below the freezing point 
of water, each humidity chamber has a freeze protec- 
tor which is factory set at 1 67®C ( +- 35®F) This con- 
trol actuates a solenoid valve to dram water/rom the 
chamber and prevent damage due to freezing of 
water 

• FORCED AIR CIRCULATION - High volume air cir- 
culation IS standard on all units to achieve uniformity 
in the work space The interior fan is connected to an 
external motor with a stainless steel shaft 

• WIRING - Coded readily accessible Irr accoMance 
with NEC 

• CONTROLS - All located in a hinged access panel 
with all relays and fuses readily accessible 

• ACCESS PORT - A 2" diameter access port and plug 
in left side wall on floor models and right side wall 
on bench tops is pre-punched in the metal - addi- 
tional ports and plugs available. 

• DIGITAL TEMPERATURE & HUMIDITY CONTROL- 
LERS 
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BENCH TOP MODELS 
BHO-403. 405. 408 



FLOOR MOOeiS 
HO-400, 500. 200 



UNIT SPECIFICATIONS 





BENCH TOP MODELS 
Temp. Range; - 1TC to - 93"C 
Hum. Range; ZOH lo 95W 


Temp Range: - 18^C to ♦ 177'C 
Hum. Range.* Vfh lo 95*/% 


FLOOR MODELS 
Temp Range: -40 Clo » 17TC 
Hum. Range: 20*>V lo 95W 


Temp. Range; - 73*C lo ▼ 177*C 
Hum. Range: 205V lo 955V 


MODEL NO. 


BHD-403 


BHD 405 


BHD-408 


HD 400 SERIES 


HD-SOO SERIES 


HD-200 SERIES 


Working 


ft J 


03 


05 


08 


08 


16 


33 


64 


08 


16 


33 


64 


08 


16 


33 


64 


Volume. 


Liters 


85 


14? 


227 


227 


453 


765 


1812 


106 


227 


453 


765 


1812 


227 


453 


765 


1812 


Power 


VAC* 


115 \ 100 


115 1 100 


230 I 200 


230 1 200 


230 [ 230 


230 1 230 


230 1 230 


230, 


1 230 


230 1 200 


230 1 200 


230 1 200 


230 1 200 


230 1 200 


230 1 200 


230 \ 380 


Required; 




t 


t 


1 


1 


1 


1 


1 


1 


1 


1 220 


1 220 


1 230 


• 220 


1 220 


1 220 


1 400 




0 


1 ' 1 


1 1 1 


1 ' 1 


1 ' 1 


1 ' 1 


I ' 1 


1 • 1 


1 ' 


' 1 


1 1 3 


3 ; 3 


3 ; 3 


' 1 ^ 


3 ' 3 




^ 1 ^ 




HZ 


60 I 50 


60 ! 50 


60 1 50 


60 1 50 


60 \ 50 


60 50 


60 I 50 


60 ; 


1 50 


60 , 50 


60 , 50 


60 , 50 


60 ; 50 


60 I 50 


60 \ 50 


60 ; 50 




Amps 


16 1 16 


16 1 16 


10 1 10 


20 1 24 


28 1 28 


30 1 30 


3? \ 3? 


26 


26 


30 1 30 


30 1 30 


35 1 35 


35 f 30 


40 1 35 


50 1 50 


64 1 34 


Shipping 


Lbs. 


400 


528 


600 


650 750 


900 


1300 


1800 


750 


850 


1000 


1400 


1900 


1050 


1200 


1725 


2375 


Weiqhl 


Xq. 


18? 


240 


273 


295 340 


408 


590 


816 


340 


386 


454 


635 


862 


476 


545 


184 


1080 


Window in 


Door; 




































Inches 


Optional 


Optional 


Optional 


8«a 


12 * 1? 


1? I 1? 


1? * 12 


8 1 


i8 


12 X 12 


12 X 12 


i? X 12 


8x8 


I2x 12 


12 X 12 


12 X 12 


Access Port' 


7' dia 151 mml port & plug center 










r dia 


151 mml port & plug lelt side wall 












right side wall 




























Nominal 




































Healer 


Xw 


0 7 


07 


10 


25 


40 


50 


60 


25 


40 


50 


60 


25 


40 


50 


60 


Refiiqeraiion 


































System. 


HP 


1/3 


1/3 


1/2 


3/4 


1 


1 


VA 


\'A 


2 


3 


5 


2x 2 


2 X 2 


5 X 5 


7Vj X 7W 


Temp. Rise Time ■ • 


































Ambieoi 10 Upper 


30 


30 


40 


40 


40 


60 


30 


40 


40 


60 


30 


40 


40 


60 


30 


Limil m Min 


































Temp Pull Down • 


































Time: Ambient to 


45 


45 


45 


60 


60 


60 


60 


60 


60 


60 


60 


90 


60 


120 


120 


Lower Limit m Mm 


































Programmable 

Controller 


Opuonal 


Optional 


Opilonal 






#86651 


PROGRAMMABLE CONTROLLER STANDARD 







•Rise & PuM Do^in Times may vary decwvJmg on amb<en( condihons and pcw*ef suopty * VAC leounes additional transformer 





BHD-403 


BHD-405 


BHO-408 


HD-400. 500. 200 SERIES 


Working Volume; 
fl.7 


03 


05 


08 


08 


’ 16 


33 


64 


Size H X W X 0 In. 
H X W X D cm. 


t6 X 20 X 16 

41 X 51 X 41 


20 X 22 X 20 

51 X 56x 5t 


24 X 24 X 24 
61 X 61 X 6t 


24 X 24 X 24 
61 X 6t X 61 


30 X 30 « 31 
76 X 76 X 79 


36 X 36 X 45 
92 « 92 * t14 


48 X 48 X 48 
122 X 122 X 122 



^ ^ i *ni y ns i© -i ^ 



OPTIONAL EXTRAS 

1. Programmable Controller 86651 (standard on all floor models) 

2. 86649 IEEE-488 interface ( required #1 option) 

3. Product Protector limit controls 

4. Honeywell 2 Pen 12* Circular Chart Recorder 
.5. 8* X 8* Window in Door for BHO-400 models 

6. Manual Window Wiper 

7. internal Light with External Switch (standard on ail floor models) 

8. Water Storage Tank External 

9. Water Re-cyding System (irK:ludes external storage tank) 

10. Water Oemineralizer System 

11. ' Extra Oemineralizer Cartridge 

12. Fully Adjustable Shelves 

13. Additional Shelves 



14 Casters 

15. Additional Ports— Specify Location 
2*. 3*. 4*. 5*. and 6* diameter 
Id Chamber Cart (for BHO-403) 

17. Floor Stand (BHO-405 & 408) 

18. Liquid CO? Cooling Capability— specify psi 

19. Liquid N2 Cooling Capability 

20. Notse Reduction Package lor floor models 

21. Water Cooled Condenser 

22. Other Power Source Requirements 

23. Desiccant Oehumidifier lor low % RH requirements (consult factory 
with specs) 
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FEATURES 



FAILSAFE 

CONTROL PANEL 

FREEZE PROTECTION 

VAPOR GENERATOR 
FAILSAFE 




ACCESS POR^- 
AND PLUG 



A 




VIEWING WINDOW 



EXTERIOR HAS 
TEXTURED 
ENAMEL FINISH 
INTERIOR IS 
STAINLESS STEEL 



WIRING IS COOED 
AND READILY ACCESSIBLE 



AIR COOLED 

REFRIGERATION. 

SYSTEM 

WATER COOLED 
(optional) 






BH 

403 


0-40Q Sei 
405 


ics 

408 


HC 

08 


1-400. 500 
16 


or 200 Sci 
33 


ies 

64 


Dimensions; 


in/cm 


in/cm 


in/cm 


in/cm 


in/cm 


in/cm 


in/cm 


Outside Height 


24/61 


28/71 


32/01 


78/198 


78/190 


84/213 


99/251 


Outside Width 


42/107 


50/127 


55/140 


46/117 


53/135 


50/147 


70/170 


Outside Depth 


24/61 


30/76 


36/91 


52/132 


60/152 


72/103 


80/203 



In line with our corporate policy of continuing product improvement and quality control. Associated reserves the right to change product witho 
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APPENDIX D. RESONANCE TRACKING SYSTEM PROGRAM 



LISTING 
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10! 


RT20C 


: Resonance Tracking System 


20! 


Written by 


: Brian Beaton 


30! 


Date 


: 18 Mar 90 


40! 


Last modified 


: 21 FEB 91 


50! 


Features added 


: 1 .Read one more DVM for Inphmag 


60! 




BDAT change to 40 records 


70! 




Measurement time 20 hours 


80! 




Layout of printer output change 


90! 




Format of data storage change(for Mac ) 


100! 




28 Jan 91 by Tan B.H. 


110! 




2. Display In phase mag vs Time on CRT 


120! 




Add plots of Inphmag vs Freq and 


130! 




Inphmag vs Temp 


140! 




13 Feb 91 by Tan B.H. 


150! 




3. Compute reduced fi^equency 


160! 




Convert old data of RT20B to RT20C format 


170! 




BDAT change to 50 records 


180! 

190! 




21 FEB 91 






210! Program initialization 

2201 ******************************************************************** 



230! 
240 
250 
260! 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 ! 
400 
410 
420 ! 
430 
440 
450 
460 
470 ! 
480 
490 
500 
510 



MASS STORAGE IS ":INTERNAL,4,0" 
RE-STORE ”RT20B:INTERNAL,4,r' 



OPTION BASE 0 

DIM Theimistor(0:240)Jlesfi:eq(0:240),Time(0:240),Inphmag(0:240) 
DIMMod(0:240)^(0:240),Y(0:240),Alphat(0:240),Redfreq(0:240) 
DIM Temp(0:240) 

DIMGtor(0:240),Elong(0:240),Eflex(0:240) 

DIM Label$(0:25)[50] 

DIM Main_title$[50],Sub_title$[50] 

DIM X_axis_name$[50],Y_axis_name$[50] 

DIMBlockl$[50]31ock2$[50],Block3$[50],Block4$[50] 

DIMCI$[25],C2$[25],C3$[25],C4a$[25],C4b$[25],C5$[25] 

DIMC6$[25],C7$[25],C8$[25],C9$[25],C10$[25] 

DIMLI$[50],L2$[50],L3$[50],L4$[50] 

Clear$=CHR$(255)&CHR$(75) ! CLEAR SCRN key 

Home$=CHR$(255)&CHR$(84) ! HOME key 

Flag=0 !Flag=I to indicate editing data 

Tmode=l ! first torsional mode 

Lmode=l ! first longitudinal mode 

Fmode=3.01 12 ! first flexural mode 



Sample_rate=5 ! 1 sample every 5 minutes 

Argl=60/Sample_rate ! samples per hour 

Pi=PI ! 3.14159265 

Default_grid$="Partial" ! Full, Partial or No 
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520! 

530 Label$(l)="Resonance Frequency vs Time" I Graph Main Titles, 
540 Label$(2)="Temperature was varied" ! Sub Titles, Axes names 
550 Label$(3)="Time, hours" 

560 Label$(4)="Resonance frequency, Hz" 

570 Label$(5)="Temperature, C" 

580 Label$(6)="Resonance Frequency vs Temperature" 

590 Label$(7)="Temperature vs Time" 

600 Label$(8)="Frequency variation was observed" 

610 Label$(9)=" Y oung's Modulus vs T emperature" 

620 Label$( 10)="Shear Modulus vs Temperature" 

630 Label$( 1 1 )="Shear modulus. Pa, x 1 0'^" 

640 Label$( 12)=" Young’s modulus. Pa, xlO'^" 

650 Label$(13)="In Phase Magnitude, DCV" 

660 Label$( 14)="In Phase Mag vs Time" 

670 Label$( 1 5)="In Phase Mag vs Temp" 

680 Label$(16)="In Phase Mag vs Freq" 

690! 



700 


Penl=l 


! White 


Default colors 


710 


Pen2=2 


! Red 




720 


Pen3=3 


! Yellow 




730 


Pen4=4 


! Green 




740 


Pen5=5 


! Cyan 


(greenish blue) 


750 


Pen6=5 


! Blue 





760! 

770!******************************************************************** 



780 Options: ! Allow user to process existing data or collect new data 
790 !******************************************************************** 
800! 

810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 



! Clear the CRT 
! Home display 
! Get screen width 
! Leading spaces for centering 

! Use CRT for displaying menu 
! Start at top with blank line 
Purpose" 

Process existing data" 

Collect new data" 

Editing existing data" 

! Off all keys 
GOSUB Comment 1 



OUTPUT KBD;Clear$; 

OUTPUT KBD;Home$; 

STATUS l,9;Screen 
Center=(Screen-28)/2 
GCLEAR 
PRINTER IS CRT 
PRINT TABXY(U) 

PRINT TAB(Center);"Key 
PRINT TAB(Center);"— 

PRINT TAB(Center);"0 
PRINTTAB(Center);"4 
PRINT TAB(Center);" 5 
FOR Keynumber=0 TO 9 
ON KEY Keynumber LABEL ' 

NEXT Keynumber 

ON KEY 0 LABEL "Existing Data" GOTO Existing_data ! Turn on Key 0 
ON KEY 4 LABEL "New Data" GOTO New_data ! Turn on Key 4 
ON KEY 5 LABEL "Editing Data" GOTO Editing_data ! Turn on Key 5 
990Blink0: WATTl 
1000 DISP "Select an option" 

1010 WATTl 

1020 DISP 
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1030 GOTO BlinkO 

1040 RETURN 

1050 !******************************************************************* 
1060 New_data; ! Main program driver for collecting new data 

1070 [************ ******************************************************* 

1080 GOSUB Clear_keys 

1090 GOSUB Init_param 

1 100 GOSUB Init_periph 

1110 GOSUB Input_mn_data 

1120 GOSUB Create_new_file 

1130! GOSUB Draw_freq_time 

1140 GOSUB Draw_mag_time 

1150 GOSUB Take_plot_data 

1 1 60 GOSUB Close_file 

1170 GOSUB Free j)eriph 

1180 GOSUB Postj>rocess 

1190 STOP 



1210 Existingjdata: ! Main program driver for process existing data 

1230 GOSUB Clear_keys 

1240 GOSUB Open_old_fUe 

1250 GOSUB Post_process 

1260 STOP 



12701 ******************************************************************* 



1280 Editing_data: ! Main program driver to compute Redfreq for old data 
1 !290 1***^****************************************^^^*^*^^*^^*^*^^*^****^ 
1300 Flag=l 

1310 GOSUB Clear_keys 

1320 GOSUB Open_old_file 

1330 GOSUB Calc_redfreq 

1340 GOSUB Create_new_file 

1350 GOSUB Close_file 

1360 STOP 

1370!******************************************************************* 
1380 Clear_keys: ! "Turns off all softkeys and clears the screen 

1390!******************************************************************* 



1400 FOR Keynumber=0 TO 9 

1410 ON KEY Keynumber LABEL "" GOSUB Commentl 

1420 NEXT Keynumber 
1430 OUTPUT KBD;Clear$; 

1440 RETURN 



1 ^ ^ 0 1 sick 9k 3k 3k % sk ^ ^ :(( 9)n|( :|o(< 9(e :(( 3|e ^ 



1460 Init_param: ! User sets the minimum and maximum values for initial plot 

2 ^70 1 sk 3k 9k 3k 3k 9k 3k 3k 9k 3k 3k 3k 3k 3k 3k 9k 3k 3k 3k 3k 9k 3k 9k 3k 3k 3k 9k 9k 3k 9k 3k 3k 9k 3k 3k 3k 9k 3k 3k 3k 3k 9k 9k 3k 3k 3k 9k 9k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 

1480 Center=(Screen-46)/2 ! Leading spaces for centering 

1490 PRINT TABXY(1 ,1) ! Start at top with blank line 

1500 PRINT TAB(Center);"Key Purpose" 

1510 PRINT TAB(Center);" " 

1520 PRINT TAB(Center);" 0 Set minimum resonance frequency" 

1530 PRINT TAB(Center);" 1 Set maximum resonance frequency" 
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1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610! 

1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720! 

1730Blinkl 

1740 

1750 

1760 

1770 



PRINT TAB(Center);’’ 2 
PRINT TAB(Center);" 3 
PRINT TAB(Center);" 5 
PRINT TAB(Center);" 6 
PRINTTAB(Center);"7 
PRINT TAB(Center);" 8 
PRINT TAB(Center);" 9 



Set minimum in phase magnitude" 

Set maximum in phase magnitude" 

Set minimum temperature in degrees, C" 
Set maximum temperature in degress, C" 
Set the default values" 

Review the assigned values" 

To proceed" 



ON KEY 0 LABEL "Min Freq, Hz" GOSUB Min_frequency 
ON KEY 1 LABEL "Max Freq, Hz" GOSL!B Max_frequency 
ON KEY 2 LABEL "Min Mag, DCV" GOSUB Min_mag 
ON KEY 3 LABEL "Max Mag, DCV" GOSUB Max_mag 
ON KEY 4 LABEL "" GOSUB Commentl 
ON KEY 5 LABEL "Min Temp, C" GOSUB Min_temperature 
ON KEY 6 LABEL "Max Temp, C" GOSUB Max_temperature 
ON KEY 7 LABEL "Default Values" GOSUB Default_values 
ON KEY 8 LABEL "Review values" GOSUB Review_values 
ON KEY 9 LABEL "Proceed" GOTO Moveon 



! Wait for user selection and 

! then take appropriate action 



WATT 1 

DISP "Select an option" 

WATT 1 
DISP 

GOTO Blink 1 
1780 Moveon: GOSUB Clear_keys 

1790 DISP "Proceeding " 

1800 WATT 1 

1810 DISP 

1820 RETURN 

1 8 3 0 i*** * 



1840 Commentl: ! Alerts user when an unassigned soft key is selected 

1850 ! ************************************************** ****’•'***’•'*******’•' 

1860 BEEP 300, .1 

1870 DISP "This soft key is unassigned" 

1880 WAIT 1 

1890 DISP 

1900 RETURN 

1910! ******************************************************************* 

1920 Default_values: ! Assigns default values for minimum and maximum T, F 

1930!******************************************************************* 

1940 User_freq_min=1650 

1950 User_freq_max=1900 

1960 User_temp_min=0 

1970 User_temp_max=25 

1 980 DISP "The Default Values are Set" 

1990 WAIT 1 

2000 GOSUB Review_values 

2010 RETURN 

2020 ! ** ******** ^ ******* *^** ************************* **^*^^*^*^**^^ ****** 
2030 Min_frequency: ! Accepts user input for minimum frequency 
2040!*****************************************************************’*'* 
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2050 INPUT "The minimum frequency for the pl< in Hz ?",User_freq_min 
2060 DISP "Ti>e minimum frequency is set at: ";Ij ,er_freq_min;" Hz" 

2070 WAIT 1 

2080 DISP 

2090 RETURN 






2110 Max_frequency: ! Accepts user input for maximum frequency 

2 20 |:|C3k3k3|<9|<3iC3i<3i(:|C3i(3|(3i(9|(3((3|<:|C3i(:k9k4C9kN(9k9i(9k)k3i<:|(3i(9k9k9k9k9|(9|(3k3k^3i<^4C9k3i(9k3(c3i(9k3k9k9k9i(3|C3k3|<3(C9|<3k^3|<3(C9i(3k9k3i<9k9i<9k 



2130 INPUT "The maximum frequency for the plot, in Hz ?",User_freq_max 
2140 DISP "The maximum frequency is set at: ";User_freq_max;" Hz" 

2150 WATT 1 

2160 DISP 

2170 RETURN 



2190 Min_mag: ! Accepts user input for minimum magnitude 
2200 !******************************************************************* 
2210 INPUT "The minimum magnitude for the plot, in DCV ?",User_mag_min 
2220 DISP "The minimum magnitude is set at: ";User mag min 

2230 WATT 1 

2240 DISP 

2250 RETURN 

2260! ******************************************************************* 



2270 Max_mag: ! Accepts user input for maximum magnitude 
2280!******************************************************************* 
2290 INPUT "The maximum magnitude for the plot, in DCV ?",User_mag_max 

2300 DISP "The maximum magnitude is set at: ":User_mag max 

2310 WATT 1 

2320 DISP 

2330 RETURN 



2340!******************************************************************* 



2350 Review_values: ! Presents currently assigned values for review 
2360!******************************************************************* 



2370 

2380 

2390 

2400 

2410 

2420 

2430 

2440 

2450 

2460 

2470 

2480 

2490 

2500 



DISP "Minimum frequency; ";User_freq_min;" Hz" 
WATTl 

DISP " Maximum frequency: ";User_freq_max;" Hz" 
WATTl 

DISP "Minimum magnitude: ";User_mag_min;" DCV" 
WATTl 

DISP " Maximum magnitude: ";User_mag_max;" DCV" 
WATTl 

DISP "Minimum temperature: ";User_temp_min;" C" 
WATTl 

DISP " Maximum temperatme: ";User_temp_max;" C" 

WATT 1 

DISP 

RETURN 



2510!******************************************************************* 



2520 Min_temperature; ! Accepts user input of minimum temperature 
2530!******************************************************************* 



2540 INPUT "Minimum temperature, degrees C ?",User_temp_min 
2550 DISP "The minimum temperature is set at: ";User_temp_min;" C" 
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2560 WATT 1 
2570 DISP 

2580 RETURN 






2600 Maxjemperature; ! Accepts user input of maximum temperature 
2610 !******************************************************************* 



2620 INPUT "Maximum temperature, degrees C ?",User_temp_max 
2630 DISP "The maximum temperature is set at: ";User_temp_max;" C" 

2640 WATT 1 

2650 DISP 

2660 RETURN 



2670 !******************************************************************* 



2680 Init_periph: ! Initializes the voltmeter and the frequency counter 

2690 1 *^^^^^^^^^************************************************^*******^ 

2700! 720 - HP5316A Universal Counter 

2710! 722 - HP3456A Digital Voltmeter 

2720! 724 - HP3478A Digital Multimeter (added on 23 Jan 91) 

2730! 

2740 OUTPUT 720; "IN" ! Default state 

2750 OUTPUT 722;"HF4R1M66STG100STI" ! 2 wire ohms, THMS degrees C 

2760 ! 100 line cycles integration 

2770 OUTPUT 724;"F1RAN5T1Z1D1" 

2780 RETURN 

2790!******************************************************************* 



2800 Free_periph: ! Frees the voltmeter and the counter from the HPIB bus 

2^ JQ|a|(9|(3fc}|(9|c3|(]f(3|(9|<3|<3|<3{c9|(3{(3|<3|c3|<3((9|<3{C9|(3{c3|(9|(3|(9)c9|(3|(9|(:f(3|(3(C]f(]f(:(C3|(9|<3|t3|(3|(3|<]f(9|<]f(]f(9((;f(3|(:)c%9(c;f(9|(9|(3fc]f(jf(}|c9|(:)(:)C3|(3(c}|t3|(9|<:tc 

2820 LOCAL 720 

2830 LOCAL 722 

2840 LOCAL 724 

2850 RETURN 

2^601******************************************************************* 



2870 Input_run_data; ! Accepts user input of selected run data 



2880! 

2890 

2900 

2910 

2920 

2930 

2940 

2950 

2960 

2970 

2980 

2990 

3000 

3010 

3020! 

3030 

3040 

3050 

3060 



aft s|c + + a(< + a|«3|ta(« a|« 3|t aft a|t + 3#« 3*e ;#t + :♦« a(< s|c + ;#c 3*« 3|c ;*t 3*« ;4ea(«J|« + 3*« + 3»c 3f« % s(c % 3|c 



Center=(S creen -42)/2 
PRINT TABXY( 1,1) 
PRINT TAB(Center);"Key 
PRINT TAB(Center);"- — 
PRINT TAB(Center);"0 
PRINT TAB(Center);" 1 
PRINT TAB(Center);" 2 
PRINT TAB(Center);" 3 
PRINTTAB(Center);"5 
PRINTTAB(Center);"6 
PRINT TAB(Center);" 7 
PRINT TAB(Center);" 8 
PRINT TAB(Center);" 9 



! Leading spaces for centering 
! Start at top with blank line 
Purpose" 



Enter rod identification" 

Enter the run number for this mode 
Enter the mode for this mn" 

Enter the date for this run" 

Enter the mass for this rod" 

Enter the length for this rod" 

Enter the diameter for this rod" 
Review the information entered" 

To proceed" 



ti 



ON KEY 0 LABEL "Rod ID" GOSUB Rod_id 
ON KEY 1 LABEL "Run No." GOSUB Run_number 
ON KEY 2 LABEL "Mode" GOSUB Mode 
ON KEY 3 LABEL "Date" GOSUB Date 
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3070 ON KEY 4 LABEL "" GOSUB Commentl 
3080 ON KEY 5 LABEL "Rod Mass" GOSUB Rod_mass 
3090 ON KEY 6 LABEL "Rod Ungth" GOSUB Rodjength 
3 100 ON KEY 7 LABEL "Rod Diameter" GOSUB Rod_diameter 

3110 ON KEY 8 LABEL "Review entries" GOSUB Review_entrys 

3120 ON KEY 9 LABEL "Proceed" GOTO Onward 

3130 Blink2; WATT 1 

3140 DISP "Select an option or proceed" 

3150 WATT 1 

3160 DISP 

3170 GOTOBlink2 

3180 Onward: Volume=(Length/100)*Pi*.25*(Diameter/100)^2 
3190 Density=(Mass/l000)/Volume 

3200 GOSUB Clear_keys 

3210 DISP "Proceeding " 

3220 WATT 1 

3230 DISP 

3240 RETURN 

^ 2 ^ Q I ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

3260 Rod_id: ! Accepts user input of rod identification 
^ 0 ! ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 

3280 INPUT "Rod identification block label (i.e. ECP4) ?",Blockl$ 

3290 WHILE LEN(Blockl$)>12 

3300 DISP "Limit Rod Identification to 12 characters" 

3310 BEEP 300,. 1 

3320 WATT 1 

3330 DISP 

3340 INPUT "Rod identification block label (i.e. ECP4) ?"31ockl$ 

3350 END WHILE 

3360 DISP "The Rod ID is set at: ";Blockl$ 

3370 WATT 1 

3380 DISP 

3390 RETURN 

3410 Run_number: ! Accepts user input of mn number 

3420 ! ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ ^ 
3430 INPUT "Run number for this mode (i.e. 4) ?",Block2$ 

3440 DISP "The Run number is set at: ";Block2$ 

3450 WATT I 

3460 DISP 

3470 RETURN 

3490 Mode: ! Accepts user input of the mode 

3500 !******************************************************************* 
3510 INPUT "Mode (i.e. Torsional, Longitudinal, Flexural) ?",Block3$ 

3520 SELECT UPC$(Block3$) 

3530 CASE "FLEXURAL","TORSIONAL", "LONGITUDINAL" 

3540 DISP "The mode is set as: ";Block3$ 

3550 WATT 1 

3560 DISP 

3570 CASE ELSE 
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3580 BEEP 300,. 1 

3590 DISP "Choices are: Torsional, Longimdinal or Flexural" 

3600 WATT 1 

3610 DISP 

3620 GOTO Mode 

3630 END SELECT 

3640 RETURN 

3650 !************************* *********’*** ********************* ********** 



3660 Date: ! Accepts user input of run date 



3680 INPUT "The date for this sample (i.e. 10 SEP 89) ?",Block4$ 

3690 WHILE LEN(Block4$)>9 

3700 DISP "Limit date entry to 9 characters" 

3710 BEEP300,.! 

3720 WATT 1 

3730 DISP 

3740 INPUT "The date for this sample (i.e. 10 SEP 89) ?",Block4$ 

3750 END WHILE 

3760 DISP "The Date is set at: ";Block4$ 

3770 WATT 1 

3780 DISP 

3790 RETURN 



3800! ******************************************************************* 



3810 Rod_mass: ! Accepts user input of the rod mass 

3820!******************************************************************* 
3830 INPUT "The mass for this rod (units: grams) ?",Mass 

3840 DISP "The mass is set at: ";Mass;" grams" 

3850 WATT 1 

3860 DISP 

3870 RETURN 



3880!******************************************************************* 



3890 Rod_length: ! Accepts user input of the rod length 

3910 INPUT "The length for this rod (units: centimeters) ?",Length 
3920 DISP "The length is set at: ";Length;"centimeters" 

3930 WATT 1 

3940 DISP 

3950 RETURN 

3960!******************************************************************* 
3970 Rod_diameter: ! Accepts user input of the rod diameter 

3980!******************************************************************* 
3990 INPUT "The diameter for this rod (units: centimeters) ?", Diameter 
4000 DISP "The diameter is set at: ";Diameter;" centimeters" 

4010 WATT 1 

4020 DISP 

4030 RETURN 



4050 Review_entrys: ! Presents currently assigned values for review 

4070 DISP "Rod ID: ";Blockl$ 

4080 WATT 1 
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4090 


DISP " Run: ";Block2$ 


4100 


WATT 1 


4110 


DISP "Mode: ";Block3$ 


4120 


WAJTl 


4130 


DISP " Date: ";Block4$ 


4140 


WATT 1 


4150 


DISP "Mass: ";Mass; "grams" 


4160 


WATTI 


4170 


DISP " Length: ";Length;" centimeters" 


4180 


WATTI 


4190 


DISP "Diameter: ";Diameter;" centimeters' 


4200 


WATT I 


4210 


DISP 


4220 


RETURN 



4230!******************************************************************* 



4240 Draw_mag _time: I Produces Magnitude vs Time graph w/o curve 



^25Q !******************************************************** *********** 



4260 Main_tide$=Label$(14) 

4270 Sub_title$=Label$(2) 

4280 X_axis_name$=Label$(3) 

4290 Y_axis_name$=Label$(13) 

4300 Xmin=0 

4310 Xmax=20 

4320 Y min=Actual_mag_min 

4330 Ymax=Actual_mag_max 

4340 GOSUB Generic_plot 

4350 RETURN 

4360!******************************************************************* 



4370 Draw_freq_time: ! Produces Frequency vs Time graph w/o curve 

4390 Main_title$=Label$(l) 

4400 Sub_title$=Label$(2) 

4410 X_axis_name$=Label$(3) 

4420 Y_axis_name$=Label$(4) 

4430 Xmin=0 

4440 Xmax=20 

4450 Ymin=Actual_freq_min ! Actual minimum frequency 

4460 Ymax=Actual_freq_max ! Actual maximum frequency 

4470 GOSUB Generic_plot 

4480 RETURN 

4490!******************************************************************* 



4500 Young mod: ! Compute Young's modulus in flexure modes 
4510 !******************************************************************* 



4520 D=Diameter/100 

4530 L=Length/100 

4540 Arg6=Density*((32*L^2)/(Pi*D*Fmode^2))^2 

4550 Eflex(I)=Arg6*Resfreq(I)''2 

4560 RETURN 

4570!******************************************************************* 



4580 Shear_mod: ! Compute Shear Modulus 

QQ f 9k ak 3k 3k 3k 3k sk 3k 3k 3k 4t 3k ak 3k 3k 3k ak 3k ak 3k sk ak :>k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k sk 3k 3k 3k 3k 3k 3k sk 3k 3k ak 3k 3k ak 3k ak 3k ak 3k 3k 3k 3k 3k ak 3k 9k 3k ak 3k 
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4600 D=Diameter/100 

4610 L=Length/100 

4620 Arg5=(Density*4*L^2)/(Tmode^2) 

4630 Gtor(J)=(Arg5*Resfreq(J)^2) 

4640 RETURN 

4650 !********************************************************************* 

4660 Lyoung_mod: ! Compute Young's modulus in longitudinal mode 

46*70 !************************************************************ ******* 

4680 D=Diameter/100 

4690 L=Length/100 

4700 Arg7=(Density*4*L''2)/(Lmode^2) 

4710 Elong(K)=Arg7*Resfreq(K)''2 

4720 RETURN 



4'730J ********************************************** ********************* 



4740 Compute_redfreq: ! Compute redfreq of new data 

4750!******************************************************************* 
4760 Temp(I)=Thermistor(I)+273.15 

4770! X(I)=-12.9*(Temp(I)-283.15)/(107+Temp(I)-283.15) ! PR 1592 only 

4780 X(I)=-21.5*(Temp(I)-211)/(43.1+Temp(I)-211) ! Plexi-glass only 

4790 Y(I)=X(I)*LOG(10) 

4800 Alphat(I)=EXP(Y(I)) 

48 1 0 Redfreq(I)=Resfireq(I)* Alphat(I) 

4820 RETURN 



4840 Calc_redfreq: ! Convert RT20B to RT20C format 

4850! ******************************************************************* 

4860 FOR 1=0 TO 20* Arg 1 

4870 Temp(I)=Thermistor(I)+273. 1 5 

4880 X(I)=-12.9*(Temp(I)-283.15)/(107+Temp(I)-283.15) IPR 1592 only 

4890 Y(I)=X(I)*LOG(10) 

4900 Alphat(I)=EXP(Y(D) 

4910 Redfreq(I)=Resfreq(I)*Alphat(I) 

4920 NEXT I 

4930 RETURN 

4940!******************************************************************* 
4950 Compute_modulus: I Computes appropriate material modulus based on mode 
4960!******************************************************************* 

4970 DISP "Computing appropriate modulus " 

4980 D=Diameter/100 ! convert to meters 

4990 L=Length/100 ! convert to meters 

5000 SELECT UPC$(Block3$) 

5010 CASE "FLEXURAL" 

5020 Arg2=Density*((32*L''2)/(Pi*D*Fmode''2))'^2 

5030 Eflex_max=Arg2*Actual_freq_max^2 

5040 Powep=LOG(Eflex_max) DFV LOG(IO) 

5050 Scale_factor=10^Power 

5060 Scale_eflex_max=Eflex_max/Scale_factor 

5070 Scale_eflex_min=(Arg2*Actual_freq_min^2)/Scale_factor 

5080 FOR 1=0 TO 20* Arg 1 

5090 Eflex(I)=(Arg2*Resffeq(I)^2)/Scale_factor 

5100 NEXT I 
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5110 CASE "TORSIONAL" 

5120 Arg3=(Density*4*L^2)/(Tmode^2) 

5130 Gtor_max= Arg3 * Actual_freq_max^2 

5 1 40 Power^LOG(Gtor_inax) DIV LOG( 1 0) 

5150 Scale_factor=10^Power 

5160 Scale_gtor_max=Gtor_max/Scale_factor 

5170 Scale_gtor_min=( Arg3* Actual_freq_min'^2)/Scale_factor 

5180 FOR J=0 TO 20* Arg 1 

5190 Gtor(J)=(Arg3*Resfreq(J)^2)/Scale_factor 

5200 NEXT! 

5210 CASE "LONGITUDINAL" 

5220 Arg4=(Density*4*L'^2)/(Lmode''2) 

5230 Elong_max=Arg4*Actual_fTeq_max'^2 

5240 Power=LOG(Elong_max) DIV LOG(IO) 

5250 ScaIe_factor=10'^Power 

5260 ScaIe_elong_max=Elong_max/ScaIe_factor 

5270 ScaIe_elong_min=(Arg4*Actual_freq_min'^2)/ScaIe_factor 

5280 FOR K=0 TO 20*Argl 

5290 Elong(K)=(Arg4*Resfreq(K)''2)/Scale_factor 

5300 NEXT K 

5310 END SELECT 

5320 DISP 

5330 RETURN 



^ j Jfe ♦♦♦♦ a|<5|<* 5|<*3|<a|<5|<>K ale** aH* ak* ak* a|<*5|< 4<5|<5|< Jle* sfe ** ak ale***** ♦ 



5350 Percent_moduli: ! Determines the percent change in modulus per degree C 
^ 

5370 Delta_temp=ActuaI_temp_max-Actual_temp_min 

5380 Fsum2pt=Actual_freq_max+Actual_freq_min 

5390 Fdiff2pt=Actual_freq_max-Actual_freq_min 

5400 Fminls=Actual_temp_max*Slope+Intercept 

5410 Fmaxls=Actual_temp_min*Slope+Intercept 

5420 Fsumls=Fmaxls+Frninls 

5430 Fdiffls=Fmaxls-Fminls 

5440 Prcnt_per_c_2pt=100*4*(Fdiff2pt/Fsum2pt)/DeIta_temp ! 2 point 

5450 Prcnt_per_c_ls=100*4*(Fdiffls/Fsumls)/Delta_temp ! least squares 

5460 RETURN 



^ Q j 3k ak ak ak ak ak ak ak ak ak a(c ak ak ak ak ak ak ak ak ak ak a(c a(c ak a(c a(< a|< afe ak ak ak ak a(e ak ak ak ak a(c ak ak a|c a(c ak a(e ak a(e a(c ak afe ak ak a(e ak ak a|c ak ak ak afe a{e a|c afe ak a(c ak afe ak 



5480 Mag_temp: ! Produces In phase mag vs Temperature graph w/ curve 

1 3(( 3k 3k ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak ak 

5500 Main_title$=Label$(15) 

5510 Sub_title$=Label$(2) 

5520 X_axis_name$=Label$(5) 

5530 Y_axis_name$=Label$(13) 

5540 Xmin=Actual_temp_min 

5550 Xmax=Actual_temp_max 

5560 Y min= Actu al_mag_min 

5570 Ymax=Actual_mag_max 

5580 GOSUB Generic_plot 

5590 CALL Generic_curve(Pen3,Thermistor(*),Inphmag(*),20*Argl) 

5600 RETURN 
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5 ^ 1 0 !***** ’^^ ******************************************** ’^^ ********* ’^^ ***** * 
5620 Mag_freq: IProduces In phase mag vs Resonant frequencies 

5630!******************************************************************* 
5640 Main_title$=Label$(16) 

5650 Sub_title$=" " 

5660 X_axis_name$=Label$(4) 

5670 Y_axis_name$=Label$(13) 

5680 Xmin=ActuaI_freq_min 

5690 Xmax=Actual_freq_max 

5700 Ymin=Actual_mag_min 

5710 Y max= Actual_mag_max 

5720 GOSUB Generic_plot 

5730 CALL Generic_curve(Pen3,Resfreq(*),Inphmag(*),20*Argl) 

5740 RETURN 

5750!******************************************************************* 



5760 Shear_temp: ! Produces Shear modulus vs Temperature gr^h w/ curve 
5770!******************************************************************* 



5780 

5790 

5800 

5810 

5820 

5830 

5840 

5850 

5860 

5870 

5880 



Main_title$=Label$( 10) 

Sub_title$=" " 

X_axis_name$=Label$(5) 

Y_axis_name$=Label$(l l)&VAL$(Power) 
Xmin=Actual_temp_min ! Actual minimum temperature 

Xmax=Actual_temp_max 1 Actual maximum temperature 

Ymin=Scale_gtor_min I Scaled minimum shear modulus 

Ymax=Scale_gtor_max I Scaled maximum shear modulus 

GOSUB Generic_plot 

CALL Generic_curve(Pen3 ,Thermistor(* ),Gtor(* ),20* Arg 1 ) 
RETURN 



5890 ! *^^*^* ********** *^^* ******************** *^** ************* ^*^**^^^^^ 



5900 Young_temp: ! Produces Young's modulus vs Temperature graph w/ curve 

5910!******************************************************************* 



5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 

6070 

6080 

6090 

6100 



Main_title$=Label$(9) 

Sub_title$=" " 

X_axis_name$=Label$(5) 

Y_ax is_name$=Label$( 1 2)& V AL$(Po wer) 



Xm in= Actual_temp_min 
Xmax=Actual_temp_max 
SELECT UPC$(Block3$) 

CASE "FLEXURAL" 

Y min=Scale_eflex_min 

Y max=Scale_eflex_max 
GOSUB Generic_plot 
CALL Generic_curve(Pen3,Themiistor(*),Eflex(*),20*Argl ) 

CASE "LONGITUDINAL" 



! Actual minimum temperature 
! Aaual maximum temperature 



I Scaled minimum Y (flexural) 

! Scaled maximum Y (flexural) 



Ymin=Scale_elong_min ! Scaled minimum Y (longitudinal) 
Ymax=Scale_elong_max ! Scaled maximum Y (longitudinal) 
GOSUB Generic_plot 

C A LL Generic_curve (Pen3 ,Thermistor(* ) ,Elong( * ) ,20* Arg 1 ) 
END SELECT 
RETURN 
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6120 Generic_plot: ! Produces design and layout of line graph 



6130! 

6140 

6150 

6160 

6170 

6180! 

6190 

6200 

6210 

6220 

6230 

6240 

6250 

6260 

6270! 

6280 

6290 

6300 

6310 

6320 

6330! 

6340 

6350 

6360 

6370 

6380 

6390! 

6400 

6410 

6420 

6430 

6440 

6450 

6460 

6470! 

6480 

6490 

6500 

6510 

6520 

6530 

6540 

6550 

6560 

6570 

6580 

6590 

6600 

6610 



GINTT ! Initialize various graphics parameters 

GCLEAR ! Clear the graphics display 

GRAPHICS ON ! Turn the graphics display on 

OUTPUT KBD;aear$; ! Clear the CRT 

IF Plot_device$="Plotter" THEN ! Route output to HP7475A plotter 
GRAPHICS OFF 
BEEP 1000,. 1 

DISP "Press CONTINUE when the plotter is ready" 

PAUSE 

OUTPUT 705;"IP 1300,1000,9000,6750;" ! Sets scaling points 

PLOTTER IS 705,"HPGL" 

END IF 

Xgdumax=100*MALX(l, RATIO) ! How many gdu's wide screen is 

Ygdumax=100*MAX( 1,1 /RATIO) ! How many gdu's high screen is 

CALL Scale(.02,.98,.10,.98,Top,Left,Xcntr,Ycntr,Xgdumax,Ygdumax) 
PEN Pen2 
FRAME 

IF LEN(Main_title$)>30 THEN ! Size main tide 

CALL Label(4,.6,0,6,Pen2,Xcntr,Top-2AIain_title$) 

ELSE 

CALL Label(6,.6,0,6JPen2,Xcntr,Top-2AIain_title$) 

END IF 



CALL Label(4,.6,0,6,Pen2,Xcntr,.9*Top,Sub_tide$) 

CALL Block_info(B lock 1 $ ,B lock2$ ,Block3 $ ,Block4$ ,Xgdumax ,Pen5 ,Pen2) 
CALL Scale(.15,.90,.25,.8,Top,Left,Xcntr,Ycntr,Xgdumax,Ygdumax) 

CLIP OFF 

CALL Label(4,.6,0,4,Pen5,Xcntr,-l 8,X_axis_name$) 

CALL Label(4,.6,90,4,Pen5,-l 2,Ycntr,Y_axis_name$) 

CUP ON 



SELECT Xscale$ 

CASE "Auto scale X" ! Auto scale X axis 

CALL Xscale(Xmin,Xmax,Xrninor,Xmajor) 



CASE "User scale X" 
SELECT Pass$ 

CASE "Initial" 

Xmin=0 
Xmax=20 
Xmajor=l 
Xminor=.5 
CASE "FoUow on" 
Xmin=Xmin_manual 
Xmax=Xmax_manual 
Xmajop=Xmajor_manual 



! User scale X axis 

! Initial pass 
! Minimum time, hours 
! Maximum time, hours 



! Use user modified values 
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6620 Xminor=Xminor_manual 

6630 END SELECT 

6640 END SELECT 
6650! 

6660 SELECT Yscale$ 

6670 CASE "Auto scale Y" ! Auto scale Y axis 

6680 CALL Yscale(Ymin,Ymax,Yminor,Ymajor) 

6690 CASE "User scale Y" ! User scale Y axis 

6700 SELECT Pass$ 

67 10 CASE "Initial" ! Initial pass 

6720! Yinin=User_freq_min ! Use user entered min fireq 

6730! Ymax=User_freq_max ! Use user entered maximum freq 

6740 Ymin=User_mag_min ! Use user entered min mag 

6750 Yiuax=User_mag_max ! Use user entered max mag 

6760 Ymajor=(Yinax-Ymin)/5 

6770 Yminor=Ymajor/5 

6780 Pass$="Follow on" 

6790 Xscale$="Auto scale X" 

6800 Yscale$="Auto scale Y" 

68 1 0 CASE "Follow on" ! Use user modified values 

6820 Ymin=Yinin_manual 

6830 Ymax=Ymax_manual 

6840 Y major= Y major_m;uiual 

6850 Yminor=Yminor_manual 

6860 END SELECT 

6870 END SELECT 

6880 WINDOW Xmin.Xmax.Ymin.Ymax 

6890! 

6900 CALL Lbl_axes(2..6,Pen4,Xmin,Xmax,Xmajor,Yinin,Ymax,Ymajor) ! Axes 

6910 PEN Pen2 

6920 AXES Xminor,Yminor,Xmin,Ymin,Xmajor/Xmiiior,YmajorA^ minor, 3 

6930 AXES Xminor,Yminor.Xmax,Ymax,Xmajor/Xminor,YmajorA" minor, 3 

6940! 

6950 SELECT Grid_type$ ! Grid 

6960 CASE "Partial" 

6970 GRID Xmajor,Ymajor,Xmin,Ymin 

6980 CASE "Full" 

6990 GRID Xminor,Yminor,Xmin,Ymin,Xmajor/Xminor,YmajorA^ minor,l 

7000 END SELECT 

7010! 

7020 PENUP 

7030 GRAPHICS ON 

7040 RETURN 

7050 ! ********************** *************’''*******’''****’*'** **************’*'* 
7060 Take_plot_data: ! Collects frequency vs temperature data during the mn 

7080 ON KEY 0 LABEL "Take Data" GOTO Take 

7090 FOR Keyiuimben=l TO 9 

7100 ON KEY Keynuml->er LABEL "" GOSUB Comment 1 

7110 NEXT Keynumber 

7120 Wait; GOTO Wait 
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7130 Take: ON KEY 0 LABEL "" GOSUB Commentl 
7140 PEN Pen3 

7150 OUTPUT KBD;"Taking data 

7160 T0=TIMEDATE 
7170 SELECT UPC$(Block3$) 

7180 CASE "FLEXURAL" 

7 1 90 FOR 1=0 TO 20* Arg 1 

7200 ENTER 720;Resfreq(I) 

7210 ENTER 722;Thermistor(I) 

7220 ENTER 724;Inphmag(I) 

7230 Time(I)=I/Argl 

7240 PLOT Time(I) Jnphmag(I) 

7250 GOSUB Compute_redfreq 

7260 GOSUB Young_mod 

7270 OUTPUT @Path_l;Eflex(I),Thermistor(I),Resfreq(I)Jnphmag 

(I),Redfreq(I) 

7280 DISP I 

7290 WATT 295.787 lAdjustment for 1 sample every 5 minutes 

7300 NEXT I 

7310 CASE "TORSIONAL" 

7320 FOR J=0 TO 20*Argl 

7330 ENTER 720;Resfreq(J) 

7340 ENTER 722;Thermistor(J) 

7350 ENTER 724;Inphmag(J) 

7360 Time(J)=J/Argl 

7370 PLOT Time(J),Inphmag(J) 

7380 GOSUB Compute_redfreq 

7390 GOSUB Shear_mod 

7400 OUTPUT (2)Path_l;Gtor(J),Thermistor(J),Resffeq(J),Inphmag(J 

),Redfreq(J) 

7410 DISPJ 

7420 WATT 295.787 

7430 NEXT! 

7440 CASE "LONGITUDINAL" 

7450 FOR K=0 TO 20*Argl 

7460 ENTER 720;Resfireq(K) 

7470 ENTER 722;Thermistor(K) 

7480 ENTER 724;Inphmag(K) 

7490 Time(K)=K/Argl 

7500 PLOT Time(K)Jnphmag(K) 

7510 GOSUB Compute_redfreq 

7520 GOSUB Lyoung_mod 

7530 OUTPUT @Path_l;Elong(K),Thermistor(K),Resfreq(K),Inphmag( 

K),Redfreq(K) 

7540 DISP K 

7550 WATT 295.787 

7560 NEXTK 

7570 END SELECT 

7580 T1=TIMEDATE 

7590 DISP "IT TOOK ";DROUND(Tl -TO,4);"SECONDS" 

7600 WATT 10 
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! Clears the CRT 



7610 OUTPUT KBD;Clear$; 

7620 BEEP 
7630 ON KEY 0 LABEL "Post Process" GOTO 7650 
7640Waitl; GOTOWaitl 
7650 RETURN 



7660!******************************************************************* 
7670 Post_process: ! Permits user to extract plots, tables, and other info 
7680!******************************************************************* 



7690 

7700 

7710 

7720 

7730 

7740 

7750 

7760 

7770 

7780! 

7790 

7800 

7810 

7820 

7830 

7840 

7850! 

7860 

7870 

7880 

7890! 

7900 

7910 

7920 

7930 

7940 

7950 

7960 

7970 

7980 

7990 

8000 

8010 

8020 

8030 

8040 

8050 

8060 

8070 

8080 

8090 

8100 

8110 



GCLEAR 

CALL Least_squares(20*Argl ,Thermistor(*)JResfreq(*) A>B.C,D,E,F,G) 

Slope=A ! Keeps call to Least_squares to one program line 

Intercept=B 

Correlation=C 

Slope_error=D 

Intrcpterr=E 

Tmean=F 

Fmean=G 



Actual_temp_min=MIN(Thermistor(*)) ! Find maximum and minimum 
Actual_temp_max=MAX(Thermistor(*)) ! temperature and frequency 
Actual_freq_min=MIN(Resfreq(* )) 

Actual_freq_max=MAX(Resfreq(*)) 

Actual_mag_min=MIN(Inphmag(*)) 

Actu al_mag_max=MAX(Inphmag(* )) 



Volume=.25*(Length/100)*Pi*(Diameter/100)^2 
GOSUB Percent_moduli 
GOSUB Compute_modulus 



OUTPUT KBD;Clear$; ! Clear the CRT 

OUTPUT KBD;Home$; ! Home display 



GCLEAR 

PRINT TAB XY(1,1) ! Start at top with bleank line 

Center=(Screen-42)/Z ! Leading spaces for centering 

PRINT TAB(Center);"Key Purpose" 

PRINT TAB(Center);" " 

PRINT TAB(Center);" 0 Plot Frequency vs Time" 

PRINT TAB(Center);" 1 Plot Temperature vs Time" 

PRINT TAB(Center);" 2 Plot Frequency vs Temperature" 
SELECT UPC$(Block3$) 

CASE ="TORSIONAL" 

PRINT TAB(Center);" 3 Plot Shear modulus vs Temperature" 
CASE ="FLEXURAL" 



PRINT TAB(Center);" 3 Plot Young's modulus vs Temperature" 
CASE ="LONGITUDINAL" 

PRINT TAB(Center);" 3 Plot Young's modulus vs Temperature" 
END SELECT 



PRINT TAB (Center);" 4 Send Information to the Printer" 
PRINT TAB(Center);" 5 Select the graph output device" 
PRINT TAB(Center);" 6 Plot In phase mag vs Temperature" 
PRINT TAB(Center);" 7 Plot In phase mag vs Temperature" 
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8120 PRINT TAB(Center);" 7 Select the type of X axis scaling" 

8130 PRINT TAB(Center);" 8 Select the t^ of Y axis scaling" 

8 140 PRINT TAB(Center); " 9 Exit this program" 

8150! 

8160 ON KEY 0 LABEL "Freq vs Time" GOSUB Freq_time 
8170 ON KEY 1 LABEL "Temp vs Time" GOSUB Temp_time 

8180 ON KEY 2 LABEL "Freq vs Temp" GOSUB Freq_temp 

8 1 90 SELECT UPC$(Block3$) 

8200 CASE="TORSIONAL" 

8210 ON KEY 3 LABEL "G vs Temp" GOSUB Shear_temp 

8220 CASE ="FLEXURAL" 

8230 ON KEY 3 LABEL "E vs Temp" GOSUB Young_temp 

8240 CASE ="LONGITUDINAL" 

8250 ON KEY 3 LABEL "E vs Temp" GOSUB Young_temp 

8260 CASE ELSE 

8270 ON KEY 3 LABEL "" GOSUB Commentl 

8280 END SELECT 

8290 ON KEY 4 LABEL "Print Info" GOSUB Dump_info 

8300 ON KEY 5 LABEL "Output Device" GOSUB Output_device 

8310! ON KEY 6 LABEL "Grid Option" GOSUB Grid_option 

8320 ON KEY 6 LABEL "Imag vs Temp" GOSUB Mag_temp 

8330! ON KEY 7 LABEL "X scale option" GOSUB Xscale_option 

8340 ON KEY 7 LABEL "Imag vs Freq" GOSUB Mag_freq 

8350 ON KEY 8 LABEL "Y scale option" GOSUB Yscale_option 

8360 ON KEY 9 LABEL "Exit Program" GOSUB Program_end 

8370! 

8380 Blink3: WATT 1 

8390 DISP "Make a DECISION" 

8400 WATT 1 

8410 DISP 

8420 GOTO Blink3 

8430 RETURN 

*^**************************************************^*****^****^*** 
8450 Grid_option: ! Accepts the user's choice for the plot grid 

^ ^ j ^ 3k sk 3k 3k 3f< 3k 3k :f< :f( 3(< 3|c 3f( 3k 3k 3k 3k ak 3k 9k 3k 3k ak 3k ak 3k 9k 3k 3k ak 3k ak SK 3|< :|c ak 3k 3f< 3k 3k 3k 3k 3k 3k ak a(c 3k ak 3k 3k 

8470 GRAPHICS OFF ! Turns off graphics display 

8480 OUTPUT KBD;Clear$; ! Clears the CRT 

8490 DISP Grid_type$;" grid is currently selected" 

8500 INPUT "Enter F - full; P - partial; N - No grid ?",Response$ 

8510 SELECT Responses 

8520 CASE "F" 

8530 Grid_type$="FuU" 

8540 CASE "P" 

8550 Grid_type$="Partial" 

8560 CASE "N" 

8570 Grid_type$="No" 

8580 CASE ELSE 

8590 DISP "No change" 

8600 WATT 1 

8610 DISP 

8620 WATT 1 
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8630 END SELECT 

8640 DISP Grid_type$;" grid is selected" 

8650 WAIT 1 

8660 DISP 

8670 RETURN 






8690 Dump_info: ! Sends selected data and table information to the printer 
8700 !******************************************************************* 
8710 PRINTER IS 701 

8720 Perfskip$=CHR$(27)&CHR$(38)&CHR$(108)&CHR$(49)&CHR$(76) 

8730 Formfeed$=CHR$(12) 

8740 PRINT PerfskipS ! Skip on Perforation 

8750 PRINT Formfeeds 

8760 PRINT USING "3/,#" ! Three line feeds 

8770 PRINT "Rod: "&Blockl$;TAB(70);"Page 1 of 2" 

8780 PRINT "Run: "&Block2$ 

8790 PRINT "Mode: "&Block3$ 

8800 PRINT "Date: "&Block4$ 

8810 PRINT USING "3/,#" ! Three more line feeds 

8820! 

8830 Cl$="Time" 

8840 C2$="Temperature" 

8850 C3$="Frequency" 

8860 C4a$="Shear modulus" 

8870 C4b$="Young’s modulus" 

8880 C5$="In Phase Component" 

8890 C6$="hours" 

8900 C7$="deg C" 

8910 C8$="Hz" 

8920 C9$="Pa*10^" 

8930 C10$="DCV" 

8940! 

8950 SELECT UPC$(Block3$) 

8960 CASE "TORSIONAL" 

8970 PRINT USING 8980;Cl$,C2$,C3$,C4a$,C5$ 

8980 IMAGE 4A,6X.11A,5X,9A.5X,13A,5X,18A 

8990 PRINT 

9000 PRINT USING 9010;C6$,C7$,C8$,C9$,Power,C10$ 

9010 IMAGE 5A,8X,5A,11X,2A,12X,6A;ZZ,14X,3A,12X 

9020 PRINT 

9030 FOR 1=0 TO 20*Argl STEP Argl 

9040 PRINT USING 9050;Time(I),Thermistor(I),Resfreq(I),Gtor(I),In 

phmag(I) 

9050 IMAGE X,3D,8X,SDD.DDD,4X, 

DDDD.DDDE,9X,D.DDD,13X,SDD.DDDDDE 

9060 NEXT I 

9070 CASE "FLEXURAL" 

9080 PRINT USING 9140;Cl$,C2$,C3$,C4b$,C5$ 

9090 PRINT 

9100 PRINT USING 9010;C6$,C7$,C8$,C9$, Power, C10$ 

9110 PRINT 
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9120 FOP T=0 TO 20*Argl STEP Argl 

9130 P/ NT USING 9050;Time(J),Thermistor(J),Resfreq(J),Eflex(J) 

,Inphmag(I) 

9140 IMAGE 4A,6X,1 1A,5X,9A,5X,15A,3X,18A 

9150 NEXT! 

9160 CASE "LONGITUDINAL" 

9170 PRINT USING 9140;Cl$,C2$,C3$,C4b$,C5$ 

9180 PRINT 

9190 PRINT USING 9010;C6$,C7$,C8$,C9$, Power, C10$ 

9200 PRINT 

9210 FOR K=0 TO 20* Argl STEP Argl 

9220 PRINT USING 9050;Time(K),Thermistor(K),Resfreq(K),Elong(K) 

,Inphmag(I) 

9230 NEXT K 

9240 END SELECT 

9250! 

9260 PRINT Formfeeds ! Advance to top of next page 

9270 PRINT USING "3/,#" ! Three more line feeds 

9280 PRINT "Rod: "&Blockl$;TAB(70);"Page 2 of 2" 

9290 PRINT "Run: "&Block2$ 

9300 PRINT "Mode: "&Block3$ 

9310 PRINT "Date: "&Block4$ 

9320! 

9330 PRINT USING "3/,#" ! Three more line feeds 

9340 PRINT "Physical properties:" 

9350 PRINT 

9360 PRINT USING "3X,12A,14X,DDDDDD";"Mass, grams:"dVlass 

9370 PRINT USING "3X,20A,06X,DDD.DDD";"Length, centimeters: ".Length 

9380 PRINT USING "3X,22A,04X,DDD.DDD";"Diameter, centimeters: ".Diameter 

9390 PRINT USING "3X,21A,07X,D.DDDE";"Volume, cubic meters: ".Volume 

9400 PRINT USING "3X,16A,09X,DDDD.D";"Density, kg/m'^3:", Density 

9410! 

9420 PRINT USING "3/,#" ! Three more line feeds 

9430 PRINT "Least-squares fit results [frequency versus temperature]:" 

9440 PRINT 

9450 PRINT USING "3X,12A,13X,DDDD.DDD";"Slope, Hz/C:",Slope 

9460 PRINT USING "3X,18A,07X,DDDD.DDD";"Slope error, Hz/C:",Slope_error 

9470 PRINT USING "3X,14A,1 lX,DDDD.DDD";"Intercept, Hz: ".Intercept 

9480 PRINT USING "3X,20A,05X,DDDD.DDD";"Intercept error, Hz:",Intrcpterr 

9490 PRINT USING "3X,12A,13X,DDDDDDDDD";"Correlation:",Correlation 

9500 PRINT USING "3X,20A,05X,DDDD.DDD";"Mean temperature, C:",Tmean 

9510 PRINT USING "3X,19A,06X,DDDD.DDD";"Mean frequency, Hz:",Fmean 

9520 PRINT 

9530 PRINT 

9540 PRINT "Other statistics:" 

9550 PRINT 

9560 Tmin=Actual_temp_min 

9570 Tmax=Actual_temp_max 

9580 Tave=(Tmax+Tmin)/2 

9590 Fmin=Actuai_freq_min 

9600 Fmax=Actual_freq_max 
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9610 Fave=(Fmax+Fmin)/2 

9620 PRINT USING "3X,23A,02X,DDDD.DDD";"Minimum temperature, C:",Tmin 

9630 PRINT USING "3X,23A,02X,DDDD.DDD";’’Average temperature, C;",Tave 

9640 PRINT USING "3X,23A,02X,DDDD.DDD";"Maximum temperature, 

C:",Tmax 
9650 PRINT 

9660 PRINT USING "3X,22A,03X,DDDD.DDD";"Minimum frequency, Hz;",Fmin 

9670 PRINT USING "3X,22A,03X,DDDD.DDD";"Average frequency, Hz:",Fave 

9680 PRINT USING "3X,22A,03X,DDDD.DDD";"Maximum frequency, 

Hz:",Fmax 
9690 PRINT 

9700 Ll$="Percent change in shear modulus per degree C: " 

9710 L2$="Percent change in Young's modulus per degree C:" 

9720 L3$="Two point max -min zq>proach:" 

9730 L4$="Multi-point least-squares approach:" 

9740 SELECT UPC$(Block3$) 

9750 CASE "TORSIONAL" 

9760 Gmin=Scale_gtor_min 

9770 Gmax=Scale_gtor_max 

9780 PRINT USING "3X,30A,ZZ,5X,D.DDD";"Minimum 

" & C4a$ ,Po wer ,Gmin 

9790 PRINT USING "3X,30A,ZZ,5X,D.DDD";"Maximum 

"&C4a$,Power,Gmax 

9800 PRINT USING "3/,#" ! Three more line feeds 

9810 PRINT Ll$ 

9820 PRINT 

9830 PRINT USING "3X,27A,8X,DDD.DDD";L3$,Prcnt_per_c_2pt 

9840 PRINT USING "3X,35A,DDD.DDD";L4$,Prcnt_per_c_ls 

9850 CASE "FLEXURAL" 

9860 Emin=Scale_eflex_min 

9870 Emax=Scale_eflex_max 

9880 PRINT USING "3X,30A,ZZ,5X,D.DDD";"Minimum 

" & C4b$ ,Po wer , Emin 

9890 PRINT USING "3X,30A,ZZ,5X,D.DDD";"Maximum 

"&C4b$,Power,Emax 

9900 PRINT USING "3/,#" ! Three more line feeds 

9910 PRINT L2$ 

9920 PRINT 

9930 PRINT USING "3X,27A,8X,DDD.DDD";L3$,Prcnt_per_c_2pt 

9940 PRINT USING "3X,35A,DDD.DDD";L4$,Prcnt_per_c_ls 

9950 CASE "LONGITUDINAL" 

9960 Emin=Scale_elong_min 

9970 Emax=Scale_elong_max 

9980 PRINT USING "3X,32A,ZZ,5X,D.DDD";"Minimum 

"«feC4b$ ,Power,Emin 

9990 PRINT USING "3X,32A,ZZ,5X,D.DDD";"Maximum 

"&C4b$,Power,Emax 

10000 PRINT USING "3/,#" ! Three more line feeds 

10010 PRINT L2$ 

10020 PRINT 

10030 PRINT USING "3X,27A,8X,DDD.DDD";L3$,Prcnt_per_c_2pt 
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10040 PRINT USING "3X,35A,DDD.DDD";L4$,Prcnt_per_c_ls 

10050 END SELECT 

10060 PRINT Formfeeds ! Advance to top of next page 

10070 PRINTER IS CRT 

10080 RETURN 

10090 ! ****************************************************************** 
10100 Freq_time:! Produces frequency vs time graph w/curve 

10120 GOSUB Draw_freq_time 

10130 CALL Generic_curve(Pen3,Time(*),Resfreq(*),20*Argl) 

10140 RETURN 

10160 Output_device: ! Permits user to route graphs to the screen or plotter 

10180 OUTPUT KBD;Clear$; ! Clear the CRT 

10190 OUTPUT KBD;Home$; ! Home display 

10200 GRAPHICS OFF ! Turn off the graphics display 

10210 SELECT Plot_device$ 

10220 CASE <>"Plotter" 

10230 INPUT "Graphs appear on the screen, OK? (Y/N) ", Responses 

10240 IF UPCS(ResponseS)o"Y" THEN 

10250 Plot_deviceS="Plotter" 

10260 DISP "Graphs wUl be sent to the Plotter" 

10270 ELSE 

10280 DISP "Graphs will remain on the screen" 

10290 WATT 1 

10300 DISP 

10310 GRAPHICS ON 

10320 END IF 

10330 CASE "Plotter" 

10340 INPUT "Graphs are sent to the plotter, OK? (Y/N) ", Responses 

10350 IF UPCS(ResponseS)o"Y" THEN 

10360 Plot_device$="Screen" 

10370 DISP "Plots wUl be sent to the Screen" 

10380 ELSE 

10390 DISP "Plots wUl stay routed to the plotter" 

10400 END IF 

10410 END SELECT 

10420 WATT 1 

10430 DISP 

10440 RETURN 

Q^^QI ^3{<9k3k3i(3t(:k3i(3t(3f(3k3k3t(3t(9k4(3|c:tc:)(:t(3te3k3k3k3i(9{(3i(9{(9tC9{(9k9i<9{(3k3k3i(9{(^9i<9{(3k9i<3{(:i(3f(^:f(3f(3k3((3k3k3k3f(3f(3k3((9i(3f(:)(:)(:f0i(9i<^:(( 

10460 Xscale_option: ! Permits user to auto scale or manual scale the X axis 

J Q I % ak * a|f :(« :|c 3|< :|c 5|c % a|f 3(c % 5|e :(c a|c a(c 3{< 5|c 3|< a(e a|c 3(e 3{< 3|< :(« a|c 3|< :{< a|f % j|< 3j< jR 3|< 3{< jK 3j< 5|t 3|< jjt a|c 3|< a|c 3|c :fc 5|< 5k * ♦ ♦ 3k * 3|c ♦ 3k 5k * ak 3fe * * ♦ >k 

10480 OUTPUT KBD;ClearS; ! Clear the CRT 

10490 OUTPUT KBD;Home$; ! Home display 

10500 GRAPHICS OFF ! Turn off the graphics display 

10510 SELECT XscaleS 

10520 CASE "Auto scale X" 

10530 INPUT "X axis is automatically scaled, OK? (Y/N) ",ResponseS 

10540 IF ResponseSo"Y" THEN 
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Xsv'alr$="l'scT sc.ilc X" 

INri'T "Minimum X axis value ?"3onin_manual 
INrU'V "Maximum X axis v;ilue ?"J'Ciua\_manual 
IXTUV "Major X axis increment ?"JXmajor_manual 
INl'in' "Minor X axis increment ?".Xniinor_manual 
DlSr "'I'lie X axis w ill Ix’ sc;iled numuallv" 
l'L5U 

Dl5r "nie X axis will remain auto scaled" 

W AIT I 

Disr 

GRAOHlCS ON 
END IF 

CASE "I'ser scUe X" 

INTIT "X axis is m.uiually sc;iled. OK? ".Responses 
IF ResponseSo"Y" THEN 
XscaleS="Auto sc;ile X" 

DlSr "The X axis will be scaled automaticallv" 

ELSE 

INTlT "Ch.mce nunu.il limits ? (Y/N) ".Responses 
IF Respv>nseS="Y" THEN 
INTlT "Minimum X axis value ?",Xmin_manual 
INFIT "M.iximum X axis value ?".Xmax_manual 
INTlT "MajorX axis increment ?".Xmajor_manual 
INTlT "Minor X axis increment ?".Xminor_manual 
DISP TTe X axis will be scaled with new values" 

ELSE 

DlSr The X axis will remain m.inuallv scaled" 

A' AFT I 
D15F 

GR.AFHICS ON 
END IF 
END IF 
END SELECT 
W.ATT I 
DISF 
RETL’RN 



ile_occion; ! Permits user m auto sc.ile or m.uuKil sc.ile the Y axis 



OITFIT KBD.OearS. » Oe.ir the CRT 

OITFIT KBD'.HomeS. ! Home display 

GR-AFHICS OFF ! Turn otT the cniphics displav 

SELECT YscaleS 
CASE ”.Auto scale 

INTLT T’ axis is automaticallv sc;il«\i. OK ' C N"' ”.Res^vn.<e 
IF Res;v>aseS^'-“Y THEN 
YscaieS= I'ser .<ca!e Y 

INTlT 'Mirar.yam Y axis % alue .'".Vnun manual 
INTlT Maximum Y axis ' .due ' .Ym.x\ nnmual 
INTlT 'MaKT Y axis uvrement ' .Ynuior_ni.uiiial 
INTlT 'Maiv't Y axis ux'remeiit '“.Ynunor m.uni.al 



$ 



1 1060 DISP "The Y axis will be scaled manually" 

11070 ELSE 

1 1080 DISP "The Y axis will remain auto scaled" 

11090 WATTl 

11100 DISP 

11110 GRAPHICS ON 

11120 END IF 

11130 CASE "User scale Y" 

1 1 140 INPUT "Y axis is manually scaled, OK? (Y/N) ",Response$ 

11150 IF Response$<>"Y" THEN 

11160 Yscale$=" Auto scale Y" 

11170 DISP "The Y axis will be scaled automatically" 

11180 ELSE 

11190 INPUT "Change manual limits ? (Y/N) ",Response$ 

11200 IF Response$="Y" THEN 

11210 INPUT "Minimum Y axis value ?",Yinin_manual 

1 1220 INPUT "Maximum Y axis value ?",Ym;ix_mjuiual 

1 1230 INPUT "Major Y axis increment ?",Ymajor_manual 

1 1240 INPUT "Minor Y axis increment ?",Yminor_manual 

1 1250 DISP "The Y axis will be scaled witli new values" 

11260 ELSE 

1 1270 DISP "The Y axis will remain mmuially scaled" 

11280 WAIT 1 

11290 DISP 

1 1 300 GRAPHICS ON 

11310 END IF 

11320 END IF 

1 1330 END SELECT 
11340 WAITl 
11350 DISP 
11360 RETURN 

1 ] 3 TO 



1 1380 Temp_lime: ! Produces Temperature vs Time graph w/curve 
I 1 ^90! 

1 1400 Main_title$=Labcl$(7) 

1 1410 Sub_title$=Labcl$(8) 

1 1420 X_axis_namc$=Labcl$(3) 

1 1430 Y_axis_name$=Labcl$(5) 

11440 Xmin=0 

11450 Xmax=20 

1 1460 Ymin=Actual_lenip_min 

1 1470 Ymax=Actual_lcmp_max 

1 1480 GOSUB Gcncric_plol 

1 1490 CALL Gcncric_curvc(Pcn3,Timc(*),Tlicrmislor(’'‘),20*Argl ) 

11500 RETURN 



Minimum time is 0 hours 
I Maximum lime is 20 liours 

I Actual minimum tcm|)craturc 
! Actual maximum tcm|x;raturc 



1 1520 Frcq, temp: ! Produces Fret|uency versus leiii|>eiature graph w/cuive 

11540 Main_title$=Label$(6) 

11550 Sub_title$="" 

1 1560 X_axis_name$=Label$(5) 
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11570 Y_axis_name$=Label$(4) 
1 1 580 Xmin=Actual_ternp_min 
11590 Xmax=Actual_temp_max 
11600 Ymin=Actual_freq_min 
11610 Y max= Actual_freq_max 



! Actual minimum temperature 
! Actual maximum temperature 
! Actual minimum frequency 
! Actual maximum frequency 



11620 GOSUB Generic j>lot 

1 1630 CALL Generic_curve(Pen3,Thermistor(*),Resfreq(*),20*Argl) 

11640! PENPen4 

1 1 650! MOVE Actual_temp_min,Actual_temp_min*Slope+Intercept 

11660! DRAW Actual_temp_max^ctual_temp_max*Slope+Intercept 

11670! PENO 

11680 RETURN 



1 159Q !*♦*♦* ************ ************* *****’** ********************* ♦****♦*♦* 



1 1700 Create_new_file: ! Creates a file for a new data run 

11720 Pass$="Initial" 

11730 Grid_type$=Default_grid$ 

1 1 740 Xscale$="User scale X" 

1 1750 Yscale$="User scale Y" 

11760 Plot_device$="Screen" 

11770 BEEP 1000,.! 

1 1780 ON ERROR GOTO Fix2 

1 1790 INPUT "Filename to store the data ?",Filename$ 

1 1800 CREATE BOAT Filenames, 50 

1 1810 ASSIGN (o)Path_l TO Filenames 

11820 OFF ERROR 

11830 OUTPUT (o)Path_l;BlocklS,BIock2S,Block3S,Block4S 

11840 OUTPUT (2)Path_l;Mass,Length,Diameter, Density 

11850 BEEP 1000,.! 

1 1860 DISP "Output will be stored under Filename: ";FilenameS 

1 1870 WATT 1 

11880 DISP 

11890 IFnag=lTHEN 

1 1 900 FOR 1=0 TO 20* Arg 1 

11910 OUTPUT (S>Path_l;Mod(I),Thermistor(I),Resfreq(I),Inphmag(I), 

Redfreq(I) 

11920 NEXT I 

11930 END IF 

11940 RETURN 



j^j^^^Q|3|(3(c:|C3fC3f(:)C3i(3f(:(0f(3((3{C3|(3{(3|0f(^3fC3((a{(3((3fC3|(3{c%3{0|(3((^3|0|(3|0|(a|0fC3f(3f<3|(3fC3|(3{ca|(3{C3|<3{(3(<%3|(3((3|(3f(3|(3{c:|(3(0|(%3|(3(c3|(:((3|(:)c3|(%3i( 

11960 Open_old_file: ! Opens and retrieves data from an existing file 



11980 PassS=’’Follow on" 

1 1 990 Grid_typeS=Default_gridS 

12000 XscaleS-' Auto scale X" 

12010 Yscale$=" Auto scale Y" 

12020 Plot_deviceS=’’Screen" 

12030 BEEP 1000,.! 

1 2040 ON ERROR GOTO Fix 1 

12050 INPUT "Filename to retrieve the data ?", Filenames 

1 2060 ASSIGN (o)Path_l TO FUenameS 



no 



12070 ENTER @Path_l;Blockl$,Block2$,Block3$,Block4$ 

12080 ENTER (3)Path_l;Mass, Length, Diameter, Density 

12090 BEEP 1000,. 1 

12100 DISP "Retrieving data stored under Filename: ";Filename$ 

12110 FOR 1=0 TO 20* Arg 1 

12120 ENTER (o)Path_l ;Mod(I),Thermistor(I),Resfreq(I),Inphmag(l),Red 

freq(I) 

12130 Time(I)=I/Argl 

12140 NEXT I 

12150 OFF ERROR 
12160 ASSIGN (o)Path_l TO * 

12170 RETURN 



12180 !****************************************************************** 



12190 Fixl: ! 

122001 ****************************************************************** 
12210 SELECT ERRN 
12220 CASE 53 

12230 DISP "Limit file names to 10 characters. No punctuation." 

12240 CASE 56 

12250 DISP "This file doesn't exist on the data disk" 

12260 CASE 58 

12270 DISP "This file is not a BDAT file" 

12280 CASE ELSE 

12290 DISP " This file can not be processed by RT20B" 

12300 END SELECT 

12310 PRINTER IS CRT 

12320 BEEP 300,. 1 

12330 WATT 1 

12340 DISP 

12350 WATT 1 

12360 GOTO Open_old_fiIe 

1 23701****************************************************************** 



12380 Fix2: ! 



1 23901****************************************************************** 



12400 SELECT ERRN 
12410 CASE 53 

12420 DISP "Limit file names to 10 characters. No punctuation." 

12430 CASE 54 

12440 DISP "Duplicate filename! Try another name." 

12450 END SELECT 

12460 BEEP 300,. 1 

12470 WATT 1 

12480 DISP 

12490 WATT 1 

12500 GOTO Create_new_file 

1 75 1 nf ****************************************************************** 



12520 Close_file: I Closes new data file after data collection is completed 

12540 GCLEAR 

12550 ASSIGN (SPath.l TO * 

12560 BEEP 1000,. 1 
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12570 DISP "Data is stored under Filename: ";Filename$ 

12580 WATT 1 
12590 DISP 

12600 RETURN 

1 2610 !***************************************************************’*'** 
12620 Program_end: ! Shuts down shop, plays a little melody 

12630!****************************************************************** 
12640 OUTPUT KBD;Clear$; 

12650 GRAPHICS OFF 

12660 CONTROL 1,12;0 

12670 BEEP 157,.l 

12680 BEEP201,.! 

12690 BEEP 178,-1 

12700 BEEP 272,-1 

12710 WATT -5 

12720 BEEP 272,-1 

12730 BEEP 178,-1 

12740 BEEP 157,-1 

12750 BEEP 201,-1 

12760 DISP 

12770 DISP "Press RUN when you are ready for another try 

12780 END 

1 27 90 

12800 SUB Label(Csize,Asp_ratio,Ldir,Lorg,Pen,Xpos,Ypos,Text$) 

^ ^ ^ 3|c sfc 3|c sfc 9|c 3|c ?|f 3|c 3|c 3|c ?| f 3|c 3|C ?| f 3|c ?| f 3|c 3|c 3|C 3|C 3|C 3|( S|C 3|c 3|C ?| t 3|c 3|c 3|( ?| f 3|c 3|c 3|c 3|C S|C 3|C 3|c Sfc 3|c 3|c ?| f 3|c 3|c 9|C 3|^ 9|C 9|C 9|C Sfc ?|f 3|c 9|C 3|C ?| f 3|c 9|c 3|C 3|c 3|C ?|f 3|c ?| f s|c 

12820! This subroutine defines several systems variables (Csize, LDIR, etc-), 

1 2830! and labels the text (if any) accordingly- 
12840 DEG 

12850 CSIZE Csize ,Asp_ratio 

12860 LDIRLdir 

12870 LORG Lorg 

12880 PEN Pen 

12890 MOVE Xpos,Ypos 

12900 IF Text$o"" THEN LABEL USING "#,K";Text$ 

12910 PENUP 
12920 SUB END 

1 2930 ! ****************************************************************** 
12940 SUB Lbl_axes(Csize,Asp_ratio,Pen,Xmin,Xmax,Xstep,Ymin,Ymax,Ystep) 

12950 ! ****************************************************************** 
12960 DEG 

12970 CSIZE Csize,Asp_ratio 

12980 PEN Pen 

12990 CLIP OFF 

13000 LDIR 0 

13010 LORG 6 

13020 Yrange=Ymax-Ymin 

13030 Yoffset=-02*Yrange 

13040 FOR L=Xmin TO Xmax STEP Xstep 

13050 MOVE L,Y min- Y offset 

13060 IF ABS(L)<-001 THEN L=0 

1 3070 LABEL USING "#,K";L 
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n080 NEXTL 
13090 LORG8 

13100 Xoffset=.02*(Xmax-Xmin) 

13110 IFYrange<=l THEN 

13120 Mmax=DROUND(YrangeA'stq3,l) 

13130 Yval=Ymin 

1 3 140 FOR M=0 TO Mmax 

13150 IF ABS(Yval)<=.001 THEN Yval=0 

13160 MOVE Xmin-Xoffset,Y val 

13170 LABEL USING "#,SD.DD’’;Yval 

13180 Yval=Yval+Ystep 

13190 NEXTM 

13200 ELSE 

13210 FOR M=Ymin TO Ymax STEP Ystep 

13220 IF ABS(M)<=.001 THEN M=0 

13230 MOVE Xmin-Xoffset,M 

13240 LABEL USING "#.K";M 

13250 NEXT M 

13260 END IF 

13270 CUP ON 

13280 PENUP 

13290 SUBEND 

13310 SUB Block_info(Blockl$,Block2$,Block3$,BIock4$,Xgdumax,Pen,Pen2) 

13330 CALL LabeI(3,.6,0,2,Pen,2,2,"Rod: ”&Blockl$) 

13340 CALL Label(3,.6,0,5,Pen,Xgdumax/3,2,"Run: "«&Block2$) 

13350 CALL Label(3,.6,0,5,Pen,Xgdumax*2/3,2,"Mode: ”&Block3$) 

13360 CALL Label(3,.6,0,8.Pen,.97*Xgdumax,2,Block4$) 

13370 MOVE 0,4 
13380 PEN Pen2 
13390 DRAW 133,4 
13400 PENUP 
13410 SUBEND 

j 2420 !********************************************************** ******** 
13430 SUB Least_squares(Imax,X(*),Y(*),Slp,Int,Cor,Slp_er,Int_er,Xmean,Ymean) 
2440 ! ****************************************************************** 

13450 DISP ’’Computing least-squares fit " 

13460 Sumx=0 

13470 Sumy=0 

13480 Sumxx=0 

13490 Sumxy=0 

13500 FOR 1=0 TO Imax 

13510 Sumx=Sumx+X(I) 

13520 Sumxx=Sumxx+X(I)^2 

13530 Sumy=Sumy+Y(I) 

1 3 540 Sumxy=Sumxy+X(I) * Y (I) 

13550 NEXT I 

13560 Delta=(Imax+l)*Sumxx-Sumx^2 
13570 Int=(Sumxx *Sumy-Sumx* Sumxy )/Delta 
13580 Slp=((Imax+ 1 )*Sumxy-Sumx*Sumy )/Delta 
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13590 Sumeirerr=0 
13600 FOR J=0 TO Imax 

13610 Sumerrerr=Sumerrerr+(Y(J)-Int-Slp*X(J))*2 
13620 NEXT! 

13630 Sigmayy=Sumerrerr/(Imax+l-2) 

13640 Sigmay=SQR(Sigmayy) 

13650 Int_er=Sigmay*SQR(Sumxx/Delta) 

1 3660 Slp_er=Sigmay*SQR((Imax+ 1 )/Delta) 

13670 Xbarsum=0 
13680 Ybarsum=0 
1 3690 FOR K=0 TO Imax 
13700 Xbarsum=Xbarsum+X(K) 

13710 Ybarsum=Ybarsum+Y(K) 

13720 NEXTK 

13730 Xmean=Xbarsum/(Imax+l) 

13740 Ymean=Ybarsum/(Imax+l) 

13750 Sigmaxxsum=0 
13760 Sigmayysum=0 
13770 Sigmaxysum=0 
1 3780 FOR L=0 TO Imax 

1 3790 Sigmaxxsum=Sigmaxxsum+(X(L)-Xmean)^2 

1 3800 Sigmayysum=Sigmayysum+( Y (L)-Y mean)''2 

13810 Sigmaxysum=Sigmaxysiim+(X(L)-Xmean)*(Y(L)-Ymean) 

13820 NEXTL 

13830 Sigmax=SQR(Sigmaxxsum/(Imax+l)) 

1 3840 Sigmay l=SQR(Sigmayysum/(Imax+ 1 )) 

13850 Sigmaxy=Sigmaxysum/(Imax+l) 

1 3860 Cor=Sigmaxy/(Sigmax*Sigmay 1 ) 

13870 DISP 
13880 SUBEND 

1 3890 !****************************************************************** 

13900 SUB Scale(L,R,B,T,Top,Left^center,Ycenter,Xgdumax,Ygdumax) 

13910 ******************************************************************* 

13920 Top=T*Ygdumax 

13930 Bottom=B*Ygdumax 

1 3940 Left=L*Xgdumax 

13950 Right=R*Xgdumax 

13960 Xcenter=(Right+Left)/2 

13970 Ycenter=(Top+Bottom)/2 

13980 VIEWPORT Left, Right3ottom, Top 

13990 SUBEND 

1 4000 1************^^^*^*****^^*^**^**^^^*^***^*^**^*******^***^^^^*^*^** 
14010 SUB Yscale(Ymin,Ymax,Y minor, Ymajor) 

I 2Q j 3k 3k 3k 3k 3k 3k 3k ak 3k 3k 3k 3k 3k * 3k ak )k ak 3k )k ak 3k 3(c 3k % 3k 3k 3f( iK 3k 3k ak 3k )k ak 3k ak 3k 3k 3k 3k 3k 3k 3k ak 3k ak 3k 3k 3k 3k 3|c 3k a|c ^ ak 

14030 DIM Diff(36)>Imor(36) 

14040 DATA .1,.005,.2,.01,.25,.01,.3,.02,.4,.02,.5,.02,.75,.05,1 ,.05 

14050 DATA 5,.2,10,.5,15, 1,20,1 ,25,1,30,2,40,2,50,2,75,5,100,5,125,5 

14060 DATA 150,10,200,10,250,10,300,20,400,20,500,20,750,50,1000,50 

14070 DATA 1250,50,1500,100,2000,100,2500,100,3000,200,4000,200,5000 

14080 DATA 200,7500,500,10000,500 
14090 FOR 1=1 TO 36 
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14100 READ Diff(I),Minor(I) 

14110 NEXTI 

14120 Yrange=Ymax-Ymin 

14130 Index=l 

14140 WHILE Y range>Diff(Index) 

14150 Index=Index+l 

14160 END WHILE 
14170 Yminor=Minor(Index) 

14180 Ymajor=Diff(Index)/5 
14190 IFYmin<0THEN 

14200 Newmin=Yrnin-Ymajor+ABS(Yrnin MOD Ymajor) 

14210 ELSE 

1 4220 Newmin= Y min-( Y min MOD Y major) 

14230 END IF 

14240 Newmax=Newmin+Diff(Index) 

14250 WHILE Ymax>Newmax 
14260 Index=Index+l 

14270 Yminor=Minor(Index) 

14280 Y major=Diff(Index)/5 

14290 IF Ymin<0 THEN 

14300 Newmin=Ymin-Y major+ABS(Ymin MOD Y major) 

14310 ELSE 

14320 Newmin=Ymin-(Ymin MOD Ymajor) 

14330 END IF 

14340 Newmax=Newmin+Diff(Index) 

14350 END WHILE 
14360 Ymax=Newmax 
14370 Ymin=Newmin 
14380 SUBEND 

14400 SUB Xscale(Xmin,Xmax,Xminor,Xmajor) 

Q I :k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 9k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 3k 

14420 DIM Diffe(40),Minorx(40),Majorx(40) 

14430 DATA 5. ,0.2, 1.0, 6.,. 2, L.7.0, 0.2, 1.0, 8.0, 0.2, 1.0, 10.0, 0.5, 2.0, 12.0, 0.5 

14440 DATA 2.0,14.0,0.5,2.0,16.0,0.5,2.0,20.0,1.0,4.0,25.0,1.0,5.0,30.0,1.0 

14450 DATA 5.0,35.0,1. 0,5.0,40.0,1. 0,5.0,50.0,2.0,10.0,60.,2.0,10.0,70.0,2.0 

14460 DATA 10.0,80.0,2.0,10.0,100.0,5.0,20.0,120.0,5.0,20.0,140.0,5.0,20.0 

14470 DATA 200.0,1 0.0,20.0,400.0, 1 0.0,50.0,600.0,20.0, 1 00.0,800.0,50.0, 100.0 

14480 DATA 1000.0,50.0,100.0,1200.0,50.0,200.0,1400.0,50.0,200.0,1600.0 

14490 DATA 50.0,200.0,1800.0,100.0,300.0,2000.0,100.0,500.0,2500.,100.,500.0 

14500 DATA 3000.,100.,500.0,3500.,100.,500.0,4000.,200.,400.0 

14510 FOR 1=1 TO 34 

14520 READ Diffx(I),Minorx(I),Majorx(I) 

14530 NEXT I 

14540 Xrange=Xmax-Xmin 

14550 Index=l 

14560 WHILE Xrange>Diffe(Index) 

14570 Index=Index+l 

14580 END WHILE 
14590 Xminor=Minorx (Index) 

14600 Xmajor=Majorx(Index) 
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14610 IF Xmin<0 THEN 

14620 Newmin=Xmin-Xmajor+ABS(Xmin MOD Xmajor) 

14630 ELSE 

14640 Newmin=Xrnin-(Xmin MOD Xmajor) 

14650 END IF 

14660 Newmax=Newmin+Diffe(Index) 

14670 WHILE Xmax>Newmax 
14680 Index=Index+l 

14690 Xminor=Minorx(Index) 

14700 XmaJor=Majorx(Index) 

14710 IFXmin<0THEN 

14720 Newmin=Xmin-Xmajor+ABS(Xmin MOD Xmajor) 

14730 ELSE 

14740 Newrnin=Xmin-(Xmin MOD Xmajor) 

14750 END IF 

14760 Newmax=Newmin+Diffx(Index) 

14770 END WHILE 
14780 Xmax=Newmax 
14790 Xmin=Newmin 
14800 SUBEND 

148 10 ! ****************************************************** ************ 
14820 SUB Generic_curve(Pen^(*),Y(*),Max j)oint) 

1 *♦♦♦*** ****************************** 3k* **♦*♦♦** 4c* *:(< *:!<:*<******* **:|<*:|< 

14840 PEN Pen 

14850 FOR Point=0 TO Max_point 
14860 PLOT X(Point),Y(Point) 

14870 NEXT Point 
14880 PENUP 
14890 SUBEND 

14910! Nothing follows. 

14920!****************************************************************** 

14930!****************************************************************** 

(4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<4<^4<^4<4<4<4<4<4<4<4<4<^3|c3|(3ic3k3ic3k3(c3k3k3k3k^3k3k3ic3f(3|c3k 

14950 SUB Calc_redfreq(Imax,Thermistor(*),Resfreq(*),Redfreq(*)) 

14970 FORI=OTOImax 

14980 Temp(I)=Thermistor(I)+273.15 

1 4990 X(I)=- 1 29*(Temp(I)-283 . 1 5)/( 1 07+Temp(I)-28 3.15) 

15000 Y(I)=Xa)*LOG(10) 

15010 Alphat(I)=EXP(Ya)) 

15020 Redfreqa)=Resffeq(I)*Alphata) 

15030 NEXT I 
15040 SUB END 
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APPENDIX E. LISTING OF HP BOAT FILE TO MACINTOSH 
TEXT FILE CONVERSION PROGRAM 
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10 PRINTER IS CRT; WIDTH 80 

20 PRj wnr 

30 PRINT "^PROGRAM HPTOMAC: Program to transfer BOAT files on the HP to" 

40 PRINT "TEXT files readable by Cricket Graph on the Macintosh." 

50 PRINT "Ref. pp. C-6 and C-17 of the Cricket Graph User’s Manual." 

60 PRINT "Author: Steve Baker Last revision date: 10 Aug 1989" 

70 PRINT '**HARDWARE REQUIRED: HP Series 200 with 98644A serial interface 



80 PRINT "Series 300 with built-in interface, Macintosh with built-in serial" 

90 PRINT "interface, Hayes modem cable for the Macintosh (DB-9 or Din-8 to" 

100 PRINT "male DB-25 connectors. Ref. p.G-5 of the VersaTerm Pro User’s" 

1 10 PRINT ’’Manual), standard male-male RS-232 cable with at least pins 1-8" 

120 PRINT "and 20 wired straight through for the HP. These cables can each" 

130 PRINT "be hooked up directly to a modem, or they can be connected to each" 

140 PRINT "other through a null modem with the following pins connected:" 

150 PRINT "Itol, 2to3, 3to2, 4and5to8, 6to20, 7to7, 8to4and5, 20to6." 

160 PRINT ’’^SOFTWARE REQUIRED: HPTOMAC running on an HP Series 200/300 
and a" 



170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

and" 

310 

320 

330 

340 

350 

360 



PRINT "communications program capable of faithfully capturing incoming" 

PRINT "text (ASCII characters) to a TEXT-type file running on the Mac," 

PRINT ’’such as MacTerm (desk accessory bundled with Borland SideKick) or" 
PRINT "VersaTerm PRO. NOTE: text captured by Microsoft Works cannot be" 
PRINT "read by Cricket Graph." 

PRINT "*INSTRUCTIONS FOR USE:’’ 

PRINT "Both serial ports must be set to 2400 baud, 8 data bits, one stop" 

PRINT "bit, no parity. The following will correctly set up both machines" 

PRINT "following power-up." 

PRINT ’’*SETTING UP THE HP: HPTOMAC asks for the select code of the" 
PRINT "serial port and sets the baud rate to 2400 baud. The other" 

PRINT "settings are automatic at startup, so no additional user action is" 

PRINT "required (Ref. p. 13-10 of Vol.2 of BASIC Interfacing Techniques)." 
PRINT ’’*SETTING UP MACTERM: launch the ’’’’Configure MacTerm’”’ utUity 

PRINT "click on the appropriate buttons to set serial port parameters" 

PRINT "(Ref. p. 107 of SideKick v2.0 User’s Manual). Clicking on "’’Save" 

PRINT ’’Setup’”’ saves these settings. The serial port parameters may also" 

PRINT "be set in the ’’’’MacTerm’"’ menu while mnning MacTerm (MacTerm is" 
PRINT "launched from the desk accessory menu)." 

PRINT ’’*SETTING UP VERSATERM PRO: Launch VersaTerm PRO. Set the 



baud" 



370 PRINT "rate to 2400 in the ’’’’Baud’’’’ menu. Enable Xon/Xoff, set Parity to" 

380 PRINT "none. Char Size to 8 bits, and Stop Bits to 1.0 in the" 

390 PRINT ’’’’’’Settings’”’ menu. Otlier features may be enabled and disabled in" 

400 PRINT "the ’’’’Extras’”’ dialog box, selected from the ’’’’Settings’’’’ menu" 

410 PRINT "(Ref. pp. F-15 to F-18 of the User’s Manual). Be sure that ""Auto" 

420 PRINT "Tek 4014 Entry’”’ is disabled." 

430 PRINT "*TRANSFERRING A BDAT FILE from the HP to the Mac is 



accomplished" 

440 PRINT "in text recording mode under MacTerm (ref. p. 109 in User’s" 

450 PRINT ’’Manual) and in Save Stream mode under VersaTerm PRO (ref. pp." 
460 PRINT’’B-28,F-8,F-13, F-15 in the User’s Manual). Refer to each User’s" 
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470 PRINT "Manual for details. HPTOMAC is self-prompting." 

480 PRINT "*NOTE that PROG files may also be transferred over the serial port" 

490 PRINT "connection simply by LOADing each PROG file into the HP's memory" 

500 PRINT "and issuing a LIST #Sc command, where Sc is the select code of the" 

510 PRINT "serial port. The setup and procedure for receiving a PROG file" 

520 PRINT "under MacTerm or VersaTerm PRO is identical to that for receiving" 

530 PRINT "a BOAT fUe." 

550 ! 

560 !*****PROGRAM BEGINS HERE************************************* 
570 ! 

580 DIMCh$[l]J^ame$(100)[31] 

581 DIM C1$[50],C2$[50],C3$[50],C4$[50] 

590 INTEGER OutdevJ,N 

600 PRINT "PROGRAM HPTOMAC: Program to transfer BOAT files on the HP to 
TEXT" 

610 PRINT "files readable by Cricket Graph on the Macintosh. Detailed" 

620 PRINT "instructions for using this program are in the comments above." 

630 BEEP 
640 PRINT 

650 PRINT MAKE SURE THE MAC SERIAL PORT IS SET TO 2400 BAUD, 

8 DATA" 

660 PRINT "BITS, 1 STOP BIT, NO PARITY *****'' 

670 PRINT 

680 MASS STORAGE IS ":,4,0" 

690 PRINT "Mass storage default has been set to :,4,0" 

700 INPUT "Enter select code of output device (CRT=1 , PRT=701 , Serial Port=9 o 
r 10):",Outdev 

7 10 IF Outdev=9 OR Outdev=10 THEN CONTROL Outdev,3;2400 
720 INPUT "Enter name of source data fiIe:",Srcfile$ 

730 ASSIGN @Pathl TO Srcfile$ 

740 INPUT "Enter number of elements per record (1 to 100):",N 
750 IF N<1 OR N>100 THEN GOTO 740 

760 INPUT "Do you want to send element names? (Y or N, default is N)",Ch$ 

770 IF Ch$="Y" OR Ch$="y" THEN 
780 Ch$="Y" 

790 FOR 1=1 TON 

800 PRINT "Enter name of ";I;"th element:" 

810 LINPUT Name$(I) 

820 NEXT I 
830 END IF 

840 INPUT "Hit CONTINUE to continue",Ch$ 

850 PRINT 

860 PRINT "Sending output to device select code ";Outdev 

870 PRINTER IS Outdev;WIDTH OFF 

880 IF Ch$="Y" THEN 

890 PRINT 

900 FOR 1=1 TO N-1 

910 PRINT Name$(I);CHR$(9); 

920 NEXT I 

930 PRINT Name$(N) 
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940 END IF 

950 ON END @Pathl GOTO 1040 

951 ENTER @Pathl;Cl$,C2$,C3$,C4$ 

952 ENTER @Pathl;Mass,Length, Diameter, Density 
960 REPEAT 

970 FOR 1=1 TON-1 
980 ENTER (3) Path 1; Data 

990 PRINT Data;CHR$(9); 

1000 NEXT I 
1010 ENTER (o)Path 1 ;Data 
1020 PRINT Data 
1030 UNTIL False 
1040 OFF END (©Pathl 
1050 ASSIGN (SPathl TO * 

1060 BEEP 

1070 PRINTER IS CRT; WIDTH 80 

1080 PRINT "Transfer complete. Close Mac file." 

1090 PRINT "Open new Mac file and enter HP source file name to transfer" 

1 100 INPUT "another file with the same attributes or hit RUN to restart: ",Srcfi 
le$ 

1 1 10 ASSIGN @Pathl TO Srcfile$ 

1 120 GOTO 840 
1130 END 
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